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/''35TRACT 


The  transverse  modes  of  high  frequency  combustion  instability  in 
liquid  propellant  rocket  motors  are  studied  both  theoretically  and  experimentally. 
The  nature  of  the  phenomenon  is  described,  and  previous  c-xper imer, la I  and 
theroretical  work  on  the  problem  's  'iisrussed.  A  prel  o''T'oratory  series 

of  experiments  demonstrates  the  need  for  additional  theoretical  analysis  in  order 
to  explain  three  effects:  (I)  the  influence  on  stability  of  the  orientation 
of  The  eier„8i-its  of  the  injection  pnitern  with  respect  to  rne  radial  direction, 

(2)  the  greater  tendency  to  instability  ol  the  spinning  forms  of  tangential 
modes  as  compared  to  the  standing  forms.  (3)  the  effect  of  the  distribution  of 
injection  across  the  injector  face. 

The  basic  theory  is  derived  for  the  case  of  a  circularly  cylindrical 
combustion  chamber  with  uniformly  distributed  injection.  Follow  ng  Crocco, 
the  concept  of  a  sensitive  time  lag  is  introduced.  However,  dependence  of  the 
time  lag  on  the  radial  and  tangential  velocity  components  as  well  as  on  the 
thermodynamic  state  prooerties  (represented  by  the  pressure)  is  assumed.  S'lall 
perturbation  theory  is  -sed,  in  conjunction  with  an  approximate  order  of 
magnitude  analysis,  fo  derive  a  characteristic  equation  relating  +h5  'requsney 
of  neutral  oscillation  at  the  stability  limits  to  the  parameleis  describing 
The  combustion  process.  Extensions  of  the  basic  tl'e^ry  o-s  made  to  cover  the 
cases  of  (I)  nonuniformly  distribured  i.njection,  (2)  a  variable  angle  sector 
motor,  and  (3)  more  general  time  lag  formulations.  The  use  of  the  sector  motor 
in  the  study  of  tangential  modes  of  instability  is  discussed.  Calculations 
for  an  'idealized  rocket  motor  are  presented  in  order  to  illustrate  the  results 
of  the  theoretical  analysis. 

A  series  of  experiments  using  both  sector  and  full  cylindrical  mofors 
is  described.  It  is  shown  that  partial  and  qualitative  verification  is  obtained 


-  i  i  i  - 

for  the  theoretical  predictions  that  (I)  tl j  standing  form  of  a  fangential 
mode  IS  more  stable  than  the  spinning  form  because  the  tangential  velocity 
6-ftects  are  inoperative  in  the  former,  (2)  the  radial  velocity  f I uctuati o.iS 
cause  a  shifi  of  the  unstaole  ranges  of  the  sensitive  time  lag,  (3)  increasing 
•tability  is  obfained  as  +^>e  propellant  injection  is  moved  toward  the  center 
of  the  chamber,  (4)  the  pressure  dependenr  combustion  parameters  correspond  to 
those  v/hics  govern  the  longitudinal  modes  of  instabili+y.  Tne  validiiy  of  the 
linearized  approach  is  assessed  and  areas  for  further  theoretical  and  experimental 
work  are  indicator 
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CHnPTER  I 


Introducvion 


A  Description  of  the  Phenomenon 

In  a  liquid  propellant  rocket  motor,  the  combustion  process  is  never 
entirely  smooth.  During  the  steady  state  period  between  starting  ana  cutoff, 
fluctuations  occur  in  all  important  properties  (pressure,  temperature,  velocity, 
etc.)  around  the  desired  operating  values.  The  amplitude  of  such  1 I uctuat'on^ 
can  vary  over  a  wide  range,  from  motor  to  motor,  and  in  one  motor  for  different 
operating  conditions.  When  the  fluctuations  are  completely  random,  the 
operation  is  classified  as  "rough"  combustion.  A  qualitative  distinction  can 
be  made  between  "rough"  and  "smooth"  combustion  on  the  basis  of  the  root- 
mean-square  amplitude  of  the  fluctuations. 

Combustion  instability,  on  the  other  hand,  consists  of  organized 
oscillations,  which  are  maintained  and  amplified  by  the  combustion  process 
itself.  It  is  generally  agreed  that  there  are  three  kinas  of  i nstabi I i fy, 
differentiated  by  the  frequency  of  the  resulting  oscillations.  These  are 
referred  to  as  low,  intermediate  and  high  frequency  instability. 

Low  frequency  unstable  combustion  has  received  much  theorerical  ?nci 
experimental  study  ( I ,2,3,4)*.  It  is  now  rather  well  understood  and  compara¬ 
tively  easy  to  ovoid  or  cure,  A  theoretical  analysis  exists  for  the  interme¬ 
diate  frequency  case  (5),  but  the  phenomenon  is  rarely  observed.  This 
dissertation  is  concerned  with  the  high  frequency  type  of  combustion  insta¬ 
bility,  which  is  by  far  the  most  dest-uctive  type  and  the  most  difficult  to 
contrc I . 


*  Numbers  in  parentheses  refer  +o  References  listed  on  np.  [20,  ',2! 


High  frequency  ins+ability  is  a  forced  oscillafion  of  +he  combustion 
chamber  gases,  driven  by  the  combustion  process  interacting  with  the  resonan-e 
effects  of  the  chamber  geometry,  with  damping  supplied  by  the  exhaust  nozzle. 
The  resulting  oscillation  patterns  are  very  similar  to  those  of  the  acoustic 
modes  of  the  chamber.  For  a  cylindrical  chamber,  one  can  distinguish  between 
longitudinal  and  transverse  acoustic  modes.  In  a  longl+udinal  mode  the  prop¬ 
erty  variations  occur  in  the  axial  direction,  while  conoitions  are  uni  form  on 
a  section  normal  to  the  axis.  Purely  transverse  modes  involve  uniformity  of 
the  properties  along  the  axis,  but  variations  in  tne  radial  and  circumferential 
directions.  Figure  I  shows  the  acoustic  oscillation  patterns  for  several  purely 
transverse  modes  for  a  circular  cross  section.  The  most  commonly  observed 
modes,  in  order  of  increasing  frequency,  are  the  first  tangential,  second 
tangential,  first  radial,  and  the  first  combined  radial  tangential.  In  the 
"P"  diagrams,  several  typical  isobars  are  drawn,  for  a  given  instant,  with 
the  solid  lines  indicating  values  in  excess  of  the  mean  pressure  iuvel.  The 
"V"  diagrams  show  streamlines  at  the  same  instant.  Each  acoustic  mode  has  a 
characteristic  frequency,  which  is  related  to  the  S  number  v.-i+h  the 

diagrams  (the  significance  of  this  number  wilt  be  discussed  mce  fuliy  in  a 
later  section).  It  has  been  observed  that  the  frequencies  of  fuliv  develop¬ 
ed  pure  transverse  mode  combustion  oscillations  fall  within  few  percent  of 
Tne  acoustic  frequencies  for  the  corresponding  modes.  Each  tangential  or 
combined  mode  can  exist  in  two  forms:  the  standing  form,  in  which  the  n^dal 
sufaces  (indicated  by  the  solid  diametral  fines  of  Figure  I)  are  stationary, 
and  the  spinning  form,  in  which  the  nodal  surfaces  rotate  at  the  angular  fre¬ 
quency  of  the  oscillation.  The  two  forms  of  the  first  tangential  mode  ara 
illustrated  in  Figure  2  which  shows  schematically  pressure  and  velocity 
patterns  ot  four  instants  during  a  period  of  oscillation. 


Purely  lonqitudindl  modes  occur  in  chambers  v/ith  large  length  to 
diameier  rafio-  kLIQ  I);  purely  transverse  modes,  for  L/D«^  combined 
I ong i tudina I -transverse  modes  may  be  observed  when  L/D"'l  ;  provided  in  all 
cases  ihat  certain  requirements  for  instability  are  met.  The  derivation  and 
discussion  of  these  requiremer.+s  is  the  substance  of  this  dissertatior. 

A  linearized  theoretical  ana (vs is  will  be  employed;  although  such 
a  theory  is  capable  of  prediciing  tiie  onset  of  instability,  it  is  entirely 
unabie  to  determine  the  characteristics  of  the  fully  developed  oscillations. 
In  particular  the  amplitude  of  the  fluctuations,  which  is  set  Ly  various 
competing  non-linear  processes,  may  be  small  or^targe:,- depepdthgf'ion-the.  - 
mode  and  the  operating  conditions.  The  amplitudes  of  purely  transverse  modes 
have  been  observed  in  most  cases  to  be  very  large;  peak-to-peak  values  of 
up  to  300^  of  steady  state  chamber  pressure  were  measured  in  the  course  of 
this  investigation.  One  of  the  most  detrimental  consequences  of  such 
oscillations  is  the  burnout  of  combustion  chamber  hardware  because  of  the 
abnormally  high  heat  transfer  rates  associated  with  instability.  In  one 
ser'es  of  tests  (6)  a  jump  in  wall  heat  flux  from  3  up  to  15  Btu/ i n^-sec , • 
was  measured  accompanying  the  onset  of  a  transverse  mode  oscillation.  Such 
effects  provide  a  strong  motivation  for  investigations  leading  to  a  better 
understanding  and  control  of  this  phenomenon. 

B,  Previous  experimental  work 

A  considerable  amount  of  experimental  work  has  been  done  on 
combustion  instability  in  rocket  motors.  However,  very  little  of  this  effort 
has  been  of  a  systematic  nature,  or  has  been  connected  with  theore+ical 
Considerations.  The  initial  studies  of  Be-man  and  Cheney  (7)  of  longitudr 
inal  modes  of  oscillation  used  slit  v/indow  photographic  methods.  They 
detected  the  buildup  of  small  disturbances  into  high  ampTitude  shock  type 
waves,  and  showed  that  the  frequencies  of  oscillation  approximated  the 
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acous+ic  frequencies  within  the  comb'  stion  chamber.  FI  Ms  and  his  asso^’ates 

(8)  correlated  measuremesits  of  pressure  fluctuations  with  the  optical  methods 
used  almost  exclusively  by  Berman  and  Cheney  in  their  initial  studio's. 

Later  work  by  Ellis  extended  the  optical  studies  to  transverse  modes 

(9) .  These  investigations  showed  high  particle  velocities  (of  the  order  of 
2000  fps)  and  spiral  particle  trajec1cr'_s  (in  a  plane  norir.w.i  to  combus¬ 
tion  chamber  axis).  The  changes  in  the  oscillation  patterns  observed  in  using 
chambers  of  different  diameters  and  different  injector  configurations  were 
explained  i i iati vei y  by  means  of  a  ’’propellant  preparation  time"  and  its 
relation  to  the  "wave  travel  time",  the  fatter  being  a  function  of  the  chamber 
size. 

Various  methods  have  been  developed  by  the  industry  to  compare  the 
stability  of  different  motor  configurations  by  pulsing  during  steady  state 
rocket  operation.  In  one  method,  a  series  of  calibrated  powder  charges  is 
fired  tangentJaf  I'j^Mrttootne-chamhep  ih'order  to  disturb-'the-c6mbust}on’'and 
flow  patterns.  Each  charge  is  stronger  than  the  preceding  one.  Shutdown  of 
the  rocket  motor  is  effected  automatically  following  the  onset  of  instabili¬ 
ty.  The  number  of  charges  fired  is  thus  a  measure  of  the  relative  stability 
of  the  motor  at  the  particular  operating  conditions  of  the  test,  ’rfhiio  such 
studies  are  of  undoubted  value  in  the  development  of  a  rocket  engine  system, 
they  contribute  •--latively  little  to  the  findamental  '.jndeistanding  of  the 
phenomenon  of  instability. 

Male,  Kerslake,  and  Tischler  (10)  at  the  NASA  Lewis  Research  Center 
made  optica!  studies  which  corroborated  those  of  other  experimenters  (7,  8). 

In  addition,  they  noted  some  inters.rion  effects  between  longitudinal  and 
transverse  modes  of  oscillation,  as  well  as  the  greatly  increased  heat  transfer 
to  the  combustion  chamber  walls  which  is  associated  with  transverse  modes  of 
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instability.  Other  interesting  aspects  of  combustion  instability  involving 
chainber  geometry  were  investigated  by  this  group  (Iv.  However,  no  attempt 
was  made  to  determine  the  underlying  mechanisms. 

Investigations  by  Osborn  and  Bonne  I  I  (!2,  13)  using  a  gas  rocket 
system  have  sought  to  assess  the  effects  of  such  parameters  as  chamber 
geometry,  chamber  pressure,  and  nroce! lant  cnemistry  on  combustion  ■■'stability. 
Working  wi+b  a  constant  diameter,  variable  length  motor,  they  have  noted  with 
some  propellants  sharp  changes  in  The  stability  behavior  for  transverse  modes 
when  the  chamner  length  is  increased  to  the  point  where  longitudinal  mode 
instability  is  possible,  thus  indicating  important  interaction  effects  between 
longituoinal  and  transverse  modes.  They  have  also  reported  widened  instability 
regions  with  propellants  of  higher  heat  release  rates. 

All  of  the  investigations  described  above  were  restricted  to  the 
purely  empirical  approach,  rather  than  attempting  fundamental,  systematic 
studies  based  on  a  quantitative  theory. 

C.  Previous  theoretical  work 

Probably  the  first  theoretical  investigation  of  the  interaction  of 
pressure  waves  and  combustion  was  that  made  by  Rayleigh  in  his  explanation 
of  Rijke  tones  (II).  He  concluded  that  in  order  for  heat  to  drive  a  stanJing 
wave,  the  heat  input  should  maximize  with  time  when  the  pressure  is  at  its 
local  maximum.  Subsequent  investigations  have  extended  Rayleigh's  analysis 
To  more  general  combustion  systems. 

The  concept  of  a  combusfion  time  lag  as  a  coordinating  influence  to 
excite  organized  osciiiarions  in  a  liquid  propellant  rocket  combustion  chamber 
originated  in  von  Karman's  group  at  the  jet  Proplulsion  Laboratory  in  1941. 
Guoder  and  Friant  (2),  Yachter  (3)  and  Summerfield  (4)  made  analyses  of  low 
frequency  instability  (which  depends  upon  the  interaction  between  the  feed 
system  and  combustion  process)  based  on  a  constant  time  lag.  However,  it  was 
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impossible  to  explain  high  frequency  combustion  instability  using  a  constant 
time  lag. 

Crocco  (15)  introduced  the  +ime  varying  combustion  time  lag  in  his 
analysis  of  high  frequency  instability.  He  represented  the  total  time  lag 
as  the  sum  of  two  parts;  The  constant  "insensitive"  time  lag  ,  and  the 

"sensitive"  time  lag  V  ,  which  varies  wi+h  time,  responding  to  fluctuations 
in  the  chamber  conditions.  Thus 
Tv  =  tr.  +  X 

Rather  than  trying  to  describe  qiian+i tativel y  the  effect  of  eacn  of  the 
chamber  conditions  on  the  time  lag,  the  simplifying  assumption  was  made  that 
the  variations  of  all  physical  factors  could  be  correlated  with  values  of  the 
local  pressure.  For  small  amplitude  per furbations,  the  rate  at  which  the 
sensitive  time  lag  is  affected  can  be  expressed  by 

■C-'?  =  f 

where  the  barred  quantities  correspond  to  steady  state.  The  factor /V  was 
called  by  Crocco  the  "interaction  index". 

The  synthetic  representation  of  the  combustion  process  by  means  of 
the  time  lag  concept  eliminates  the  need  for  information  on  any  of  the  specific 
processes  occurring  in  the  chamber,  such  as  atomization,  vaporization,  .T.ix'ng, 
chemical  kinetics,  etc.  This  fact  constitutes  the  primary  advantage  of  the 
time  lag  type  analysis,  since  very  little  quantitative  information  is  available 
concerning  the  details  of  these  various  processes. 

The  application  to  the  sensitive  time  lag  theory  to  the  longitudinal 
mode  was  first  made  by  Crocco  in  1951  (15).  The  expanded  and  general izea  theory 
of  longitudinal  mode  combustion  instaii I ity  was  published,  but  without  experi¬ 
mental  verification,  in  1956  by  Crocco  and  Cheng  (I).  The  validity  of  thi  theo¬ 
ry  was  shown  in  a  series  of  experiments  reported  by  Crooco,  Grey  and  Harrje  (16). 
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In  th<'  iaiier  paper,  a  simple  means  for  measuring  the  combustion  parr-meters  ; 
and  through  stability  limits  tests  was  presented.  The  first  appiice^.ion  of 
the  sensitive  time  lag  theory  to  transverse  modes  was  made  by  Sea  I  a  (5),  who 
consiaered  the  s+ability  of  general  three  dimensionai  perturbations  on  a  one 
dimensional  steady  state  flow  situation,  ho  was  able  io  present  only  iimi+ed 
results  calculated  from  his  theory,  and  no  experimental  results,  although  he 
was  able  to  show  theoretically  the  destabilizing  effect  of  increasing  the 
nozzle  enfraPiC-s  I'.Kich  riumber.  ihe  above  rneriticned  investigations,  of  Croaco 
and  his  associates  (1,5,15,16)  provided  an  important  f''undation  for  the  present 
study  of  transverse  mode  combustion  instability. 

Maslen  and  Moore  (17)  investigated  theoretical ly  The  effects  of 
viscous  damping  in  a  chamber  without  combustion.  Their  results  showed  that 
large  amplitude  non-shock  type  waves  could  exist  in  the  spinning  form  of  the 
tangential  mode,  but  that  the  standing  mode  was  subject  to  damping  proportional 
to  the  amplitude.  The  damping  of  the  fundamental  longitudinal  mode  was  shown 
to  exceed  that  of  the  first  tangential  mode  for  low  values  of  the  '■atio  of 
chamber  length  to  diameter. 

Recently,  Pickford  and  Peoples  (18)  have  attempted  to  analyze  the 
stability  of  liquid  rocket  donuuSitloil  iby  means  of  .an  .avaiilable  enqrgy,, 
concept,  in  an  attemf*  at  an  a  priori  calculation  of  stsi .lity  behavior  from 
knowledge  of  the  physical  and  chemical  nature  of  the  propel lants, the  injection 
process,  etc.  However,  it  appears  that  the  accomplishment  of  this  objective 
still  lies  in  the  future. 

D.  Preliminary  exploratory  testing 

At  the  time  of  the  initiation  of  the  experimental  phase  of  this 


investigation  of  transverse  mode  instability,  the  theoretical  analysis  had 
not  been  carried  beyoru  that  of  Reference  2,  and  because  of  the  complexitity 
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of  the  computations  required,  s/ery  little  guidance  was  available  from  that 
source.  Therefor^  the,  experimental  program  began  with  a.  series  of  exploratory 
tests  designed  to  reveal  the  most  significant  parameters,  on  whici.  the  later, 
more  systematic  tests  were  to  concentrate.  The  first  series  of  exploratory 
tests  was  in  the  form  of  a  parameter  survey.  Later,  a  second  group  of 
experiments  using  half  chambers  ano  aiametral  barriers  was  made  in  order  ro 
evaluate  a  proposed  method  of  stability  limits  testing. 

The  experimental  apparatus  consisted  of  a  sectioned,  uncooled, 
copper  rovKot  motor  and  its  associated  hardwoi ond  ins-irumentation  (Figures 
3  and  4).  The  three  parts  of  the  motor — the  spud  type  injector,  the  four  inch 
long  cylindrical  chamber  section,  and  the  nozz I e-?-were  held  together  by  a 
hydraulic  cylinder  arrangement  to  facilitate  parameter  changes.  This  system 
had  been  used  successfully  on  the  longitudinal  program.  Tne  length  of  the 
cylindrical  section  was  fixed  by  the  so-called  "critical  length"  of  the  pro¬ 
pellant-injector  combination  (the  length  below  which  longitudinal  mode  inst- 
bility  is  not  possible),  since  it  was  desired  to  study  purely  transverse  tfodes. 
The  nozzle  consisted  of  a  30°  half-angle  convergent  part  blending  into  a  cir¬ 
cular  arc  throat  and  followed  by  a  15°  half-angle  divergent  cone.  Tha  throal 
diameter  and  the  subsonic  length  wore  determined  by  the  combination  of  "Ihrust 
level,  chamber  pressure,  and  chamber  diameter.  Al!  surfaces  exposed  :o  the 
hot  combustion  g^^ses  were  ceramic  coated  (0.012"  zirconium  oxide  over  0.003" 
nicket-chromium)  to  reduce  heat  transfer.  Sfarting  was  accomplished  by  means 
of  a  solid  propellant  squib  which  ignited  a  gaseous  hydrogen  preflow.  The 
hydrogen  flow  was  terminated  at  the  time  of  fuel  introduction.  The  pressuri¬ 
zed  feed  system  made  use  of  cavitaiing  venturis  for  flow  regulation.  Run 
durations  varied  between  one-half  and  three  seconds. 

Instrumenfation  included  the  usual  steady-state  equipment  (propellant 
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flow  rates,  chamber  and  injector  pressures,  etc.)  and  high  response  water 
cooled  strain-gace  pressure  transducers  (19)  flush  mounted  in  the  chamber. 

Six  locations  were  p'^ovided  for  the  high  frequency  transducers  in  the  cyliii- 
drical  chaniber  section,  in  two  axial  and  three  meridianal  planes,  in  order  to 
investigate  the  oscillation  wave  patterns.  The  steady  and  f •  uc^'-ating  compo¬ 
nents  of  each  pressure  signal  were  senarated.the  steady  state  portion  appear¬ 
ing  on  a  recording  potentiometer  and  the  fluctuating  portion  on  a  seven  channel 
FM  magnetic  tape  recorder.  Amplitudes  and  frequencies  of  the  recorded  pressure 
oscillations  were  then  measured  electronically  from  the  tape  rccofding.  A 
special  run  cutoff  device  provided  rapid  shutdown  in  thS  event  of  a  large  ampli¬ 
tude  instability,  in  order  to  minimizeburnout  of  the  hardware.  A  detailed 
discussion  of  the  instrumentation  is  found  in  (20i. 

The  exploratory  tests  were  designed  to  investigate  the  following 
parameters,  each  of  which  was  to  be  varied  independently  of  the  others; 

(1)  Chamber  dianietei — ^-In  the  transverse  mode  theory,  it  was  shown 
that  the  chamber  diameter  (or  radius)  ig  the  characteristic  length  associated 
with  transverse  modes  of  oscillation.  For  the  exploratory  tesis^  three  dia¬ 
meters  were  c.ioren,  in  a  range  practical  for  testing  at  the  forrestal  Research 
Center  rocket  test  facility.  These  diameters  were  3,  6,  and  9  inches. 

ever.  It  was  found  that  the  three  inch  diameter  equipment  could  not  designed 
with  the  flexibility  required  by  the  exploratd'",;  rc-ies.  Most  of  tie  experi¬ 
ments  were  conducted,  therefore,  with  the  two  larger  diameters. 

(2)  Chamber  pressure— The  effect  of  this  parameter  was  largely 
unknown  at  the  start.  Although  some  effect  on  the  combustion  parameters  (time 
lag  and  interaction  index)  was  expected,  no  clear  trends  had  been  shown  by 
longitudinal  mode  testing.  The  levels  selected  for  testing  were  150,  500,  aid 
900  psio. 

(3)  Thrust- — The  effecS  of  the  thrust  ievei  was  also  unknown. 
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Ir.jec+ion  veioci-tias  were  kept  cc  istanr,  and  the  thrust  was  varied  by  means 
of  the  injector  hole  size.  Some  effect  on  the  combustion  parameters  could 
be  expected  from  the  hole  size  changes  because  of  differences  in  fhe  atomi¬ 
zation  characteristics  of  the  injectors.  In  addition,  the  thrust  level, 
chamber  pressure,  and  chamber  diameter  determine  the  cltamber  Kach  number, 
which  had  been  shown  to  be  o:  -ortcf'ce  in  the  limited  res-j!'’.,  presented  in 

Reference  2.  Tests  were  made  at  thrust  levels  of  lOO,  500  and  SOOO  pounds. 

(4)  Mixture  ratio - The  longitudinal  stability  limits  testing 

demons tra+ed  that. the  combustior  parameters  were. strongly  apoendent.Qn„rtiJ;?tUP8 
ratio.  The  exact  physical  nature  of  this  dependence  remains  to  be  determined. 
Mixture  ratios  from  1.0  to  2.2  were  surveyed.  Included  in  this  range  were 
the  maximum  characteristic  velocity  (c*)  ratio  (1,4  for  the  propellants  used) 
and  the  stoichiometric  ratio  (2.09). 

(5)  Design  mixture  ratio— this  parameter  was  defined  as  that  ratio 
at  which  the  bipropellant  injection  velocities  were  equal  and  tfee  net  momen¬ 
tum  of  the  spray  produced  by  the  doublet  inj€5tor  was  axial.  Varying  the 
design  mixture  ratio  over  the  values  of  1.0,  1.4,  .1,8,  and  2,2  resulted  ir, 
changes  in  injector  hole  size  as  well  as  in  the  angles  of  the  Injecfor  pas¬ 
sages  with  respect  to  the  combustion  chamber  axis,  with  important  sfreda  on 
atomization  and  mixing  to  be  expected.  However,  the  influence  of  such  affects 
could  not  be  predicted.  Most  of  the  tests  wei;e  trade  using  the  1.4  design  mix¬ 
ture  ratio  injector Si 

(6)  Injector  t>^e - A  simple  unlike  impinging  doublet  type  injector, 

with  impingement  at  the  injector  face  .and  12  orifice  pairs,  was  chosen^ for  the 
initial  tests  of  the  exploratory  s-.-.rles.  Although. the  combustion  parameters 
had  been  clearly  shown  to  be  strong  functions  of  injector  type,  the  greatest 
amount  of  longitudinal  experience  had  been  obtained  with  such  an  Injector. 


(7)  Propellants — As  in  the  case  of  the  type  of  injector,  the 
propellant  combination  was  known  to  have  a  strong  influence  on  stability.  The 
combination  of  liquid  oxygen  and  95  percent  ethyl  alcohol  (Solox)  was  chosen 
for  the  first  series,  again  on  the  basis  of  greatest  instability  tasting  ex¬ 


perience. 

(8)  Injection  diameter — -Tn  the  exploratory  test  series,  vhe  12 
injection  points  were  equally  spaced  around  the  circumference  of  a  single 
circle.  Other  experimenters  had  noted  a  destabilizing  trend  as  the  diameter 
of  the  inject':;!'  circle  approached  tne  chamber  diameter.  Injectors  were 
designed  for  2  -  1/4,  5,  and  8  inch  circles.  However,  the  smallest  of  these 
proved  to  be  impractical. 

(9)  Orientation - For  maximum  experimental  flexibility,  spud  type 

injectors  were  used  (Figures  5  and  6),  Each  doublet  pair  was  contained  in  a 
single  spud,  which  could  be  oriented  in  either  of  two  ways,  90°  apart  (see 
Figure  7).  In  the  "radial"  orientation,  each  spud  was  arranged  so  that  the 
line  of  centers  of  the  injection  orifices  lay  along  a  radius.  Kith  the 

spud  rotated  through  90°  to  the  "tangential"  orientation,  the  line  of  centers 
was  tangential  to  the  injection  circle.  The  spud  or ienta  i  I  ..n  proved  To  be 
one  of  the  most  irhportant  parameters. 

The  exploratory  testing  was  begun  at  the  100  pound  thrust  revel. 
However,  the  results  were  somewhat  indeterminate  When  oscillations  were 
present,  the  amplitudes  were  very  small.  In  addition,  reliable  starting  was 
a  problem  because  the  small  throat  diamerers  necessitated  the  use  of  a  spark 
■  gn  I  I  er  r  a  1 1  lei  i  hati  1 1 so  I  i  o  ^r  opo  1 1  an  i  squ  i  b  used  in  I  ong  i  tud  i  no  I  mode 
stabiliTy  tests  (16).  Also,  with  the  small  propel  I  ant  flows  involved  in 
these  tests,  considerable  trouble  was  had  vrith  clogging  of  the  caviliating 
venturis.  Therefore,  after  a  limited  number  of  tests,  the  100  pound  level 

i 

was  discarded  in  favor  of  the  .mere  reasonable  500  pound  thrust  level. 


Six  series  of  tests  comprised  the  exploratory  testing  at  the  500 
pound  thrust  level.  These  series  were  as  follows: 

1  -  9"  chamber,  5”  injector,  spuds  tangential 

2  -  9“  chamber,  S*'  injector,  souds  tangential 

3  -  6"  chamber,  5”  irij&rtor,  spuds  tangentiaj 

4  -  9"  chamber,,  5”  injector,  spuds  radial 

5  -  9'i  chamber,  8"  injector,  spuds  radial 

6  -  6"  chamber,  5"  injector,  spuds  radial 

In  each  scries,  a  mixture  ratio  survey  was  made  (the  desicn  mixture  ratio  was 
1.4  for-  all  series)  at  each  of  three  chamber  pressure  levels,  150,  500,  and 
900  psia.  The  results  of  these  tastsare  shown  in  Figures  8  and  9. 

In  Figure  8,  the  test  points  of  series  I  to  3  are  plotted  on  the 
chamber  pressure,  mixture  ratio  plane,  with  the  nature  of  the  oscillation 
shown  by  the  shading  of  the  points  -It  can  be  seen  that  all  tests  with  the 
tangential  orientation  were  unstable  in  the  first  tangential  mode.  Many  of 
the  tests  caused  destruction  of  the  pressure  transducer  or  of  chamber  or 
injector  parts,  in  spite  of  the  run  shutdown  device.  The  peak-to-peak  amp¬ 
litudes  of  the  pressure  flutuations  varied  between  I /J  and  330  psi  on  these 
tests.  Figure  10  shows  the  wave  shapes,  and  Indicates  the  use  of  multiple 
pickups  to  determine  the  wave  pattern  within  the  combustion  chamber.  It  was 
observed  that  the  oscillations  in  these  tertr  were  all  of  the  spinning  type. 
This  fact  can  be  seen  in  Figure  10  in  that  the  90°  phase  lag  of  trace  2  with 
respect  to  trace  I  corresporas  exactly  to  the  spacing  of  the  two  pressure 
transducers  around  the  circumference  of  the  chamber.  In  addition,  it  can  be 
seen  that  the  pressure  oscillation  near  the  nozzle  is  smaller  and  less  steep¬ 
ened  than,  and  lags  slightly  behind,  the  oscillation  near  the  injdctor.  This 
result  is  to  be  expected,  since  it  has  been  shown  that  the  combustion  zone, 
with  this  particular  injector  and  propellant  combination,  is  concentrated 
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within  2"  of  the  injector  face.  None  of  the  waves,  even  close  to  the  injector, 
showed  steepening  into  shock-type  waves,  in  contrast  with  those  observed  on 
longitudinal  tests  (Figure  II).  Although  all  tests  with  the  tangential  orientation 
were  unstable,  the  amplitudes  of  the  fully  developed  oscillations  were  greater 
when  the  propellants  were  injected  near  +he  chamber  wall  (Figure  8),  thus  showing 
one  effect  of  the  injection  diameter. 

Quite  different  results  were  obtained  with  the  radiai  spud  orientation 
as  shown  in  Figure  9.  The  first  observation  which  can  be  made  is  that  the  radia! 
orientation  te^is  were  more  stable  than  the  tangential.  Only  th?  9"  chamber  with 
the  8"  injector  gave  well  defined  zones  of  instability.  The  first  tangential 
mode  was  observed  at  the  lowest  pressure  level  for  mixture  ratios  in  the  vicinity 
of  the  design  point  (1.4),  The  disappearance  of  the  first  tangential  mode  can 
be  explained  qualitatively  by  the  fact  (determined  by  the  analysis  of  Reference 
5)  that  decrease  in  the  chamber  Mach  number  is  stabilizing.  Since  all  tests  of 
this  series  were  made  at  the  500  pound  thrust  level,  the  higher  pressures  were 
obtained  by  decreasing  the  size  of  the  nozzle  throat,  thus  decreasing  the 
chamber  Mach  number,  /it  higher  pressures,  the  second  tangential  mode  was  ob¬ 
served,  but  with  very  much  smaller  amplitudes.  No  calculations  were  made  in 
Reference  2  for  the  second  tangential  mode.  However,  the  presence  of  the  ref-ond 
mode  in  these  tests  agrees  well  with  theoretical  results  presented  in  a  la+er 
section.  It  can  be  seen  that  only  for  this  series  (5'  were  well  defined  sta¬ 
bility  limits  observed,  to  which  a  linear  theory  could  be  applied. 

Injection  of  the  propel  lent  nearer  to  the  center  of  the  combustion 
cr-?mber  was  seen  in  series  4  tests  (9"  chamber,  5”  injector,  radiai  orientation) 
to  have  a  stabilizing  effect,  a  result  which  agrees  qualitatively  with  the  re¬ 
duced  amplitudes  of  the  corresponding  tangential  orientation  tests  (series  I). 

The  results  of  the  6"  chamber,  5"  injector  tests  (radial  orientation)  were 
not  ciear.  At  the  low  pressure  ievei  stability  was  observed.  At  higher 
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pressures  -The  traces  showed  several  modes  present  intermittently,  including 
the  rirst  and  second  tangential  and  first  radial  modes.  The  stabili+y  of  the 
6"  motor  as  contrasted  to  the  instauility  observed  for  similar  operating  con¬ 
ditions  wifh  the  9"  motor  can  be  explained  qualitatively  by  noting  that  while 
the  resonance  properties  of  the  9"  chamber  were  matched  to  the  time  lag  of  the 
combustion  process,  the  decrease  in  diameter  destroyed  the  resonant  coupling. 

At  this  point,  the  need  for  more  systematic  experiments, 
involving  the  determination  of  stability  limits,  which  could  be  compared 
quantitatively  with  the  theory  being  developed,  became  clear.  In  the  longi¬ 
tudinal  program,  the  primary  experiment  is  tnat  of  determining  limits  of 
staeility  on  the  mixture  ratio-chamber  length  plane.  Because  of  the  fact 
that  the  length  of  the  chamber  (above  a  certain  minimum)  has  little  effect 
on  the  combustion  pattern  but  a  grea-i-  effect  on  the  resonance  properties  of 
the  chamber,  such  limits  can  be  used  to  deduce  values  of  the  combustion  para- 
meiers,  the  sensitive  time  lag  "C  and  the  interaction  index  /7U  (16). 

For  the  case  of  the  tangential  modes  of  transverse  instability, 
the  "sector  motor"  concept  was  evolved.  In  this  type  of  test,  only  one 
sector  of  the  motor  is  operated,  the  rest  of  the  cha'-'he'*  and  no2;2!s  (to  the 
throat)  are  filled  in  with  a  copper  insert.  Such  a  motor  is  shown  schemati¬ 
cally  in  Figure  12(b)  and  in  Figure  13. 

A  180°  or  half  chamber,  formed  by  means  of  an  insert,  can  oscillate 
only  in  the  standing  form  of  a  iangential  mode.  On  the  other  hand,  a  full 
360°  chamber  can  sustain  either  the  standing  or  spinning  form  of  the  mode. 

As  in  the  tests  discussed  above,  the  spinning  form  is  the  one  normally  ob¬ 
served  in  an  unstable  full  chamber.  According  to  the  pressure  sensitive  time 
lag  theory  of  (5),  the  stability  behavior  should  be  the  same  for  both  thr 
full  and  half  chambers,  as  long  as  the  operating  conditions  are  the  same.  A 
series  of  tests  was  made  in  order  to  test  this  cone  I  u« 'op,  ■''hs  resul+s  are 
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presented  in  Table  I.  All  tests  were  made  at  the  same  conditions  of  chamber 
pressure  and  mixture  ratio,  but  it  is  clear  that  the  stability  behavior  of  the 
two  forms  of  first  tangential  mode  is  not  the  same. 

Further  tests  in  this  series  made  use  of  barriers,  of  1/4”  thick 
steel  and  of  various  lengths,  placed  'n  the  chamber  along  a  diametral  surface, 
instead  of  the  180°  sector  insert  (see  Figure  12(a)).  Such  barriers  were 
found  to  be  capable  of  producing  the  same  stable  operation  as  the  full  180° 
insert.  It  was  also  determined  that  stable  operation  resulted  with  a  barrier 
only  5/8”  long,  if  the  barrier  was  placed  against  the  injector  face.  If  a 
space  was  left  between  the  injector  and  barrier,  instability  was  observed> 
even  though  a  longer  barrier  was  used.  The  conclusions  to  be  drawn  are  that 
the  pressure  sensitivity  theory  of  (5)  is  not  sufficient  to  explain  completely 
transverse  modes  of  instability,  and  that  the  additional  effects,  which  cause 
the  spinning  form  to  be  much  more  unstable  than  the  standing  form,  are  con¬ 
centrated  near  the  injector  face,  in  the  early  combustion  zone. 

In  order  to  investigate  the  transition  from  the  stable  standing 
form  of  the  tangential  modes  to  the  unstable  spinning  form,  the  permanent 
steel  barrier  was  replaced  with  a  destructible  one.  The  ia^.  far  noS  of' 
aluminum,  1/8"  thick  by  I"  long,  and  was  placed  against  the  injec  fop  face  # 
Experiments  showed  that  the  aluminian  barrier  remained  essentially  intact  for 
the  first  second  of  operation,  but  disappeared  quickly  at  that  time.  Thus 
each  aluminum  barrier  test  could  show  motor  operation  under  both  stable  and 
unstable  conditions. 


A  comparison  of  the  stability  of  the  standing  and  spinning  modes 


using  the  tangential  spud  orientation  is  shown  in  Figure  14.  The  first  and 
second  lines  represent  the  whole  and  half  chamber  tests,  respectively,  which 
were  discussed  above.  The  third  and  fourth  lines  present  the  results  of  the 
aluminum  barrier  tests.  It  can  be  seen  that  the  early  portion  of  these  tests 
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corresponcfed  to  the  180°  sector  tests,  and  that  as  soon  as  the  constraint  to 
The  standing  irode  was  removed,  the  rocket  operation  became  unstable.  In  the 
transition  period  a  smooth  increase  in  amplitude  from  the  combustion  noise 
level  io  that  of  the  fully  developed  oscillation  was  observed. 

Results  of  a  similar  series  of  experiments  using  the  radial  spud 
orientation  are  presented  in  Figure  15.  Just  as  in  the  tangential  orientation 
tests,  the  unstable  region  observed  with  full  chamber  operation  was  not  pre¬ 
sent  in  the  half  chamber  tests.  The  aluminum  barrier  tests  reproduced  the 
half  chamoer  test  results  during  the  early  portion  of  each  run  it  <  t  second). 
However,  unlike  the  corresponding  tangential  orientation  tests,  when  the 
barrier  burned  away,  the  rocket  motor  operation  remained  stable.  Such 
anomalous  behavior  may  indicate  the  existence  of  non-linear  effects  during 
the  starting  transient  period  of  the  radial  orientation  tests. 

The  exploratory  test  program  demonstrated  three  effects  not  pre¬ 
dicted  by  previous  theories.  They  were:  (I)  the  marked  difference  in  sta¬ 
bility  behavior  between  the  tangential  and  radial  Foud  orientations,  (2)  the 
greater  tendency  to  instability  of  the  spinning  mode  compared  to  the  standing 
mode,  and  (3)  the  effect  of  the  injection  diameier  ii-  ar,  injei-i'ion  pattern 
concentrated  on  a  single  circle.  It  then  became  necessary  to  r.iodif.'  tne 
existing  theory  in  an  attempt  to  explain  the  effects  noted. 

These  modifications  are  developed  in  the  following  two  chapters. 

In  Chapter  II,  the  stability  analysis  is  made  for  a  one  dimensional  steady 
state  flow,  which  implies  that  the  injection  is  uniformly  distributed.  The 
geometry  of  I  tie  combustion  chamber  is  lhai  of  a  circular  cylinder,  cor:  e- 
sponding  to  the  full  chamber  used  in  the  exploratory  tests.  The  analysis 
is  extended  in  Chapter  III  to  cover  cases  of  non-uni formly  distributed  .n- 
jection  and  is  also  applied  to  the  sector  motor.  Chapter  IV  presents  ex¬ 
perimental  data  which  is  compared  with  the  theoretical  predlctionsj  the 
validity  of  the  modified  theory  is  then  assessed. 
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C^  PTER  II 

Instability  Analysis  for  a  Circular  Cylindrical 
Combustion  Chamber  with  Uniform  Injection 

A.  Governing  Equations 

The  analysis  which  follows  is  based  on  tne  consideration  of  the 
stabilify  of  small  pert Jrbations.  Each  of  the  quantifies  describing  the  ‘low 
in  the  rocket  motor  is  assumed  to  oscillate  about  a  given  steady  state  value. 
For  certain  operating  conditions,  the  oscillation  amplitudes  will  increase 
with  time,  even  for  arbitrarily  small  initial  amplitudes.  Since  random  fluc¬ 
tuations  of  snail  but  finite  amplitude  (generally  referred  to  as  "combustion 
noise")  are  always  present  in  a  rocket  motor,  such  unstable  operating  condi-  . 
tions  must  be  avoided.  The  approach  which  will  be  used  in  the  determination 
of  the  stability  of  a  rocket  combustion  system  is  the  following.  Conditions 
will  be  established  far  the  existence  of  neutral  oscillations,  the  amplitudes 
of  which  neither  increase  nor  decrease  with  time.  The  assemblage  of  the 
operating  points  of  the  system  which  satisfy  the  conditions  for  neutral  oscil¬ 
lation  form  the  "stability  limits"  which  divide  the  '.'ns+3>'i<;  from  the  ftable 
regions  of  operation.  The  complete  stability  behavior  of  the  rocket  is  des¬ 
cribed  by  specifying  the  stability  limits  and  indicating  the  unstable  regions. 

The  rocket  motor  to  be  considered  consists  of  a  circular  cylindrical 
combusfon  chamber  of  length  and  radius  r^,  followed  by  a  converging- 
diverging  expansion  nozzle.  The  walls  are  assumed  to  be  rigid,  impervious, 
and  adiabatic.  The  medium  within  iiie  chamber  consists  of  gaseous  products 
of  combustion  in  which  are  distributed  burning  droplets  of  the  liquid  pro¬ 
pellants.  The  volume  occupied  by  the  liquid  is  take.,  to  be  negligibly 
small.  In  this  two  phase  flow,  the  burning  droplets  act  as  distributed 
sources  of  mass,  momentum,  and  energy.  It  is  assumed  tiat  con’.bus‘'ion  is 


completed  within  the  combustion  chamber  and  that  the  gas  velocity  in  the 
chamber  is  so  small  that  the  square  of  the  Mach  number  is  negligible  compared 
to  unitv'i  Thus,  in  this  Ideal  rocket  motor,  the  flow  is  divided  into  two 
parts,  a  low  speed  combustio)^  process  and  an  expansion  to  supersonic  velocity 
without  combustion. 


It  (-3  convenient  to  carrv  out  the  analysis  in  terms  of  the  follow¬ 


ing  dimensionless  variables: 

pressure  = 


d^sl  ty 
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s'jbscript  o  denotes  >/alues  at  the  injector  face,  *  indicates  a  dimensional 
quantity,  a  superposed  arrow  indicates  a  vector,  and  barred  quantities  corre¬ 
spond  to  steady  state.  Droplet  properties  will  be  denoted  by  the  subscript 
L.  It  should  be  noted  that  the  droplet  den'.ity  is  the  total  mass  of 
droplets  per  unit  volume  of  gas.  Since  the  steady  state  gas  velocity  must 
vanish  at  the  injector  face,  j  and  5^'  are  stagnation  quan¬ 

tities.  The  cyiindrical  coordinate  system  is  shown  in  Figure  l6. 

The  equations  stating  the  conservation  of  mass,  momentum,  and  erergy 
for  the  system  of  burning  droplets  and  combustion  products  will  be  presented 
here  without  derivation.  For  a  more  complete  discussion,  the  reader  is 
referred  to  Reference s( I )  and  (5).  The  continuity  equations  are 


In  the  energy  conservation  equation,  terms  expressing  the  work  done  by 
viscous  stresses  and  heat  transferred  by  conduction  or  diffusion  will  be 


neglected.  Then  the  energy  equation  in  terms  of  the  stagnation  enthalpy. 


0 


0.V-) 

ro-.  zr  V 


is 


L  St 


{2,20 


-20- 


where-^2  includes  the  chemical  e'ergy  of  the  propellants. 

The  momentum  intercnange  between  the  liquid  and  gaseous  phases  can 
be  expressed  by  the  equc+ion 

(2.2d) 

where  the  coefficient  -ft.  is  related  to  the  frictional  drag.  The  equation 
describing  the  energy  transfer  between  phases  will,  following  Crocco  (I), 
be  replaced  by  the  simple  relation 


0 

The  equation  of  state  of  the  combustion  gases  is  taken  as 

^  c  />T 


(2.2e) 


(2.2f) 


The  set  of  equations  to  be  solved  must  be  completed  by  an  expression 
for  the  burning  rate  Q,  to  be  derived  in  a  later  section,  and  appropriate 
boundary  conditions  on  all  of  the  variables. 

Each  quantity  is  now  written  as  the  sisn  of  ?  s+eady  state  part  and 
a  time  varying  perturbation.  Thus  the  pressure  is  expressed  as 


(2.3) 


where  is  the  sicady  staie  pressure  and  “Jk*  is  the  complex,  space  dependent 
perturbation  amplitude.  (For  simplicity  in  the  analysis  below  the  term 
"pressure  perturbation"  will  be  understood  to  mean  -ff  only). 
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In  the  time  dependent  factor  in  which  the  real  part  X 

is  the  amplification  factor  and  the  imaginary  part  is  the  frequency 
of  oscillation.*  Tne  condition  X=0  corresponds  to  neutral  oscillation  at 
a  stability  limit,  while  X<o,A>d  represent  stability  and  instability, 
respectively.  As  mentioned  above,  the  analysis  will  be  directed  toward 
the  determination  of  the  conditions  which  prevail  when  or  It  is 

assumed  that  the  perturbations  are  sufficiently  small  that  the  squares  and 
products  of  perturbations  can  be  neglected  in  comparison  to  the  steady  slate 
va  I  ues . 

The  components  of  the  vectorial  gas  and  droplet  velocities  V 
and  Vl  ,  are  represented  as  follows: 

axial  components  U-^  u.j_ 

radial  components  TTJ  n;*. 

tangential  components  itf )  ‘*^1 

The  steady  state  flow  is  assumed  to  be  one  dimensional,  although  the 
perturbations  will  be  considered  in  three  dimensions,  thus 

V  -  n  ,  =  u,  7 

{ 

and  the  time  and  transverse  space  derivatives  vanish  in  th^  steady  state 
governing  equations. 


*In  the  complex  representation  of  oscillating  quantities,  used  in  this 
analysis,  the  actual  physical  value  of  a  cuantity  is  given  by  the  real 
part  of  the  corresponding  complex  expression. 
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The  s/s+em  of  steady  stat'-  equations.- i s  then 


(2.5a) 

(2.5b) 

(2.5C) 

^  ^  f  f 

(2.5d) 

u. 

u-e 

(2.5e) 

di 

(2.5f) 

with  the  iol lowing  boundary  conditions  at  i  =  0  ; 

u  =  o  ,  Dl  =  \o  »  Ao 

>  =  =  T  =  I 

"-^so  »  "^Ls  ~  "^Lso  '  9  ~  0 

The  continuity  equations  are  integrated  to  give 

-ft 

0 

where  fto^  is  the  propellant  injection  rate.  The  momentum  equation 
yields 


^  ^  1-  ^((O  ) 


Then, 

by 


since 


f= 


- 1 


f  the  steady  state.'.gas,  density-jjs:  given 


(2,8) 
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Ths  dimenEionless  gas  vsloci+y  u,  which  i?'  very  nearly  equal  to  the  Mach 

number,  increases  from  zero  at  the  injector  face  to  a  maximum  value  Ug  at 

I  ^ 

the  chamber  exit  (or  nozzle  entrance)  station,  z  =  z„  =  .  It  is 

e 

assumed  that  the  exit  velocity  Ug  is  small  enough  thar  Its  square  is  negligi¬ 
ble  compared  to  unity.  That  is,  the  analysis  will  be  carried  ouf  to  include 
terms  of  the  same  order  of  magnitude  as  ex!-*"  velocity  (written  symbci’- 
cally  as  0(Ue)  ),  but  terms  of  fhe  order  of  Ue^  will  be  neglected. 

The  order  of  magnitude  comparisons  to  be  iTade  In  the  following 
analysis  are  not  .n+ended  fo  be  mathematically  exact,  but  rather  to  express 
relations  which  are  approximately  true  in  the  range  of  parameter  vatues 
considered. 

Liquid  injection  velocities  are  typically  somewhat  less  than  the 
chamber  exit  velocities.  Initially,  the  droplet  velocity  decreases  because 
the  gas  velocity  is  very  nearly  zero  in  the  vicinity  of  the  injector.  For 
injection  into  stagnant  gas,  equation  (2.5e)  gives 


The  distance  which  the  droplet  penetrates  into  gas  is  then 


A  reasonable  assumption  is  that  the  "penetration  distance"  2^  is  0(1), 
so  that  the  droplet  coefficient  -H,  is  0(Ug),  Therefore,  the  droplet 
velocity  decreases  gradually,  as  the  gas  velocity  increases  due  to  combustion, 
until  the  two  velocities  are  equal.  Thon  15^  increases  again  as  the  droplets 
are  carried  along  with  the  gas  stream  toward  the  nozzle.  The  droplet  velocity 
is  always  less  than  the  gas  velocity  and  approaches  it  asymptotically.  It 
is  clear  that  the  momentum  interchange  coef^ic'ent  will  Increase  as  the 
droplet  vaporizes  and  burns.  Rather  than  take  this  variation  into  account 
specifically  a  constant,  effective  value  will  be  used. 


The  steady  siate  values  of  pressure,  temperature,  and  density  are. 


to  0(Ug)  ,  _ 


(2.9) 


The  relation  between  the  gas  valocity  and  the  burning  rate  is 


3 

Q.(i'  )oc  ’ 


or 


—  ^ 

^  "  df 


Integration  of  equation  (2,5f)  yields 
=  constant  = 


'LSi, 


Introducing  this  result  into  the  energy  equation  (2.5c)  gives 


-  „  dS/jS' 

whicn  can  be  written  as 


^  ]  -  O 

and  integrates  to 

tl  ^LSo  )  ”  C:0hS+O«  + 


(2.!0a) 


(2.10b) 


The  only  constant  >o:i;8  which  will  satisfy  this  relation  is  zero.  Therefore 

*^5  -  ~  ^LS  -  (2.11) 

In  order  to  investigate  the  stability  problem  +or  a  pa-ticatar 
rocket  motor,  if  is  necessary  to  have  specified  the  steady  state  combustion 
distribution  function  Q  (Z)  ,  equivalently,  the  gas  velocity  distribution 
u(£)  .  The  liquid  droplet  velocity  can  be  found  in  terms  of  the  gas  velocity 
from  equation  (2.5e),  as  discussed  in  Appendix  A.  Then,  since  it  is  assumed 
that  combustion  is  complete  at  +he  chamber  exit,  the  liquid  density  Is 
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obtained  from  the  continuity  equation  by  the  relation 

tTe  -  -.I 


Pl  " 


Uj-ir} 


(2.12) 


In  order  to  obtain  the  solution  for  the  unsteady  flow,  the  steady 
state  solution,  equations  (2.9)  to  (2.12),  is  perturbed  by  writing 
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(2.13) 


as  discussed  above.  The  exr  s  (2.13)  are  then  inserted  into  the 

equations  governing  the  flow,  equations  (2.2a)  to  {2.2f),  the  indicated 
differentiations  and  multiplications  are  carried  out,  and  the  terms  involving 
squares  and  products  of  perturbations  are  neglected.  Then  subtraction  of  the 
steady  state  equations  and  cancellation  of  the  factor  e^^  give  the  equations 
to  be  satisfied  by  the  perturbations  (i.  e.,  the  space  dependent  pei  turbation 
amplitudes).  These  steps  are  presented  in-  detail  in  (5)  and  will  not  be  re¬ 
peated  here.  The  resulting  perturbation  equations  (in  cylindrical  coordinate'-) 
are  the  fol lowing: 

conti nui ty. 


(2.14a) 
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axial  component  of  moment  m, 
c(a'+ap')+  u 


rad  i  a  I  component  of  moeme'ntum 


■»  0,  •«-'  =  .-;  §>' 


3£ 


tangential  component  of  momentum, 

S  (.  w'  «r,'  ^  If  +  3^'  t^' 


3-j 


j  _  I  a-^' 

^  ^  a? 


energy, 

.  -  3<J 


droplet  dynamics. 


.r  au/r' 


state. 


i?'  =/’'■+  T  ' 


and  droplet  energy. 


■?  ^Ls  ^  •- 


as 


rr  <5 


This  last  equation  is  satisfied  by 
i  /  _  /)/  /  , 


« 

y 


{2.14b) 


(2. 14c) 


(2.l4d) 


{2.l4e) 


(2.l4f) 


(2.l4g) 


(2.l4h) 
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Assuming  that  the  enthalpy  (or  energy)  of  the  liquid  propellants  supplied 
to  the  rocket  motor  is  constant, =0,  and  thus  =  0  .  Introducing 

this  result  into  the  energy  equation  and  replacing  the  stagnation  enthalpy 
perturbation  by 

(since  ■^'  =  for  the  nondimensional izarion  scheme  used,  and  The  equarior; 
I  ^  f  ' 

of  state  gives  T  -  y  ~  f  yields  a  modified  energy  equation. 

sfy  -(o'+  /is'-Oixu'l  4-^  ^1b‘ 

The  equations  describing  the  dynamics  of  the  liquid  droplets  are 
discussed  in  Appendix  A.  Solutions  are  obtained  of  the  form 


where 

oscillations  are  an  order  of  magnitude  smaller  than  the  corresponding  gas 
velocity  oscillations. 

B.  Burning  Rate  Perturbation 

The  system  of  perturbation  equations  is  not  complete  without  an 
expression  for  the  burni  q  rate  perturbation,  Q' .  The  relc.  lon  for  Q'  derived 
by  Crocco  (I)  and  later  used  by  Sea  la  (5)  made  use  of  the  concept  of  the 
combustion  time  lag.  That  is,  the  gradual  evolution  of  combustion  products 
from  an  element  of  propellant  mixture  traveling  through  the  conibustion  chamber 
after  injection  was  approximated  by  a  discontinuous  conversion  from  reactants 
to  products  after  a  certain  time,  called  the  time  lag,  had  elapsed.  By  this 
approximation,  the  impossible  task  cf  describing  the  combustio'’  process 


U. 
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'  S  'IT  ’ 


(2.I5) 


(i-^  ®  ^  ).  Since  and -j?/  are  0(Ug),  the  liquid  velocity 


(2.i4l 
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throuqh  the  quanti rative  knowledge  f  r  its  intermediate  history  was  greatly 
simplified.  It  was  recognized  that  if  the  combustion  chamber  co’^ditions 
fluctuate,  the  time  lage  must  also  vary.  An  additional  approximation  was 
introduced  to  simplify  the  description  of  the  effect  of  the  chamber  condi¬ 
tions  on  the  time  lag.  The  total  time  lag  'Cj  was  taken  as  the  sum  of  an 
insensitive  time  lag  and  a  sensitive  :ime  lag  f  .  The  various  ,,hyslca! 
conditions  were  assumed  to  affect  the  combustion  processes  only  during  the 
sensitive  portion  of  the  total  time  lag.  Thus,  during  fluctuating  operation 
of  the  rock's:  motor,  the  time  variation  of  the  total  combustion  tims  lag  is 
taken  up  completely  by  the  sensitive  tims  lag;  the  insensitive  time  lag  is 
constant. 

The  burning  rate  perturbation  can  be  related  to  the  rate  of  change 
of  the  time  lag  by  considering  a  small  fractional  element «A!i^  of  the  propel¬ 
lants  burning  at  a  given  location  at  the  time  t.  The  burning  rate  of  such 
an  element  is  related  to  the  local  burning  rate  per  unit  volume  by 

6^  t)  =  9,i)' JO  di 

This  element  was  injected  )  seconds  before  combustion,  as  a  frac¬ 

tion  of  the  total  injection  rate 

Then,  since  the  element  burning  during  the  time  intp-va!  dt  was  injected 
during  the  inter''=*l  ,  continuity  -equires  ^raT 


or 
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Introducing  the  division  of  the  total  time  lag  into  its  sensitive  and 
insensitive  parts. 


Ty  =  T;  4  t 


CTi- 


where,  by  definition,  ^'’  =  0  »  gives 

Under  steady  state  conditions,  this  last  relation  reduces  to 


In  the  following  analysis,  attention  will  be  concen+rated  on 
intrinsic  instahi!:+y,  That  is,  instability  occurring  when  the  injection  rate 
is  constant.  Then 


and 


(T 

<r 4*7^  :=  S 


In  terms  of  the  burning  rate  per  unit  volume. 


Mi 


—  ^  (  I  “■  ^  Order  ‘fcrynS 


or 


_  Q.Q  -- 
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Spatial  effects  due  to  the  fluctuations  of  the  total  time  lag  and  liquid 
droplet  velocity  contribute  terms  of  the  same  order  of  magnitude  as  u  . 
WHile  ^  may  be  0(1)  locally,  the  effect  of  this  factor  in+agreted  over 
the  entire  chamber  is  on!>  0(Uq  ).  This  result  is  demonstra+,d  Appendix  B. 

_  7 

The  terms  introduced  by  spatial  variations  are  therefore  0(Ug“)  and  have  been 
neglected. 


In  the  original  formulation  by  Crocco,  the  time  lag  wa‘^  defined 
by  noting  that  the  transformation  for  burnt  gases  takes  place  only  when  the 
processes  occurring  during  the  time  lag  have  accumulated  up  to  a  .'ell 
determined  level,  This  conditio'  was  expressed  by  the  eq-uation 


.16) 
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f 

t-r, 

where  is  the  overall  rate  of  the  various  processes  involved  in  the  com¬ 
bustion  of  the  liquid  bi propel  I ants,  and  £  is  the  time  a+  which  +he  trans¬ 
formation  to  burned  gases  occurs.  In  the  integrals,  the  rate-|^  must  be 
computed  using  the  instantaneous  values  of  the  physical  conditions  at  each 
instant  -fc'  at  the  location  X' (representing  the  spatial  variables  ) 

where  ftp  element  is  at  that  instant.  The  sensitivity  of  the  rate  ^  to 
variations  in  chamber  conditions  (correlated  to  the  local  pressure)  was 
expressed  by  the  interaction  index /rw,  defined  by 


i-Cr 


(2.17) 


3J-  91 

a-f*  dT  ^t> 


(2.18) 


where  the  parenthetical  expression  is  evaluated  at  the  steady  state  operating 
conditions.  Tiie  assumption  was  made  that  the  sensitive  portion  of  the  time 
lag  occurrs  Just  before  the  conversion  to  combustion  products.  Then  during 
the  insensitive  time  lag  the  unsteady  process  rate  is  equal  to  the  steady 
state  rate. 

Since  the  unsteady  pressure  is  of  the  form  "jr  CT  ,  the  process  rate 

during  the  sensitive  time  lag  is  given  by 

)  >  t'Z 


(2.20) 
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Introducing  thesv.-  expressions  for  ■ ;  into  equation  (2.17)  gives 


^-'r  <-rr  i-t 

pt  i 


The  first  term  in  this  equation  can  be  wri 

.t 


j 


which,  by  equation  (2.17)  is 

ff‘“'  *  L 

i-Xr 

Then  equation  (2.21  becomes 

t-Xy 


v-tjr 


(2.21) 


(2.22) 


Assuming  the  spatial  nonuni formity  of  to  be  negligible,  the  integral  on 

the  left  side  of  equation  (2.22)  can  be  written  as  t 

l5r  *' 

=  T-I  =  - 


neglecting  higher  order  terms.  Differentiation  of  equation  (2.23)  wi iii 
respect  to  time  gives 


dt  ’-'P  (k  d-tJ)  e 

which,  can  be  rewritten  as 

<  >Yx'4Ct-i:3  -5T 


'  L  yc’^av  ] 

Finally,  combining  equations  (2.16  and  (2.24),  the  burning  rate,  perturbation 


(2.24) 


is  given  by 


f  i-  -ST  ?  ,  /  .  -  ... 


(2.25a) 
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which  i?  of  the  form 


a'= 


(2.25b) 


In  the  preliminary  exploraiory  tests,  the  o'"'entation  of  the  spray 
fans  (produced  by  unlike  impinging  donhle+s  arranged  around  fhe  circumference 
of  a  circle)  was  found  to  be  critical  with  regard  to  transverse  mode  instabil¬ 
ity.  Such  a  result  could  not  be  predicted  using  the  burning  rate  expression 
of  equation  (2.25a),  in  the  derivation  of  which  it  was  assumed  that  all 
quantities  affecting  fhe  combustion  process  could  be  correlated  to  the  local 
pressure.  While  this  assumption  can  be  applied  to  thermodynamic  state  propei — 
ties  such  as  temperature  and  density,  it  can  not  be  applied  to  the  components 
of  the  gas  velocity,  which  may  be  of  equal  importance  in  the  combustion  pro¬ 
cess.  The  effects  of  tne  velocity  perturbations  must  therefore  be  taken  into 
account  separately  from  the  effects  of  those  quantities  which  can  be  correlated 
with  the  pressure. 

Only  the  transverse  (i.  e.,  radial  and  tangential)  cornponents  of 
the  velocity  perturbation  can  be  significant  in  the  explanation  ot  tne 
orientation  effect.  For  purely  transverse  modes  of  oscillation,  the  lonq’- 
tudinal  component  is  always  much  smaller  than  the  transverse  ccmponeni » , 

In  addition,  the  longitudinal  velocity  perturbafon  vanishes  at  the  injector 
face  and  has  its  smallest  magnitude  in  the  early  combustion  region,  the  zone 
in  which  the  diametral  barriers  wore  seen  in  the  preliminary  testing  to  have 
the  greatest  influence  on  stability. 

Of  the  various  intermediate  prf' :esses  occurring  during  the  combustion 
of  liquid  bipropellants,  those  most  sensitive  to  velocity  are  the  vaporizaticn 
of  the  liquid  droplets  agd  the  rnl>(ing  of  thd  vapgrUsd  propellants  which  must 
precede  chemiqaj  reaction.  In  the  proe§§8  of  vapopUation,  the  important 
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quantity  is  the  absolute  nvignitude  of  the  ve  ocity  of  the  hot  gases  relative 
to  the  liquid  droplet.  The  differences  'n  the  veiocity  paTterns  of  the  cfand- 
ing  and  spinning  forms  of  tangential  modes  of  oscillation  (which  were  obser-ed 
to  have  strikingly  different  stability  characteristics)  could  result  in  differ¬ 
ences  in  the  local  and  overall  vaporization  rates.  However,  since  the  direction 
of  the  relative  velocity  is  immaterial.  The  oi  ic..;ation  effect  is  not  e.v,;,.ii- 
able  by  means  of  the  vaporization  process.  In  addition,  the  theoretical 
study  of  Wieber  and  Mickelsen  (21)  indicates  that  th--  osclllatihg  velocity,; 
effect  on  drople+  'aporization  is  essentially  nonlinear,  because  of  the 
dependence  of  the  evaporation  rate  cn  the  abccluts  magnitude  of  the  relative 
velocity.  On  the  other  hand,  the  mixing  of  the  propellants  by  the  oscillating 
velocities  is  very  strongly  dependent  on  the  orientation  of  the  liquid  spray 
fans  (produced  by  the  unlike  impinging  doublet  injector).  Although  no  de¬ 
tailed  description  of  such  a  complex  phenomenon  is  now  possible,  the  follow¬ 
ing  schematic  discussion  illustrates  one  process  by  which  the  burning  rate 
may  be  caused  to  oscillate  by  an  oscillating  transverse  velocity. 

Consider  the  mixture  of  gaseous  combustion  products,  vaporized 
propellants  (oxidizer  and  fuel),  and  liquid  propellant  droplets  at  seme 
axial  sfation  downstream  of  an  injector  doublet.  For  an  unlike  impinging 
injector,  liquid  mixing  is  imperfect  and  some  degree  of  stra+i f icati '^n  exists 
in  the  mixture.  For  concr 3+eness,  the  area  of  interest  is  sb-.wn  in  Figure 
17  (a)  for  a  doublet  spud  oriented  with  the  line  of  centers  of  its  orifices 
tangential  to  the  injection  circle.  Then  the  stratification  is  almost 
entirely  in  the  tangential  direction,  as  shown  by  the  lines  of  ernsiant 
mass  fraction  of  vaporized  oxidizer  in  the  ojS  mixture,  Y^.  The  exact 
shape  of  +he  constanf  contours  will  be  dependent  on  the  design  and  oper¬ 
ating  mixture  ratios,  the  propellant  characteristics,  etc.  Exause  of  the 
turbulence  in  the  combustion  chamber,  the  stratification  pattern  shown 
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schematical iy  in  Fi^iure  17  (a)  represen+s  only  a  mean  condition. 

A  typical  fuel  rich  droplet  is  shown  at  position  A  in  the  region 
which  is  relatively  rich  in  vaporized  fuel.  A  corresponding  oxidizer  rich 
droplet  is  shown  at  B  in  the  oxidizer  rich  region  of  the  mixiure.  As  a 
droplet  evaporates,  the  vapors  diffuse  away  and  must  mix  with  the  othgr 
vaporized  propellant  in  the  proper  pr^'pcrtions  for  chemical  racticui.  In  a 
rockel'  combustion  chamber,  the  transport  and  mixing  of  the  vapors  is  most 
likely  to  be  carried  out  primarily  by  turbulence  rather  "than  molecular 
diffusion.  The  overall  burning  rate  cf  a  fuel  rich  droplet, ,  wi 1 1 
therefore  be  a  function  of  the  amount  of  oxidizer  vapor  near  the  droplet, 

=  where  subscript  A  indicates  conditions  in  the  vicinity  of  the 

fuel  rich  droplet  at  position  A  (the  subscript  B  wi I i  be  used  similarly  to 
denote  conditions  near  the  oxidized  drop  atB). 

In  the  presence  of  small,  periodic,  transverse  gas  velocity 
oscillations  ,  the  gaseous  mixture  will  be  transported  rela¬ 

tive  to  the  droplets,  causing  oscillations  of  the  local  mass  fractions  of 
oxidizer  and  fuel.  The  perturbation  in  the  rate  can  then  be  written. 

where  the  proportionality  factor  depends  on  CtC and  or.  the  p'-opel- 

lant  physical  and  chemical  characteristics.  The  dependence  is  not  neces¬ 
sarily  that  determined  by  investigators  of  laminar  droplet  burning  (22). 

The  change  in  the  local  oxidizer  mass  fraction  due  to  the  tangential  velocity 


(2.26) 


perturbation  is 


where  is  the  relative  velocity  of  the  gases  with  respect  to  the  droplet. 


(2.27) 


-35- 


The  motion  of  +he  larger  droplets  is  small  cc  npared  to  the  gas  motion,  while 
Tp.e  sma i  I er  dropiets  fend  to  move  wi1ii  the  gas.  Therefore  the  discussion  here 
pertains  to  the  larger  droplets.  As  a  droplet  evaporates,  the  efffict  on  it 
of  the  transverse  velocity  decreases.  Thus  the  velocity  effect  described  here 
can  be  expected  tc  be  most  important  in  the  early  combustion  region,  in 
agreement  with  the  results  of  the  exploratory  roofing. 

Combining  equations  (2.26)  and  (2.27)  gives  an  expression  for  the 
rate  perturbation  for  the  fuel  rich  droplet  at  A, 


(2.28a) 


For  the  corresponding  oxidizer  rich  droplet  at  B,  the  expression  is 

iijjt 

Tr*  ”  2^,4  V  A  icj  ^  ( 

Since  the  fuel  rich  droplet  is  in  an  oxidant  deficient  region,  an  increase 
in  the  oxidizer  mass  fraction  will  increase  its  contribution  to  the  overall 
Durning  rate.  That  is,  ^  ?0  •  The  opposite  is  true  for  the  oxidizer  rich 
droplet;  <,0  since  an  increase  in  the  local  oxidizer  fraction  is  accompan¬ 
ied  by  a  decrease  in  the  local  fuel  fraction.  A  perturbation  vuioci ly  in 
the  positive  0  -direction  will  increase  the  oxidizer  rich  droplet  burniiig 
rate  and  decrease  the  fuel  rich  droplet  burning  rate.  Therefore  the  effects 
on  the  fwo  droplets  will  tend  to  cancel. 

The  net  burning  rate  perturbation  is  given  approximately  by 

i:.  -  [t  ^ ""  , 


(2.28b) 


(2.29-) 


where  the  subscript  o  denotes  conditions  s-.  the  center  of  the  spray.  Then 

>  0  as  i  c  I  I  •  That  is,  for  a  positive  velocity  perturbation, 
the  sign  of  the  net  rate  perturbation  depends  on  the  relative  magnitudes  of 
and  .  In  order  to  determine  the  effects  of  the  velocity  perturbation  on 
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The  entirs  spray.,  the  summation  must  inc'jde  ali  of  the  propellant  droplets. 

If  one  of  the  propellants  vaporires  at  a  greater  rats  than  +he 
other,  most  of  the  droplets  in  the  mixture  will  be  of  the  slower  evapor^’Ting 
propellant.  Therefore,  there  will  be  no  cancellation  of  opposing  effects  and 
the  overall  velocity  effect  will  be  determined  by  the  characteristics  of  the 
less  volatile  propellant.  For  the  liquia  oxygen-ethyl  alcohol  come  .'nor -.on 
used  in  the  experimental  phase  of  the  present  investigation,  it  can  be  esti¬ 
mated,  using  the  spray  evaporation  calculations  of  Priem  and  Heidmann  (23) 
and  the  relative  vaporisation  rate  data  of  Pass  and  Tischlsr  (24),  tha+  at 
any  station,  twice  as  much  oxidant  as  fuel  will  have  been  vaporized.  Thus 
an  important  velocity  effect  can  be  expected.  For  propellant  combinations 
using  non-cryogenic  oxidizers,  the  evaporation  rates  of  the  fuel  and  oxidizer 
should  not  be  so  different,  and  the  velocity  mixing  effect  should  be  consider¬ 
ably  reduced. 

The  analysis  of  liquid  droplet  dynamics  presented  in  Appendix  A 
shows  that  the  relative  velocity  perturbation  is  proportional  to  the 

gas  velocity  perturbation  uf*.  Then  the  net  rate  perturbation  can  be  written 
in  the  form 


=■  -M.  ur’e 


i-tut 


(2.29b) 


where 

is  a  velocity  interaction  index  analogous  to  the  pressure  interaction  index 
defined  originally  by  Crocco  (equation  2.18).  A  spray  thickness  ..  can  be 
defined  as  the  distance  between  two  reference  contours  of  constant  oxidizer 
fraction  (e.  g.  lines  I  and  2  in  Figure  17  (a)).  Then 


( 
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ar.'j  ths  velocity  Index  £.  is  seen  to  be  inverse!,  proportional  to  the  spray 
thickness.  In  general,  the  net  combustion  process  rate  perturbation  wi i i 
depend  on  both  the  radial  and  langenrial  velocities, 

iivt 


•fnvi:  '  i<f  )<2'' 


(2.29c) 


where /.r  and are  ♦'he  radial  and  the  tangential  velocity  indices. 

If  the  injector  doublet  which  produces  the  mixture  discussea  above 
is  rotated  through  an  angle  from  the  tangential  orientation,  the  resulting 
mixture  will  also  be  rofated,  and  will  appear  as  shown  in  Figure  17  (b).  The 


spray  thicknesses  in  the  radial  and  tangential  Id irectionr  are  thenr 
d  ,  qI 


—  £,  -  <5^e 


Each  velocity  index  is  Inversely  proportional  to  the  corresponding  spray 
thickness,  or 

For  the  tangential  spud  orientation,  and  ;  for  the  radial 

or ientation,^[-|^  so  that/^  -o  .  For  a  given  injector  spud,  producing  a 
spray  with  a  given  thickness  d,  the  sum  of  the  squares  of  the  velocity 
indices  wiil  be  a  constant,  independen+  of  the  angle  of  orientation. 

The  simple  argument  presented  above  is  intended  to  suggest  how  the 
transverse  velocity  oscillations  can  affect  the  mixing  of  the  'aporized 
propellants  and  produce  oscl.iations  in  the  overall  burning  ra'o.  For  the 
unlike  impinging  doublet,  it  is  seen  that  ttie  perturbation  of  the  combustion 
process  rate  is  a  linear  function  of  the  radial  and  tangential  components  of 
the  velocity  perturbation,  assuming  that  the  perturbation  is  very  anail. 
Other  injector  types  will  not  exhibit  this  lii-tar  velocity  effect;  some  wiil 
be  insensitive  to  oscillating  velocities  while  other  will  show  a  velocity 
dependence  which  cannot  be  linearized. 

At  the  present  time,  the  magnitudes  of  the  velocity  indices  and 
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their  variations  with  axial  distance  cannc  i"  be  calculated  because  of  the  lack 
of  quantitative  knowledge  of  the  processes  involved  in  liquid  bipropellant 
coniDJSfion  under  turbulent  rocket  combustion  chamber  conditions. 

In  this  chapter,  it  will  be  assumed  that  the  velocity  effects  occur 
during  the  same  rime  interval  (the  sensitive  time  lag)  as  the  pressure  effects. 
Although  this  is  the  simplest  manner  in  which  fhe  velocity  effects  can  ce 
introduced,  it  is  not  necessarily  the  proper  one.  Consideration  is  given  to 
alternative  formulations  and  their  implications  in  Chapter  III  -  .  It  should 
be  noted,  however,  that  the  alternative  burning  rate  expressions  involve  addi¬ 
tional  parameters  which,  like  ■  X  ,  cannot  at  present  be 

calculated  a  priori  from  a  consideration  of  physical  and  chemical  processes. 
Resort  to  experimental  techniques  will  be  necessary  in  order  to  determine 
the  proper  description  of  the  unsteady  combustion  process,  including  both 
pressure  and  velocity  effects.  Assuming,  then,  a  linear  dependence  of  the 
perturbation  of  the  overall  rate  of  the  processes  occurring  during  the  sensi¬ 
tive  time  lag  on  the  perturbations  of  pressure  and  radial  and  tangential 


velocity,  equation  (2.20)  becomes 


4-  e-r  <  t'  <-t 

The  time  rate  of  change  of  the  sensitive  time  lag  is  found,  by  the  same 
procedure  as  used  previously,  to  be 


(2.30) 


V  ,  -  '  ''  ^ 


-  1  - -  e  { 


-  i/' 


— _ _ _ p  L 


(2.3I) 


w;.  10^5  i  VtxrL-tl;'' 

Then  the  burning  rate  perturbation  can  be  written  as 


:  dw 


-  .  r  ■  V  „  .. . . f ,  V.  .u  ^  ;7  „r'/v.  O  ]  ) 


(2.32) 
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where 


9=  ^  I  s- 

L 


y&.c-ii) 

Ar'C^c  -f  3) 

v-rx,ro)  ^ 


ur '  ( X,;  Cti )  j 


In  general,  T=Z  Ci)  varies  from  one  propellant  element  to  another. 
The  effect  on  combustion  instability  of  such  a  spread  in  the  sensitive  time 
lag  has  been  evaluared  by  Crocco  and  Cheng  (!'.  In  this  analysis,  it  is 
assumed  that  all  propellant  elements  have  equal  mean  sensitive  time  lags, 
so  that  f-  constant. 

C.  Solution  for  the  Perturbations 

The  addition  of  velocity  terms  to  the,  burning  rate  perturbation 
expression  makes  separation  of  the  variabjss  impossible.  Therefore  a  differ¬ 
ent  approach  is  used.  The  solution  is  asuumed  to  take  the  form  of  a  series; 
the  pressure  perturbation,  for  example,  is 


i>’  ^  'Pi 


t 


(2,33) 


2 

where  ■^>,  is  assumed  to  be  0  (Ugll^),  p2  ^  ^“e  Since  the 

--  2 

governing  equations  w^jra  derived  neglecting  terms  of  0(Ug  ),  the  solution 
will  be  carried  only  as  far  as  the  first  order  correction,  .  As  dis¬ 
cussed  previously,  it  is  assumed  that  the  combustion  process  proouces  a 
through  ''elocity  in  the  combustion  chamber  which  is  small  compared  to  the 
velocity  of  sound,  and  that  the  droplet  drag,  and  hence  the  droplet  oscil¬ 
lation  amplitude,  is  also  small.  Therefore,  the  zeroth  order  solution  cor¬ 
responds  to  acoustic  oscillations,  and  the  first  order  correction  'fx;|iides 
the  effects  of  the  combustion  process,  the  liquid  gas  interaction,  and  the 


exhaust  nozzle. 
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The  zeroth  order  equations  are 


Q-i  B  r  o 

(2.34a) 

5  Uo 

-  _  1  5^0 

^  2  i 

1 

“  r 

(2.34b) 

(2.34c) 

_L 

~  “  y-r 

(2.34<i) 

V 

c 

II 

(2.34e) 

The  boundary  conditons  at  the  combustion  chamber  walls  are  that  the  velocity 
perturbation  normal  to  the  wall  must  vanishj  that  is,  u’  =  o  at  2  =  o  and 
v'  =  o  at  r  ®  I.  In  addition,  it  is  necessary  to  impose  conditions  of 
boundedness  and  periodicity  in  order  to  obtain  a  solution.  At  the  combustion 
chamber  exit,  2=2.,  the  proper  boundary  condition  is  a  relation  be- 
tween  the  perturbations  of  pressure,  velocity,  and  entropy.  (This  "nozzle 
admittance  relation"  is  discussed  more  fully  later  in  this  chapter  and  In 
Appendix  C.)  However,  the  oscillations  in  a  closed  end  cyiir.dsr  are  quits 
similar  to  those  in  a  rocket  combustion  chamber.  Therefore,  It  is  instructive 
to  consider  the  zeroth  order  solution  with  the  solid  wall  boundary  condillon 
imposed  at  the  nozzle  entrance, 

=  O 


Equations  (2.34a)  to  (2,34s)  can  be  solved  for  to  give 

=  O 


(2.35) 


where  4  is  the  Lapfacian  operator 
A 

S'-* 

The  solution  is  obtained  by  assuming 


^  r»  B9*- 

in  the  form 


(2.36a) 
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so  that  tho  other  perturbation  quantities  become 

f  (-■>  ©,  w 


U, 


J(r)  <5)/6; 

YS  dl:  ^ 


-  “ 


Pfi)  if- 


d  r 


,  i-  P(2;i 

Substitution  ot  the  c;xpression  {2.36a)  into  equation  (2.35)  yields  the 
following  three  ordinary  differential  equations: 

+  Sy^  )  -  C> 


(2.36b) 

(2.36c) 

(2.36d) 

i2.36e) 


^  *  r  _ 

^  r 


+  (4Hr-’-vV5c  =  « 


=  O 

1.  z 

where  Syi,  and  v  are  separation  constants. 

Equarion  (2.37b)  is  a  Bessel  equation  of  order  ,  for  which  tho 
solutions  are 

^  X/ 

The  boundary  condition  at  r  =  I  is  the  solid  well  condition,  :z  Q, 

At  r  =  0,  which  is  a  singular  point  for  the  differential  er.-atioji,  the 
condition  is  that  the  perturbations  must  be  finite;  I.  e.,  is  to  be 
regular.  Therefore  the  Bessel  function  of  the  second  kind, is 
excluded  because  of  its  logarithmic  singularity  at  r  =  0.  The  first 

condition  determines  the  value  of  the  constant  5Vu.  That  is,  Syt,  is  the 
solution  of  the  equation 

For  each  value  of  the  parameter  ,  there  are  an  infinite  number  of  zeros 
Syi,  ,  counted  by  the  index  h.  Physically,  h  is  the  number  of  cylindrical 


(2.37a). 

{•2v37b) 

(2.37c) 
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surfaces,  including  the  cuter  solid  wall,  on  which  the  radial  velocity 
perturbation  vanishes;  h  =  I  correspond^  to  a  purely  tangential  mode,  h  =  2 
to  the  first  radial  mode,  h  =  3  to  the  second  radial  mode,  etc. 

The  solutions  of  the  harmonic  oscillator  equation  (2.36c)  are 


iv0 


-  lye 


The  first  of  these  solutions,  »  represents  a  wave  moving  in  the  '.=gative 

e  d.rection  since 

i’o  ~  e 

u, 

and  the  phase  velocity  is  .  Simi  larlv,  tne  second  solution.  €  , 

represents  a  wave  moving  at  the  same  speed  in  the  posit iva  direction.  The 
combination  of  the  two  oppositely  moving  waves  produces  a  standing  oscillation 
pattern,  with 

'V'  cos  T^Q-e  ^  dig  ^  1/0  e 

Both  the  spinning  and  standing  solutions  must  be  periodic,  since  the  direc¬ 
tions  represented  by  0  =0  and  0  sUf  are  physically  ioentica).  Thus 


Moreover,  the  standing  mode  tangential  velocity  perturbation  vanishes  on  V 
stationary  diame+ral  surfaces  (nodal  surfaces)  which  can  th.er.  be  rspiaced  .by 
solid  surfaces  (e.  g.,  barriers  or  sector  inserts)  w^ithout  disturbing  the 
oscillation  pattern.  In  a  spinning  mode  the  nodal  surfaces  rotate  ar  fhe 
oscillation  frequency.  For  V-^O  »  ■*'ha  solution  of  equation  (2.36c)  is  simply 
(3^  =  constant 

which  represents  a  purely  radial  nede,  since  there  is  no  variation  in  the 
tangential  direction.  The  first  tangential  mode  cor.-esponds  to  V  =  1, 

the  second  tangential  mode  to  V=  2,  and  so  on. 

The  solution  of  equation  (2.36a)  which  satisifies  the  solid  wall 
boundary  condition  at  the  injector  (E  =  0  )  is 

r^(E)  =  Poo  ? 
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where  -  Pq(0) 
the  nok;zle  end  of  t 

Ai 

can  be  satisfied  if 


is  an  arbitrary  constant, 
he  combustion  chamber, 

(?-=&)  =0 

s  =  iw  .  Then 


The  solid  wall  condition  at 


p6  (Z)  =  Poo  CoS  V ^ 
and  the  boundary  condition  becomes 

Sin  -  O 

wiiich  implies  that 

or 

CO  =, 

'2e  (2.38) 

thus  the  acoustic  oscillations  in  a  closed  bylinder  are  neutral  (A.=0)  with 

discrete,  well  defined  frequencies.  The  purely  longitudinal  modes  correspond 

to  So,  =r  o,  i  •••  so  that  and 

i>o=  =  Poo  ^  =  O 

Purely  transverse  modes  are  obtained  when  •  Then  the  frequency 

of  oscillation  is  ^=rs^l,.  As  discussed  above,  Svi,  is  a  func^-ion  ot  the 
indices  V  and  which  specify  the  mode.  Isobar  and  streamline  patterns 
for  several  transverse  (standing)  mode^with  the  corresponding  values  of  Syi,, 
are  shown  in  Figure  1.  For  these  modes, 

Foo  U^Ci)  =  o 

The  most  general  three  dimensional  oscillations  Involve  combinations  of 
longitudinal  and  transverse  modes,  in  which  both  the  axial  and  transverse 
velocity  perturbations  are  of  the  same  Order  of  magnitude  as  the  pressure 
perturbation.  The  oscillation  frequency  in  the  general  case  is  given  by 
equation  (2.38)  with  both  s>i,  and having  nonzero  values. 

The  present  investigation  is  concerned  only  with  the  purely  trans¬ 
verse  modes.  By  this  restriction  of  attention,  considerable  reduction  in 
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complexi+y  is  obtained  and  the  propei  ties  peculiar  to  the  transverse  oscil¬ 
lations  may  more  easily  be  seen.  However,  important  effects  can  be  expected 
to  result  from  the  interaction  of  longitudinal  and  transverse  modes,  and 
ultimately  the  analysis  of  general  three  dimensional  oscillations  should  be 
made. 

ihe  oscillations  in  a  rocket  combustion  chamber  differ  from  the 
acoustic  vibrations  just  discussed  for  two  reasons.  The  first  is  the  pres¬ 
ence  of  the  liquid  propellant  combustion  process  with  the  resulting  small 
nonzero  through  velocity.  Second,  the  boundary  condition  at  the  chamber  exit 
imposed  by  the  exhaust  nozzle  in  supercritical  operation  is  significantly 
different  from  the  solid  wall  condition  used  in  the  acoustic  solution.  The 
nozzle  boundary  condition  for  transverse  modes  of  neutral  oscillations  was 
derived  by  Crocco;  an  analytical  discussion  of  this  derivation  is  given  in 
(25). 

Let  the  perturbations  of  pressure,  velocity,  and  entropy  be  given 

by 

a  ^ 

w'  ivd’’) 

=  5fe) 


where  or  e  are  the  transverse  distribution 

LLot 

functions  discussed  above,  and  the  time  dependence  is  harmonic,  e  .  From 
the  transverse  momentum  conservation  equations  it  follows  that  V  =  W. 


Then  the  boundary  condition  at  the  nozzle  entrance,  which  by  analogy  with 
acoustical  terminology  was  called  by  Crocco  tne  nozzle  admittance  relation 
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takes  The  form 


'P(o}  ?(o} 


+  <3B  Syl, 


Tm 


-f-e 


ys^) 

?(o) 


(2.40) 


The  admittance  coefficients,/!  ,  (B  ,  and  (2  are  complex  functions  of  the  mode, 
frequency,  chamber  cxi-r  velocity,  and  nozzle  convergent  section  geometry,  A 
discussion  of  the  calculation  of  the  admittanct  coctficients  for  an 
axisymmetric  nozzle  is  included  as  Appendix  C  of  this  dissertation. 

The  computations  which  have  been  made  indicale  that  the  admittance 
coefficients  are  app.  oximately  0(ue)>  or  slightly  larger,  for  frequencies 
in  the  vicinity  of  the  acous+ic  oscillation  frequency.  For  +'-ansverse  modes, 
and  are  0(1).  It  will  be  shown  later  that  if  combustion  is  complete 

P(o^ 

at  the  nozzle  entrance  the  en+ropy  perturbation  is  at  most  O(ue^).  Then 

from  the  admittance  relation  it  can  be  seen  that  the  axial  velocity  perturba¬ 
tion  is  not  zero  as  in  the  acoustic  case,  but  is  really  approximately 

0(Ue). 

By  definition, 

£' waw  +  o(".j 

*0  a  ^Oo 

and  from  the  discussion  above 

^  V' ■X'A.t 

That  is, 

5^^  O  ':«i  ) 

Expanding  the  funciion  and  neglecting  higher  order  terms  gives 

the  not  unexpected  -esult 

s‘4.SvA  =  OecTe^ 

since  Zq  is  6(1).  Thus,  for  purely  transvese  modes,  the  combustion 


process  and  exhaust  nozzle  give  rise  to  a  small  axial  velocity  perturbation 
and  a  small  shift  in  frequency  from  fhe  corresponding  acoustic  frequency. 

In  addition,  the  osciilatior.a  in  a  rocket  motor  may  be  unstable,  whereas 
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the  acouslic  oscillations  are  always  nei  ral. 

It  should  be  noted  that  the  series  expanbion  used  for  the  pertur¬ 
bations, 

cannot  be  applied  to  the  admittance  coefficients  ,  <S ,  and  (?.  An  expansion 
of  the  form 

Jt  ~  Jlo 

requires  that =  OCUg).  Since,yi^=  0(Ue),  the  requirements  is  that  Jk, 

W 

must  be  0( Ug"; .  The  admittance  coefficients  are  functions  or  tno  frequency 
Cd;  therefore  the  terms  in  the  series  expansion  can  be  wri .ten 

■  ("'M)  , 

Thus,  since  is  0(Ug),  the  derivative  must  he  OCUg)  also.  The 

curves  presented  in  Figure  18,  which  show  a  frequency  dependence  typical  of 
the  admittance  coefficients  in  gene,  al,  demonstrate  that  can  be  0( I ),  so 
that  the  series  expansion  breaks  down.  In  this  analysis,  therefore,  the 
application  of  the  nozzle  boundary  condition  will  be  deferred  until  the  first 
order  solution  is  obtained. 

The  equations  governing  the  first  order  solution  tor  the  pertur¬ 
bations  will  now  be  derived.  In  this  deriv^ition  the  simplifying  assurfption 
will  be  made  that  th'  gradient  of  a  quantity  is  of  ihe  samo  order  of  magni¬ 
tude  as  the  quantity  itself.  However,  an  exception  must  be  made  in  the  case 
of  the  steady  state  axial  gas  velocity,  u(E)  ,  In  order  to  allow  for  the 
case  of  a  rtrlatively  concentrated  combustion  zone,  the  veloci+y  gradient  ^ 
must  be  allowed  to  be  0(1)  locally,  ra-vier  than  0(Ue).  The  chamber  length 
2e  is  assumed  to  be  0(1);  therefore  the  axial  distance  over  which 
can  only  be  0(Ug).  The  determination  of  the  first  order  solution  is  much 
simplified  by  requir.ng  the  steady  state  velocity  gradient  to  be  of  the  same 
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order  OT  magnitude  as  the  exit  velocity.  11  is  demonstrated  in  Appendix  B 
that  the  sane  solution  results  when  the  magnitude  restriction  on  the 
velocity  gradient  is  removed,  although  ■‘‘he  analysis  is  considerably  more 
comp  I icated . 

Assumi ng, 

obtained  as  follows.  The  per turoations  are  given  by 

j  -  ' 


an,  ihaT  7^  is  0(5  ),  the  first  order  equations  are 
'  4^  e  ’  ^ 


r  = 

=  OCue  h  ) 


4r 


ys  3^  ' 


ur 


a'  '  Ai  =  oK  k) 

u,l  =  I 

Vj  .  '  r.  |±»  *  O 

Where  the  first  order  solution  has  been  found  to  be 


(2.41) 


'pc  — Fioj^l-f-  (S^+Syl,)^  J 

with  standing  .mode)  or  e 

mode),  and  S^f-Syi,- OC^J.  The  expressions  (2.4|)  are  inserted  into 
equations  (2.Ma),  (2.14b),  (2.l4d),  and  (2.14i)  and  Tie  i.0!-oth  o'-d«r 
relations  are  subtracted  off.  The  resulting  system  of  partial  differential 
equations  is: 


J  1 

.  au.,  :  9UI 
^  STg  ^  J7 

T-T  u.  j.9ir_ 

r  "rae  = 

du  j 
/i  1 

YSrTs  Y  J 

(2.42a) 

Sa, 

- 

(2.42b) 

S-v; 

-  '  ^ 
"  ~  r  §r 

/•a<- 

(2.42c) 

SW| 

=  -  -i-  ^  . 
jfiraa 

<■  5s?  + 

die  ' 
- 

^ae 

(2.42d) 

Sj£, 

Y 

-sp, 

ol.iZ  /  2(-i  ',-j 
di  v  Y 

(2.42e) 

(2.4|a) 


(spinning 


(2.42e) 
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The  boundary  v,onditions  on  thes  equations  are  those  which  were 
applied  to  the  zeroth  order  solution,  namely,  the  solid  wall  condition  at 
z  =  0  and  r  =  I,  the  condition  of  regularity  at  r  =  0,  and  the 
periodicity  condition  in  the  circumferential  direction,  including  the  same 
distinction  between  the  standing  and  the  spinning  forms  of  tangential  modes. 
The  nozzle  boundary  condition  will  be  to  the  overall  sc'ution,.  +hat 

is,  to  the  sum  of  the  zeroth  and  firsf  order  solutions. 

The  systems  of  first  order  equations  can  be  solved  for ,  giving 


dr  t'fe 


(2.43) 


.  nZ  gl  I  a  I  3i^ 

where  A .  Since  equation  (2.43)  is  linear,  the 
solution  can  be  found  by  superposition.  Let 


where  the  partial  solutions  and  •fp  are  solutions  of  the 

following  equations 

Poo^Vfc  (Si, 

*■  %  if' 

Afc-s'lp  =  [*'P.+  ]  s  f/K  ©./ 

The  partial  solution  can  be  obtained  by  letting 

0^(6) 

Then  equation  {2.45a)  becomes 


(2.44) 


{2.45a) 

(2.45b) 

(2.45c) 

(2.45d) 


(2.46) 
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The  transverse  distribution  functions  satisfy  fhe  relation 

t  ^  r  ar  ^  -  —  s>;j,  0^ 

Making  use  of  this  relation,  equation  (2.46)  becomes 

42  ^  ^  3i- 

However,  since  5^+5^  and  “  are  both  C(Me),  the  differential  equation 

“co 

for  reduces  to 


(2.47) 


(2.48) 


This  equation  is  easily  integrated  to  give 

..2 


[a  da'  t- 


(2.49) 


where  the  integration  constand  represents  the  value  of  at  the 

injector  face,  2=0;  therefore  =  0(Uq). 

The  partial  solution  -p^  is  obtained  in  the  same  manner  from 
equation  (2.45d)  with  and 


p  f^\  f'  f'’ 

ufe'Mj'  +icj, 

<»  ■'o 

where  again  \4-^®  =0('<i)- 

TOO 

A  different  procedure  must  be  adopted  in  o-dtr  io  obtain  +hs 
partial  sol utions  *1^  and because  equations  (2.45b)  ar  (2.45c)  are  not 
separable.  These  equations  are  of  the  form 


(2.5J) 


(2.51) 


The  solutions  (eigenfunctions)  of  the  corresponding  homogeneous  equa+ion  wc  a 
determined  previously  in  the  discussion  of  the  zeroth  order  perturbation 
solution.  The  transverse  distribufion  functions  for  these  solutions  were 
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found  to  be 


'P  (r)  =  r)  0  ro;  i 

^  f  s  hi  j  Wp  V  /  , 


p=o,  1,^,  3,  ...  ^  cf=  ... 

These  transverse  eigenfunctions  form  a  complete,  orthogonal  set.  A  solution 

of  equation  {?.5I)  can  be  obtained,  therefore,  by  expanding  Gc'iJ©?  in  a 

* 

series  of  the  transverse  eigenfunctions. 


OO  <30 


?=o  <1  =  1 


Write  the  solution  of  equation  (2.5!)  in  the  form 

^  oo 

a  ...  P_.  (t-)  U, 


4  =  Z  Z  Pp<, 

f*0  9** 


(2.52a) 


The  transverse  eigenfunctions  must  satisfy  the  relation  stated  previously 
in  equation  (2.49).  Then  equation  (2.5i)  yields 


^p<)  ~  CsV-S’m  J  =  AfCf  Fff  )  for  all 

Let  Clp^  =  Af(^  •  Then  the  partial  differential  equation  (2.5!)  h2"  boen 
reduced  to  the  set  of  linear,  second  order,  inhomogeneous,  ordinary 
differential  equations 

The  expansion  coefficients  are  oiven  by 

2ir  ( 

~  J  r  G ippj  0p  rcir  4$ 

6*0  '^rso 

where  H  is  a  normalizing  factor,  defined  by 

I 

I  J  @p  rJrolO 

and  the  complex  conjugate  of  (3)p. 


(2,53a) 


'2.54a) 


*  The  validity  of  the  expansion  depends  on  the  existence  of  Green's  function 
for  the  boundary  condition  ~  o  ,  which  follows  the  existence  of 

Gi-een's  function  for  the  boundary  condi  fion  •p(''^O~0’  (26,27) 


-51- 


Consider  first  +he  partial  solution  ,  given  by 


%‘ll  Bp,  Pw  W  0f 


(2.52b) 


P^o  <j=l  '  ' 

.<her-e  the  axial  distribution  functions  satisfy  the  differentia!  equations 


(2. Sib) 


The  mode  of  instaiiiity  under  consideration  is  specified  by  the  pair  of  indices 
T.’,  h.  If  p  =V  ,  q  =  h,  equation  (2.53b)  reduces  to 

\  —  /j)  ^ 

^2^1  PoJ 

exactly  as  for  the  partial  solutions  and  .  Then 

-  (X  I  d-  (2.55a) 

I  <*0 

as  before.  For  all  other  values  of  p  and  q  ,  equation  (2.53b)  becornss 

\  ( ^  _c^  ")  rp  c'fi 

fC  ;  "  p  ST? 


\  Poo  J 


since  s^r 4.  O  ^  ^ e  )  '  solution  of  this  equation  is  obtained 

by  the  method  of  variation  of  constants  as 

/O  I  )  .  I - - - 


Po,  j‘ 


S'/”?  -V„ 


|Cgp<,Cail,  Vs^-S^^  ? 


The  expansion  coefficients  Bnn  are  given  by 


5p«j®p0p'‘ , 

for  u  -pinning  mode,  Q^^^^-ipQ  .  Then,  denoting  by 

Q  ^  (svixf)  Jp  <■<?!  r 

15  po  ^  - - -f _ _ _ 

rJr  - - 

From  this  expression,  since  the  functions  £  and  e  ‘  are  orthogonal,  it 


can  be  seen  that 


_  Sy;„  J"  CSvi,  I")  J}/  C-^y^  t' 


/  •rAr 


(2.56) 
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The  £3n,e  result  is  obtained  for  the  standing  mode,  because  and 

Co5  p  0  (O$0  iai'  )  are  also  orthogonal. 

The  partial  solution  is  obtained  by  the  same  precedure  as  used 
for  .  Thus  let 


oo  fTO 

"tc  ^  X  2  ^  ^cfx? 


where  the  axiai  distribution  functions  are 


r  ^ 

^  ^  J 


fCfx»  —  I  du,^„ -SiihvSp^-SM, 

^  -i  ®  » 

The  expansion  coefficients  for  are  given  by  the  expression 

^  %'(3;  ,-drJe 

r  i'  ®p  ®*  ’^'1",''®  -.v» 

For  a  spinning  mode  Then  fiJ/s-CVe  orthogonal  to  all  of 

the  ^unctions  except  0^/*  ,  so  that 

c„-o  ,  r*-^  1 


(2.52c) 


(2.57a) 


{2.57b) 


Cvo  =-tV 


j[  Jif(^ykr)  r)  <ii 


(2.53a) 


In  particular,  for  q  =  h  ,  the  second  of  equations  (2.58a)  becomes 

p  _  -  ;  £  ^ 

TTTT; - J-  =  -ICv^ 

where  Cyi,  is  a  positive  real  number.  Thus,  the  expansion  coefficient 
corresponding  to  the  mode  of  instability  under  consideration  !s  a  negative 


(2.58b) 


purely  Imaginary  number. 

However,  for  a  standing  mode  of  osc ; I lation, ®y=CoS p6 ,  and  the 
derivative  6^^--2/S(nV^  is  orthogonal  to  ©^Wflspd  for  all  p,  over  the 
interval  Os  0  £  2ir  .  in  order  to  obtain  expressions  for  the  expansion 
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coefficients,  the  interval  definition  of  f  .e  eigenfunction  expansion  is 
changed  from  0^6<2.-7r  to  OS.0£tr  .  Since  the  standing  mode  oscillation 
pattern  is  symmetrical  about  9  =  0,  (or&^n'  ),  the  other  half  of  the 
solution  can  be  obtained  by  reflection.  Then  for  a  standing  mode. 


^0cosfeci&  f  T^.lSyi,r)Tp(S^r)  J  i 

Cv.  =  -  V  -1. _  _ _ _ 

J  co5*p0  J  JpCSf^'T) 


The  most  important  result  here  is  that  the  eigenfunction  expansion  for 
does  not  include  a  term  corresponding  to  the  mode  the  stability  of  whjch  is 
being  considered. 

It  should  also  be  noted  that  whereas  the  expansions  for  and  for 
in  the  case  of  a  spinning  mode  are  uniformly  convergent,  the  standing  mode 
series  for  , 

■i>c=  21  2  c„ 

is  not  uniformly  convergent.  This  result  follows  from  the  tact  that  the 
tangential  part  of  the  nonhomogeneous  term  of  equation  (Z.45c)^ 

-J*.  Sb 


f2.58c) 


■yf  .JinVP 


vanishes  at  0  =  0  and  Tf  •  On  the  contrary,  the  eigenfunctions 

®p  =  uo  not  vanish,  but  have  vanishing  derivatives  <3^^  at  the 

endpoints.  Thus  the  expansion 

oo 

$}n^r)i/SinV0  =  ^  ^  Cm  csJTp© 

^  T' 

does  not  converge  at  the  boundaries,  9  =  0,  St"  .  However,  the  sigenfunc-* 

tion  series  for  "fr  does  satisfy  the  tangential  boundary  condition 

’  '  ifrvTr  ' 
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and  does  represent  !ne  partial  soluti  n  ^  away  from  the  tangential 


boundaries. 


In  summary,  the  solution  for  the  pressure  perturbation  h?'’  been 


found  to  be 


1*'=  I (i>]  j;  r) 
r  fit/,  Faw,  <»;  +  Pcv^  fif)  J  Ov,  r)| 


for  the  spinning  mode,  and 


[  fpo^j  ^  4-  Po  r J 

+  ^  fti/,  y  to  vtf 

■*■  21  2?  vto  pd 

?2| 

for  the  standing  mode.  Both  of  these  solutions  are  of  the  form 

-^,'=  [i(*)  -f-  I  I 

r»tr*54 

so  that  the  velocity  perturbation^  can  be  written  as 

The  entropy  perturbation  must  be  evaluated  before,  tna  nozzle 
boundary  condition  can  be  Imposed,  T>'c.  fundamental  identical  relation  of 
thermodynamics  can  be  written  in  dimensionless  form  as 

Ti4,  a  <«.  -  ti  ^ 

f  JP 


(2.59a) 


(2.59b) 


(2.60a) 


(2.60b) 


(2.60c) 


(2.60d) 
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In  terns  of  steady  state  and  perturbation  quntities,  this  relation  becomes 

fV  =1'- 

X  f,  - 

It  has  been  shown  that  TnfnjS.  I  a„dt?,'=T/.  Then  the  entropy 
perturbation  is  given  by 

T‘  - 

or,  making  use  of  the  equation  of  state  (2.l4g)to  eliminate  the  temperature 
perturbation. 

The  perturbations  of  all  unsteady  quantities  have  been  assumed 

to  be  represented  by  series  expansions.  Thus  for  the  entropy 
/ 

^  y4'^  +/^,  •f'  •  »  • 

Where  from  equation  (2.61) 

••• 

The  zeroth  order  energy  equation  (2.34e)  shows  that 


(2,61) 


-  O 

From  the  first  order  energy  equation  (2.42e)  it  can  be 


seen  that 


or 

>  =-3:  -f  0(kc) 

i 

Thus,  the  entropy  perturbation  is  largest  in  the  region  where  the  combustion 
IS  nest  intense.  At  the  nozzle  entrance,  where  combustion  is  assuned  com> 
plete,  the  entropy  perturbation  vanishes. 

It  is  now  possible  to  examine  the  stability  behavior  of  the  rocket 


(2.62) 
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rnotor  by  applying  the  cnamber  exit  bounds  -y  condition,  which  is  the  nozzle 
admi ttanpe  relation.  The  conditions  at  the  stabi I i ty  I imi ts,  which  bcund 
the  regions  of  unstable  operation,  will  be  derived  fay  setting  the  amplin- 
caticn  factor  ?<  equal  to  zero.  Then  the  complex  time  parameter  s  is 
purely  imaginary,  s  =  i  ,  where  oJ  is  the  frequency  of  oscillation. 

E.  Characteristic  Equation 

The  boundary  condition  imposed  bv  the  exhaust -nozzle  upon  the 
oscillatory  flow  in  the  rocket  combuslion  chamber  has  been  presented  in 
equation  (2.40;  in  a  form  appropriate  to  the  zeroth  order  solution.  Because 
of  the  distortion  of  the  wave  pattern  introduced  by  the  effects  on  the 
combustion  process  of  the  transverse  velocity  fluctuations,  it  has  been 
necessary  to  make  use  .of  eigenfunction  expansions  in  the  perturbation  solu¬ 
tion,  which  has  been  found  to  take  the  following  for 

=  I  2  UfAV  (3>pC&) 

nr  I 

pco 

A  2? 

Then,  imposing  the  nozzle  boundary  conditions  results  in  not  one  relation  as 
in  equation  (2.40),  but  infinitely  many  relations  of  the  forr.i 


Too  TOO  r»&  rVa 

There  is  thus  an  admittance  relation  corresponding  to  each  term  in  the 


(2.64) 


eigenfunction  expansions. 
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Si  nc3  combustion  is  assumed  to  be  complete  at  the  entrance  to  the 


nozzle. 


Hence,  the  axial  distribution  functions  become,  using  equations  (2.59)  and 
(2.60), 

~  ^  fcvi,('2e)+ 

Vp<j  (ie_)  =  Ppaiie) 


(2.65) 


»h.re  F.,(t}=p„[n.(s>^^fj 

/  (2.£ 


I 


Pp  (i)  -  foo[^  +Kc] 

Pff,®)  =  P„  Wlf 

R  =  D  r -1 

-  Koo  ]_  J  J  ^  -  -Jnr-TT^ -  -f-  *^P9  «<«<>  VjA"^  ?  / 

"  -J 

Substituting  the  relations  (2.65)  into  the  nozz*'  i  admittance  equation  (2.64), 
and  defining  a  combined  admittance  coefficient  €.  by 

6'  e.  +  ifi  +t<8 


(2.67) 
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there  results  the  equation 

fJi 


(2.68) 


Each  of  the  partial  solutions  making  up  the  pressure  perturbation 

I  rfy'^r 

-b  includes  a  constan-j^;es‘ul  ting  from  the  integration  in  the  axial  direc- 

y 

tion.  For  each  pair  of  '/id ices  p,  q  (Thar  is,  for  each  term  in  The  eigen¬ 
function  expansion),  the  several  integration  constants  can  he  combined  into 
one.  Then  the  axial  distribution  functions  Ppq^2)  take  the  following 


form; 


Pfxf  ~  fao  ^  ^  3 


(2.69) 


where  the  function  does  nor  include  any  integration  constant  and  the 

constants  Kp^  are  ail  0(ue).  In  particular,  the  axial  function  corre¬ 
sponding  to  the  mode  of  oscillation  under  consideration  can  be  written  as 


—  Poo  {  I  } 


(2.70) 


The  function  Pjvi,  /  ich  from  equations  (2.65  and  (2.66)  is  given  by 

o 

is  0(Ug).  The  integration  constant  ICyj,  can  be  absorbed  into  the  constai  t 
as  can  be  seen  i.-'m  the  following  discussion.  Equuiion  (2.69)  is 
divided  through  by  ^  ^ )  giving 


(2o7l) 


I '  ~  ='  1+  +•  0(4; 


or 


=»  1  -¥■  ir  0(^  ) 


—  —  2 
Since  both  and  are  0(Ug  ), their  product,  is  0(Ug  ),  and  thus 

=  C  (  1+  -h  Q(^l) 


(2.70a) 
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Similarly,  replacing  Pqq  with  in  equation  (2.69)  results  in  an 
error  of  0(Ue^)»  which  in  this  analysis  has  been  assumed  to  be  negligible. 

The  solution  tor  ihe  pressure  perturbation,  equation  (2.60a)  can 
thus  be  written  ^ 

0)  ^  I  0  +  Q^(&)  + 

4-  I  T  ^  Qf(e)  | 

The  constant  ^  ,  which  represents  the  perturbation  amplitude  le\/el,  will 

not  affect  the  stability  solution.  For  each  f,*!  ,  the  nozzle  boundary 

condition 

+  Kp<jT^*^  Sioli  ^  - 

—  I  +  Kff  coslt ?e  1  =  O 

can  be  used  to  determine  the  integration  oonstant  .  Application  of  the 

remaining  boundary  condition. 


(2.72) 


(2.73a) 


(2.73b) 


results  in  an  eigenvalue  problem.  That  is,  equation  (2.73b)  is  the  character¬ 
istic  equation  for  the  eigenvalues  5=vt+L<a  •  The  following  discussion  will 
be  concerned  with  the  neutral  oscillations  occurring  at  the  stability  limits, 
for  which  • 

By  use  of  the  expression  for  obtained  from  equations  (2. 70) 

and  (2,71),  the  character  I s-ic  equation  (2.73b)  can  be  written  =s 


(P-i-i(|g)(R+(^)7=  K 


(2.74) 


where  the  right  hand  side  is  of  the  form 
U  —  -I-  ^ 

n  -  ^  h. 


h)  —  y  We  — t  ysvi,  2.  J 


(2.75a) 


(2.75b) 


(2.75c) 
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(2.76b) 


In  the  expressions  for  hj  and  ho  a  reduced  frequency  has  been 

introduced.  Since  Q(Si)  for  a  putely  transverse  mode^  the 

reduced  frequency  takes  on  values  near  unity. 

It  is  convenient  in  discussing  the  characteristic  equation  to 
define  several  auxiliary  functions: 

a.  Sf-  1 

'^=  J 

The  functions  hj  and  h2  in  equation  (2.76)  can  be  vritten  in  terms  of 
thest'  auxilia'-y  functions  as 

K ( i-  ) 

As  written  in  equation  (2.74)  above, the  characteristic  equation 
exhibits  the  balance  which  exists  at  the  stabllliy  limits.  In  order  for 
neutral  oscillations  to  be  maintained,  the  energy  supplied  by  the  combustion 
process,  represented  by  the  left  side  of  equation  (2.74),  must  be  completely 
absorbed  by  the  fluid  mechanical  processes,  which  are  represented  by  the 
function  h  on  the  right  side  of  the  characteristic  equation.  If  more  Than 
sufficient  energy  is  supplied,  the  oscillation  amplitude  will  grow.  Such 
an  csci'lation  is  li,jp?-|y  unstable.  Eventually,  nonlir.ut-r  energy-absorbing 
P'-ocesses  will  restore  the  energy  balance  at  some  high  level  of  amplitude. 

The  importance  of  the  oscillation  ro  the  proper  operation  of  the  rocket  motor 
depends  upon  the  amplitude  level  finally  attained.  However,  ii  has  been 
observed  that  nearly  all  transverse  mo>.e  instabilities  grow  to  very  high 
levels,  causing  rapid  combustion  chamber  destruction.  Stable  operation  results 
in  the  case  that  the  energy  release  is  insufficient  to  balance  ihe  energy 
absorbed  by  the  osciliation.  The  fact  that  the  characteristic  equation  is 
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complex  implies  that  there  is  a  phase  condition,  as  well  as  the  amplitude 
condition,  relating  the  oscillating  rates  of  energy  release  and  energy 
absorption  for  neutral  osci i lations. 

Assuming  that  the  pressure  ana  veirvcity  effects  are  character  i zed 
by  the  same  mean  sensitive  time  lag  'C  •  the  definitions  of  t‘  burning  rate 
parameters  are  .  _ 

(S’ =  ^  (I 

ce  6-9-'''^) 

^  ,  -ttot', 

J  -  e  ) 

Then  the  characteristic  equation  takes  the  form 

[l  +  ‘  ^  +  I (2.78) 


This  complex  equation  is  equivalent  to  two  real  equations.  Therefore,  for 
neutral  oscillations  of  a  given  rocket  motor.  The  rour  combustion  parameters, 

independent.  If  any  three  are  known,  then 
the  fourth,  together  with  the  osci Nation  frequency,  can  be  determined  from 
equation  (2.78).  Because  the  frequency  dependence  of  h  is  quite  com¬ 

plicated,  it  is  most  convenient  to  regard  the  frequency  as  one  of  the  indepen¬ 
dent  variables.  In  order  that  the  present  investigation  may  be  related  more 
easily  to  previous  studies  of  both  longitudinal  and  transverss  modes,  the  two 
dependent  variables  will  be  taken  to  be  The  mean  sensitive  tim^s  lag  V  and 


the  pressure  interaction  index  /)u  .  Thus, 


w' 1 1  be  determined  from  the  characteristic  e-iuation,  which  can  be  written 


At  presen"^  none  of  the  combustion  parameters  can  be  calculated 


(2.79) 
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a  priori,  although  a  method  has  been  developed  for  measuring  and  H 
experimentally  (16).  Since  very  little  is  known  about  the  velocity  parameters 
their  effect  will  be  shown  by  calculating  /n,  and  f  for  assumed  reai  values 
of  the  ratios  and  .  The  modifications  introduced  by  the  assumed 

velocity  effecis  will  then  be  noted  and  compared  with  the  results  of 
appropriaie  experiments. 

The  real  and  imaginary  parts  of  equation  (2.79), 

cos  «r)=:  K’,.  ■) 

r  (  <2. 

/n,  sth  uyv  =  nt  j 

comprise  a  simultaneous  pair  of  equations  form,  and  T  .  The  solution 


(2.80) 


is  found  to  be 


_  W  +  K; 


(2.81) 


where  V  is  determinea  modulo  ‘iT'  ,  When  uV—lf  ,  S/n  tdT*  ■  O  and 
(  I  -  CoS  wT  ) «  2  •  Then 

/n  =  tir  =  ) 

For  physically  realizable  systems,  >  o  »  ^ 

It  can  be  seen  that  /n^.,  is  the  minimum  value  of  t'^e  'nterartion  index,  for 

/V/ 

when  ,  K4o  and 

f 

^  •  I'  -t-J- 

which  is  greater  than  since  Kr  and  hi  are  positive.  Fo" 
there  are  two  solutions  for  f  ,  name'y. 

r,  =a  s!.-'  (-  ^  A  ) 

with  Ti,  >  r„  >  X,  and  <  u>^  ^  . 


(2.82) 
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A  typical  solution  for  ai(w)  and  for  assumed  values  of 

■^7^  and  is  shown  in  Figure  19.  This  solution  applies  at  the 

stability  limits^  where  X  0.  It  can  be  seen  that  for  any  given  value 
of  the  time  lag,  only  one  value  of  the  pressure  interaction  index  is 
consistent  with  neutral  oscillations.  For  the  same  'C'  ,  'j  larger  value 

of  /)v  corresponds  to  instability  (X>o  ),  a  smaller  value  to  stability 
{  X-<o  ).  Thus  the  plane  is  divided  by  the  loci  of  0  Into 

stable  and  unstable  regions,  as  shown  in  Figure  20. 

A  given  rockc-t  motor  system  at  given  operating  conditions  Is 
assumed  to  be  characterized  by  a  set  of  values  of  the  combustion  parameters 
and  .  It  the  point  on  the  /rt',T-plane  representing  this 

motor  falls  within  one  of  the  unstable  regions  (i.  e.,  above  the  stability 
limit  curve)  corresponding  to  one  of  the  transverse  niodes  and  to  the 
appropriate  values  of  and  ,  then  the  operation  of  the  rocket  will 
be  unstable,  and  large  amplitude  oscillations  can  be  expected  to  result. 
Therefore,  any  effect  which  increases  the  size  of  the  unstable  areas  on  the 
-UjT'plane  (e.g.,  by  shifting  the  stability  I imi fs  to  smaller  ^  fcr  all 
T,  or  by  flattening  the  limit  curves  without  changing  >rvf^  is  a  destabi¬ 
lizing  effect. 

Since  the  function  in  equation  (2,79)  depends  on  s>>>,  , 

There  Is  a  different  stability  pattern  for  each  inode  of  oscillation.  In 
Figure  21,  stability  limits  are  presented  fcr  several  purely  transverse 
modes.  For  a  given  interaction  index  value,  the  value  of  the  time  lag  at 
the  stability  limit  is  nearly  inversely  proportional  to  the  acoustic  frequency 
.  Thus  higher  modes  are  associated  with  smaller  sensitive  time  lags. 

By  definition. 


Hence,  a  shorter  time  lag  is  obtained  by  a  decrease  in  the  stagnation 
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sound  velocity  of  the  combustion  gases,  by  an  increase  in  the  combustion 
chamber  diameter,  or  by  a  decrease  in  th®  dimensional  time  lag  itseif. 

The  nature  of  the  functional  relationships  involved  in  the 
characteristic  equation  precludes  a  general  analytical  solution.  Therefore, 
it  is  necessary  +o  rely  on  numerieal  solutions  for  specific  cases.  A  large 
volume  of  calculations  is  thus  requiree  to  explore  the  dependence  of  rhe 
stability  patterns  on  each  physical  parameter.  In  this  situation,  many 
important  insights  can  be  gained  by  considering  a  very  much  simplified 
rocket  motor  “'ystem,  for  which  some  analytical  results  are  also  possioie. 

The  simplest  steady  state  combustion  distribution  is  that  in  which 
all  of  the  combustion  occurs  at  a  single  axial  station.  This  ideal  combus¬ 
tion  pattern  Is  actually  a  good  approximation  to  that  produced  by  an  efficient 
injector  (such  as  the  unlike  impinging  doublet  type  used  in  the  experimental 
phase  of  this  investigation),  which  provides  good  atomization  and  rapid 
mixing  of  the  liquid  propellants.  The  steady  state  axial  velocity  distri¬ 
bution  for  concentrated  combustion  is  given  by 


0 

II 

0  i  Zm 

^  -ie 

where  Zm  is  the  location  of  the  combustion  front  and  Zq  is  the  nozzle 
entrance  station.  The  other  simplification  to  be  introduced  Is 

U.^  ^  'H-i-k.  —  ConS'tO-rff 

which  is  applicable  when  the  droplet  momentum  coefficient  -fe.  is  not  too 
large. 


For  this  simple  model,  the  auxiliary  functions  defined  in 
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equation  (2,75a)  become 
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Since  the  frequency  corresDond'"^  to  the  minimum  interaction  index  is  close 
to  the  acousric  frequency,  the  reduced  frequency  ^  can  be  approximated 
by  t-i-f  ,  where  I 

Consider  first  the  case  in  which  there  are  no  velocity  effects 
on  combustion,  only  pressure  effects.  Then  ^  and  v.hcrc 

hp  and  hj  are  the  real  and  i.maginary  parts  of  the  f-.n-tion  h  defined 
in  equation  (2.75)  as 
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Introduc i r.g  the  expressions  (2.83)  into  nc  definitions  of  h,  and  h2 
of  equation  (2.75)  gives 

=  ytTe  (  ) 

h,i  =  — y  «^/c€r 

These  expressions  can  be  simplified  by  noting  the  order  of  magnitude  of  each 
term  and  nP2*''c  fi  "'•g  +he  hicher  order  terms.  It  can  be  seen  from  equation 
(2.83)  that,  since  and  are  QCM0)v<h]\e  and  ?-*€  b(i), 

Q.I- r  > /«•  0(0 

The  real  part  of  tne  combined  admittance  coefficient,  ,  is  at  most  0(ae) 

(however,  the  imaginary  part  can  take  on  somewhat  greater  values,  so  that  6-^ 
is  0(i)  in  the  range  of  exit  velocities  considered  here. 

Straightforward  calculation,  making  use  of  the  above  order  of 
magnitude  observations,  yields 

J 


^ _ -  f  f-r 

2yo*Afj  I  J 

where  from  the  condition  V?;  =0  the  frequency  i '■  r!>'en  by 

/ 

'f' - r~Tjr 

The  important  quantities  governing  the  stability  I imi  rs  are  the 
chamber  exit  velocity  u„,  through  the  semi ‘■ta nee  coefficients  Sy- and 
(which  are  also  fmquericy  dependent),  the  droplet  drag  factor  ^  ,  and  the 
length  parameter-;  .  The  quaniity  ^  ,  which  measures  the  volume 

occupied  completely  by  gaseous  products  of  combustion,  has  a  very  small 


(2.86) 


(2.87) 
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infiuence  on  .  The  derivative 

^  'i-i'  f  li-/:?;^)^ 

is  at  most  OCUg).  The  function  (jj.  ,  which  is  given  approximately  by 


oscillates  very  rapidly  with  respect  to  2^  (in  general ,  is  at  least 

35).  Thus,  +he  small  spreading  of  the  combustion  front  which  must  be 
present  in  any  ac a: i  rocket  chamber  will  cause  the  effect  of  5"  to  be 
cancel  led. 


Figure  22  shows  the  variation  of  the  nozzle  admittance  coefficients 
with  the  combustion  chamber  exit  velocity  for  two  frequencies,  (Vs  Svk 
and  *>  =  1.05  Xvii  .  1^°'’  longitudinal  modes,  the  nozzle  admittance,  defined 

by 


yyjie) 
F  Ote) 


is  positive.  It  has  been  shown  (I)  that  the  nozzle  exerts  a  strong  damping 
effect.  On  the  other  hand,  for  transverse  modes,  the  admittance  is,  from 
equation  (2.68), 
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r  ^  ^  C^) 

^  Pi^e) 


so  thar 
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For  all  but  the  smallest  nozzle  entrance  velocities,  the  transverse  mode 
nozzle  admittance  ratio  is  positive.  Therefore,  in  transverse  modes,  the 
exhaust  nozzle  exerts  negative  damping,  i.e.  it  has  a  destabilizing  effect. 

At  the  lower  velocities  (u„C,  0. 1C),  £i  is  nearly  constant, 
but  %  increases  rapidly  with  increasing  velocity.  Therefore,  the 
minimum  interaction  index  shows  a  sharp  decrease  with  increasing  exit 
velocity  in  this  range,  as  can  be  seen  from  equation  (2.86).  Figure  23 
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shows  ths+  The  eniire  stabili+y  limit  curve  is  shifred  downward  when  the 
exit  velocity  is  increased  trom  0.05  to  0.10.  For  larger  velocities, 
levels  off  and  +hen  decreases,  bur  the  stabilizing  eftec+  of  the  ^  decrease 
is  otfssf  by  the  destabilizing  effect  of  the  decrease  of  2j, .  The  resulr 
is  th  /tv continues  to  decrease  with  increasing  ijg,  al+hough  muci  more 
slowly.  Some  combustion  chamber  configurations  show  in  increase  in 
with  velocity  for  ii*  y  0.30.  However,  these  resufts  cannoT  be  trusted 
because  rhe  0{w^#  terms  which  have  been  neglected  may  be  very  important- 
Hence,  in  the  range  of  exit  velocities  appropriare  to  the  small  Mach  number 
assumption,  increasing  rhe  combustion  chamber  exit  velocity  (decreasing  the 
nozzle  cc<T f  action  ratio)  always  has  a  destabilizing  effect. 

rne  nozzle  admirtonce  coefficients  shown  in  Figure  22  correspond 
to  a  nozzle  v/ith  a  linear  velocity  profile  in  rne  subsonic  section.  For  a 
nozzle  with  a  conical  converging  part,  the  coefficients  show  similar  trends, 
but  the  values  of  ^  for  tho  "conical"  nozzle  are  much  larger  than  those 
for  the  "linear"  nozzle,  ecpacially  the  range  C(«  <  0,15.  For  higher 
exit  velocities,  the  coefficienfs  for  the  two  types  of  nczrie  tend  to  con¬ 
verge.  Therefore,  the  minimum  interaction  index  is  smaller  for  the  conical 
nozzle  than  for  the  linear  for  corresponding  combustion  chamber  condi cions, 
as  shown  in  Figure  24. 

The  dropler  Oiag  coefficient  |5  is,  by  equation  (2.83),  the  product 
o*  tne  moementum  interchange  coefficient  ,  which  represents  The  rate  af 
which  energy  is  transferred  in+o  the  oscillatory  motion  of  the  liquid 
propellant  droplets,  and  the  droplet  I  i  fetime,  ,  that  is,  the  time 

interval  during  which  the  droplets  can  absorb  energy  from  tfa  oscillations 
of  the  gases  in  the  combusf ion  chamber.  From  equation  (2.86), 
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which  is  positive  and  0(1).  Thus,  increasing  ^  has  an  appreciable 
stabilizing  effect,  as  shown  in  Figure  25.  The  quantity  jS  can  also  be 
considered  as  the  ratio  of  the  distance  of  the  combustion  front  from  the 
injector  face,  Z^,  to  the  distance  of  droplet  penetration  into  stagnant 
combustion  gas,  Zp  =  ,  based  on  a  mean  droplet  velocity  and  an 

effj^rtive  momentum  coefficient.  Then  jS  can  be  seen  to  depend  on  the 
injector  spray  characteristics,  including  the  droplet  size  distribution 
and  the  effectivenesc  of  liquid  mixing,  and  on  the  droplet  vaporization 
characteristics.  The  ^  range  for  these  calculations  (0.2^  <  0.7)  was 

chosen  so  as  to  include  values  which  are  likely  to  be  found  in  actual  rocket 
motors.  However,  further  basic  information  is  needed  to  establjish  the  value 
of  ^  more  definitely. 

The  effect  of  the  chamber  length  parameter,  ^  9^  •  Zg,  can 

be  seen  from  equation  (2. 36)  by  calculating  the  derivative. 


_ 


For  most  cases  the  positive  terms  dominate,  so  that  an  increase  in  combustio" 
chamber  length  has  a  stabilizing  effect.  The  only  exceptions  are  cases  witli 
low  jS  values  and  very  low  exir  velocities  (ug'^O.Oo),  lor  which  is 

negative  (linear  nozzle).  However,  is  onfy  OCug)  in  any  case,  so  that  a 

large  change  in  length  is  required  to  p.-oduce  a  significant  shift  in  the 
stability  limits.  In  Figure  26,  stability  limits  are  shown  for  the  ♦irst 
tangential  mode  for  two  lengths  which  are  in  •‘i.e  ratio  of  2  to  I.  Because  of 
the  factor  Svw  »  changes  in  length  are  more  significant  for  higher  transverse 
modes  than  tor  the  first  tangential.  If  the  length  is  increased  beyond  the 
critical  length,  longitudinal  and  mixed  longitudinal -transverse  modes  will 
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also  be  unstable,  ar-i  the  present  tnec  y  will  not  be  applicable. 

The  preceding  discussion  has  been  concerned  only  with  pressure 
effects.  If  now  the  tangenfiai  (but  not  radial)  velocity  fluctuations  are 
also  assumed  to  affect  the  combustion  process,  the  character i stic  equation 
becomes 


/  _  C  n>  T  .  - 


The  ratio  of  the  minimum  interaction  index  including  tangential  veloci+y 
ef fects,''i^^  ,  +o  thar  for  pressure  ef.'cCTS  aione,<ri,^p  ,  is  Then  £ffT- 

^MT  _ 

i 

fherefore,  since  >  0,  a  positive  value  of  the  tangential  velocity 
interaction  index  causes  a  downward  shift  of  the  stability  limit  curve  on 
the  plane  by  contracting  the  ordinates  of  the  limit  points,  which 

is  a  destabilizing  rrend.  Since  ^  is  real,  there  is  no  effect  on  the  tim|^^7; 
lag  value  at  the  minimum  point.  Figure  27  shows  the  Influence  o+  the  5=^5" 
tangential  velocity  effect  on  the  stability  limits  for  the  first  tangentia 1 
spinning  mode.  For  other  modes  of  oscillation,  the  effects  w! ! !  be  qual- 
itatively  the  same.  The  magnitude  of  the  velocity  effect,  for  a  given 

y»v  -t./' 

value,  depends  on  the  ratio  of  the  expansion  coefficient  to  the 

«  ■’_’V 

frequency  .  This  ratio  is  approximately  et^uai  to  ,  values 

•Jvt. 

of  which  are  given  in  Table  II  for  several  modes.  The  second  tangential 
mode  is  seen  to  be  slightly  more  influenced  by  the  iangential  velocity 
osc illations  than  the  other  modes.  It  should  be  recalled  that  the  expansion./  ^ 
coefficient  Cy^,  vanishes  for  a  standing  tangential  mode,  so  that  there 
are  no  tangential  velocHy  effects  for  such  an  oscillation. 

The  velocity  index  JIq^  is  positive  in  the  case  tnar  the  burning  ^ 
rate  is  increased  by  a  positive  velocity.  For  the  tangential  velocity,  the 


positive  direction  is  not  fi-^ed  vvirh  respect  +o  the  rocket  combustion  chamber, 
but  is  the  direction  in  which  a  spinning  wave  travels.  Assuming  that  the 
injector  is  characterized  by  a  positive  direction  (determined  by  the 
orientation  of  the  injection  pattern),  the  velocity  effect  of  a  wavs  moving 
in  the  positive  directiori  will  augment  the  pressure  effec"^.  However,  lor  a 
wave  moving  "n  the  o+her  direction,  the  o-^fect  will  opnose  +hal  of 

the  pressure.  It  is  apparent,  tnerefore,  that  the  preferred  wave  travel 
direcfion  with  respect  fo  the  physical  frame  of  reference  is  the  same  as  the 
positive  dir&'.ticn  of  ins  injector.  •^hc.ut  more  detailed  knowledge  of 

the  combustion  process  than  is  presen+ly  availaole,  this  direction  cannot  be 
determined  a  priori,  but  must  be  establisned  by  experimental  evidence. 

Consider  now  the  effect  of  the  radial  velocity  coupled  with  the 
pressure,  i.  e.,  ^  0  but  0.  Then  the  characteristic  equation 

i  s 
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Since  the  ratio  is  complex,  one  can  expect  an  influence  on 

well  as  on  ^  .  At  .'he  minimum  pressure  interaction  index  point, 


rr^ 


as 


-  £7 


Thus,  wnen  the  radial  velocity  effect  is  included,  the  condition  at  the 
minimum  is 


H,=:( 


3ri*i 


rather  than  h. 


0.  Then  =  ity 


or 


~z 


,  as  before. 


Prom  Figure  28  it  can  be  seen  that  a  positive  radial  velccity  index  implies 
a  higher  frequency  at  any  'Vv  than  that  corresponding  to  pressure  effects 
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dlone.  The  flahness  of  The  curve  indicates  that  the  iiuniinum  value  of 

the  pressure  inaex  will  not  be  significantly  altered  by  the  addition  of 
radial  velocity  effects*  Since  ~  ,  a  positive  value  of  j?,-  corre¬ 

sponds  to  a  shift  Q<-  fhe  trinimum  point  toward  siiialler  'C  .  However,  the 
(nagnitude  or  the  shift  is  /ery  small,  as  a  result  of  the  stroncj- variotfon  of 
Inr  with  Ul  and  the  low  vali!;;'^  +?.K“n  bv  the  ratio  ,  listed  in 

Table  II.  For  the  sh'ft  of  the  stabiilry  limits  reverses,  because 

a  positive  radial  velocity  index  brings  about  an  increase  in  uiV  which  is 
slightly  y! eater  than  the  increase  in  .  Thus,  except  fc;  the  near 

vicinit")  of  the  minimum  point,  the  limit  curve  is  shifted  in  the  direction 
cf  longer  time  lags  for  a  positive  index  •  These  results  are  shown  in 

Figure  29.  In  the  cylindrical  coordinate  system,  the  posirive  radial  direc¬ 
tion  is  fronri  the  center  of  the  chamber  toward  the  outer  wail.  If  a  positive 
velocity  increases  the  burning  rate,  is  positive  and  the  stability 

limit  curve  is  shifted  toward  longer  "C  ,  and  vice  versa.  However,  ft  can 
be  seen  that  the  magnitude  of  the  radial  velocity  effect  is  much  smaller 
than  that  of  the  tangential  velocity  effect  for  equal  velocity  interaction 


indices. 
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CHAPTER  III 

Extensions  oi  Instability  Analysis 


A,  Instability  Analysis  with  Nonun' form  Injection 

In  the  preceding  analysis  it  was  assumed  that  the  stead'^  state  flow  in 
the  combustion  chamber  was  one  dimensional.  That  is,  the  injection  of  the  liquid 
propellants,  their  conversion  to  combustion  products,  and  all  fluid  properties 
were  taken  +o  be  uniformly  distributed  across  each  transverse  section.  The 
anaiysis  of  longitudinal  modes  of  instability  developed  in  (I)  has  demonstrated 
that  the  axial  distribution  of  the  combustion  relative  to  the  pressure  pertur¬ 
bation  distribution  is  of  great  importance  in  determining  the  stability  behavior 
of  the  rocket  motor.  According  to  the  combustion  mode!  used  in  that  anaiysis, 
the  fluctuating  energy  release  is  proportional  +o  the  amplitude  of  the  pressure 
fluctuation.  At  a  pressure  node,  the  pressure  oscillation  vanishes.  Therefore, 
if  all  of  the  combustion  is  concentrated  at  an  axial  station  corresponding  to  a 
node  for  a  given  longitui  , al  mode,  no  oscillation  in  the  burning  rate  will  result 
and  the  rocket  motor  operation  will  be  stable  with  respcc+  to  tl.e  c! /en  mode. 

For  a  combustion  front  located  a  short  distance  from  the  nodal  point,  the  pres¬ 
sure  perturbation  amplitude  is  not  zero,  and  instability  is  possible  if  the 
sensitivity  of  the  combustion  process  to  pressure  variations  (measured  by  the 
interaction  index)  is  sufficiently  large.  As  the  combustion  front  is  moved  away 
from  the  node,  the  magnitude  of  the  interaction  index  required  for  instability 
decreases,  reaching  a  minimum  value  when  the  combustion  front  is  located  at  the 
pressure  anti  node.  If  the  combustion  is  not  concentrated,  but  is  spread  over  a 
certain  axial  distance,  then  the  stability  behavior  depends  on  the  mean 
sensitivity  over  the  combustion  zone. 


ror 


pui  c  I  y 


irai(^vc;i  9^  v/i 


instability,  the  most  important 


variations  in  the  oscillating  quantities  occur  in  the  lateral  rather  than  in  the 
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axiai  direction.  Therefore,  the  transve  se  distribution  of  the  combustion 
relative  to  the  transverse  pertubation  di si ri butions  is  of  great  significance. 

Of  necessity,  some  nonuniformity  in  +he  injection  distribution  is  presecr  in 
any  liquid  propellant  rocket  motor.  Although  the  flow  in  a  rocket  vhrix^t 
chamber  is  highly  turbulent,  there  is  definite  evioenr-e  of  strati ficaiion  of 
the  combustion,  which  persists  througn^ut  irs  flow,  t-ven  past  Mie  ,c  "ihroat 
This  stratification,  observed  optically  using  both  rectangular  (26)  and 
circular  sector  motors  (9)  wirh  transparent  l•.•aiis,  lias  neen  siearly  identified 
with  a  nons.n :  propellant  injection  di  str  i  butior  , 

The  exact  treatment  of  unsteady,  axi symmetric,  two  phase  turbulent 
flow,  with  combustion  distributed  arbitrarii,'  in  three  dimemsions,  is  quite  a 
formidable  proplem.  The  steady  stai.e  solution  appears  to  cf+er  the  groatest 
difficulties,  since  the  solution  for  ths  perturbations  can  be  based  on  an 
expansion  method  similar  to  tnat  used  in  Chanter  11.,  Little  quantiiative 
data  is  presently  aval laole  concerning  the  proce, 'es  which  are  important  in 
rocket  motoi  combustion  or  the  three  dirmensional  distribution  of  combustion 
wliich  exists  in  a  typical  iro*or.  Thus  't  is  likely  that  the  exact  oolutic.n 
of  the  general  problem  wi I }  be  difficult  to  interpret  in  terns  uf  actual  rnotor 
desi gns. 

Therefore,  because  of  the  problems  associated  '-'ith  the  excel  sol'  tion 
and  in  view  of  the  observed  strati freat ion  of  fombust i or' ,  Ihe  following  approxi¬ 
mate  approach  is  used  in  the  present  analysis.  The  one  dimensional  rrature  of 
the  steady  state  flow  is  retained.  The  effect  of  the  nonuniform  transverse 
distribution  of  cemtusTion  appears  only  as  a  variation  in  th?  sensitivity  of 
the  ccmbuc+;on  process  to  the  effects  o*  pr^s..siire  and  velocity.  This  variaticn 
is  described  by  a  distribution  functic,n  (r,4'  ).  Thus  the  burning  rate 
perturbation  for  nonuniform  inject'on  is  relcted  to  that  for  uniTorm  injection  by 


(3.1) 


The  distribution  function  is  identified  with  the  .njection  density  0'(v^  p-j  , 
defined  by 

I  'in' 
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where  /»ri^  is  the  total  mass  inject  ion  rat'^  Ther:  can  be  defined  as 

the  ratio  of  the  local  injection  density  ^  to  the  mean  injectioii  density  S'yrx. 
where 

J  TV  -  I  S(rje)  rJrdS 

^  ^  O  o 

That  is,  the  burning  rate  distribution  function  is  defined  by 
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with  the  normalizing  condition 
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For  uniform  i n  ject  ion,  ^  (r j  #  j  =.  j  , 

It  is  clear  that  this  approach  invalves  a  contradiction  in  that 
the  steaay  state  solution  assumes  that  all  fluid  properties  are  uniform  across 
every  orct-cr,  ■■■herccs  the  pei  i -•-bation  solution  is  based  on  a  no  >.niform 
tra''sversc,  combustion  di  srribution.  However,  it  was  found  previously  that  the 
major  part  of  the  oscillation  corresponds  to  that  of  a  homogeneous  gas,  because 
the  Mach  number  ir.  the  chamber  is  small  compared  to  unity.  The  presence  of 
the  co.Tibustion  introduces  a  higher  oraer  carre.,.ion  which  is  primarily  due  to 
the  sensitivity  of  the  combustion  process  to  fluctuations  of  oressure  and 
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The  identification  of  the  burning  rate  distribution^  with  the  in¬ 
jection  density  ff  neglects  the  spreading  of  the  streams  of  injected  propellants 
as  they  are  atomized,  vaporized,  and  undergo  turoulent  mixing  with  the  hot 
combustion  products  prior  to  the  final  chemioal  reaction.  Unfortunately,  it 
is  not  possible  with  present  knowledge  to  calculate  the  spreading  of  the 
propellant  streams.  However,  the  effect  of  spreading  on  stability  can  be 
estimated  by  introducing  appropriate  burning  rate  distribution  functions  Into 
tne  analysis. 

Acccraingly,  the  steady  state  solution  for  the  case  of  nonuni  form ly 
distributed  combustion  is  taked  to  be  that  given  by  equations  (2.9)  through 
(2.12).  The  zeroth  order  perturbation  solution,  equation  (2.4Ia;,  is  also  un¬ 
changed,  since  it  does  not  involve  the  combustion  at  all.  However,  the 
differential  equation  for  the  first  order  pressure  perturbation  becomes 

The  boundary  conditions  are  unchanged  by  the  i nt.^oduct i on  ot  the  hurpino  rate 
distribution.  Therefore,  as  before,  the  solution  is  obtained  by  superoosition, 
so  that 

'Pr  K  •*'■1  -*■%  (3.7) 


where 

At*  -i*  f,  Poo  iPvAJ 


(3.oc) 
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s-  pc  -  J  ^  ~  /^(■'<'j&) 


A C 1^-^ 'p^  ;- j -S  /'/y  ^ ) 

(3. 

Sines  squation  (3,Sa)  is  iaenticai  to  equation  {2.45d) ,  the  solution 
is  that  given  by  equation  (2.50). 

Secause  of  the  presence  of  the  aistribution  function  /4  ( <;& ), 
it  is  necessary  to  upply  the  eigenfunction  expansion  methoo  to  as  xell 

33  tc  -^3  2hd  .  Thus,  let 


8d) 


=  £  2  Q  c&) 

p=o 


(3,Sa) 


Vila 

‘P  ^ 


(3.9b) 


9* 


Pio  ?*•• 


{3.9c) 


Kners 


A  _ 


{ 

) 

4 


1  I  ^'br 


0^(2^*v-dr4S 


-r‘f 


.n>  A;®)  £>, 


8 


r  ro^ 


!r  6^0  ^ij,r 

31  3fl  '-'P 


(3.i0e 


(S.IOb) 
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C  = 


^  I 

J  j i//,  ^n  (Sl'^^drde 

O  o  _ 

r  j*  ®'  ®*  '■“' 


T^'e  (inferential  equation  to  be  satisfied  by  the  is 

^  ^  Oq 

wltich  has  the  solutions 

^  i  -  Ys^  I  U.fe')el2-'  + 

To.  ‘J* 


^  _ 


(3.10c) 


(3.11) 


(3.12a) 


■  l^Af^  CK^  (3_  125) 


The  solutions  for  Pe^  and  are  of  this  same  form,  and  are  given  by 
equations  (2.55)  and  2.57).  Thus  the  pressure  perturbation  becomes 

-i'  r  [Poii)  +  +6vfc  PbvJ’-)  +C*  (3. 13) 


■*■  C^6)»,c*>  t-  Psp<fa) 

Equation  (3.13)  is  of  the  same  form  as  equation  (2.60a),  so  that  the  veior  ity 
perturbations  are  given  by  equations  (2.60b),  (2.60c),  and  (2.608)  i.-.  this  case 
also. 

The  only  modifications  Introduced  by  the  nonuniformity  of  the  Injection 
and  combustion  are  the  Introduction  of  the  expansion  coefficients  snd  the 
-edeflnlng  of  and  to  include  the  distribution  function  r]  ©  ). 
Therefore,  the  characteristic  equation  takes  the  form 


in  which  ihe  function  is  unchanged  from  the  u.niform  case  and  is  given  by 
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the  expressions  of  equation  (2.75). 

The  effect  of  the  nonunf formities  can  be  seen  immediately  by  consider¬ 
ing  the  case  in  which  only  pressure  variations  affect  the  burning  rate.  Then 


(  (-  €  )  -  T 


(3.15) 


The  coefficient  >L  is  always  real,  so  that -only  the  value  of  the  interaction 
index  is  changed.  Thus,  for  ail  w  , 

/n.  I 

and 


(3.1 6a) 


t:  r. 


(3.16b) 


For  a  combustion  distribution  c)»arar.terlzed  by  Ayi,  >  r:^iv-r. .. 

That  is,  the  ordinates  of  stabUlty  limit  curve  on  the plane 
are  contracted,  increasing  the  area  of  the  unstable  region  ebove  the  curve. 
Therefore,  an  injection  distribution  for  which  Is  TOre  unstable 

than  a  uniform  one,  for  corresponding  mean  combustion  chamber  conditio.ns. 

OtjfnsiTler  tlr^  the  case  in  which  the  only  nonuniformlt*/  is  >h  the 
radial  direction,  so  that  ^  ■* yu.(ir).  This  case  is  of  the  greatest  Interest 
in  this  investigation  Because  in  the  rocket  motors  used  for  the  experimental 
phase  the  injection  of  the  propel !ants-W^  concentrated  at  a  single  radius, 
fhus  providing  a  highly  hohuniform  distributloh'.  Then 


'Vh 


The  integral  in  the  denominator  gives 

.1 


S r)  r  Jr  ^  1  Xl,  >  1 , 


-SO¬ 


SO  that  the  pressure  coefficient  A  .  b -'Comes 

fi  ^ 

^  ~  ^ i — -  ■—  -■ 

2  V^I,,  CSyh)  J 

The  iiniiting  case  of  concentration  corresponds  to  a  distribution 
function  in  the  form  of  a  puise  at  a  single  radius  irt 

/t(W=o  \  r4  7 

/u.Cr)^o  ,  r.  j 

The  normalizing  condition,  equation  (3.5),  Is,  for 

j  /Ur}  rdc  ^  ^ 


(3.17) 


(3J8) 


(*»*,)  **  ) 


(3.19) 


Figure  30  shows  the  radial  variation  of  the  pressure  coefficient  (r^)  for 
the  first  and  second  tangential  modes  and  Jhe  first  radial  inode.  It  can  be 
seen  that  a  rocket  motor  with  Injection  concentrated  on  a  cirtJe  the  raai'is  of 
which  is  greater  than  approximately  six  lenths  of  the  chanter  radius  is  ii-;'‘e 
unstable  in  -^he  tangential  modes  than  a  corresponding  motpr  with  iJoifot.r.!y 
distributed  combustion,  fhe  wst  unstable  conf igurh+ion  for  these  inodes  is 
that  for  which  the  injection  radius  equals  the  chamber  radius.  For  this  case, 
far  the  first  tangential  mode>  *  ,  ' 

0l7 

On  t.he  other  hand,  the  first  radla}  mode  iV  absolutely  sf^le  when  f-aTa63 
(pressure  node);  and  concentrated  combustion  Is  more  unstable  than  uniform 
only  for  n<04y. 
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An  estimate  of  the  effect  of  the  spreading  of  the  propellant 
streams  on  the-  stability  behavior  can  be  obtained  by  considering  a  rectangular 
injection  distribution, 

ucr)=o  ,  0^ >  < 

fj.  Cr  )=  -^  ,  r:  ^  -  -  r  4  4  A-Tj 

r  2r.A<r^ 

/<  o  ^  r  £  1 

where  4fj,=0  corresDunds  to  the  pulse  distribution.  In  Figure  31,  the  coeffi¬ 
cient  An  >  which  corresponds  to  the  first  tangential  mode,  shown  as  a 
function  of  the  extent  of  the  spreading,  A^  ,  for  three  values  of  the  mean 
radius  c  .  For  this  distribution,  coefficients  greater  than  unity  are 
reduced  and  those  less  than  unity  are  increased  as  A«  is  increased.  For 
moderate  amounts  of  spreading,  the  changes  in  the  coefficients  are  not  large. 
Therefore,  since  the  actual  burning  rate  distribution  Is  not  known  for  a 
rocket  motor  with  radially  concentrated  Injection,  a  good  approximation  to  the 
stability  limits  In  such  a  case  can  be  obtained  by  using  a  pulse  type  distribu¬ 
tion  function. 

A  more  uniform  combustion  distribution  can  be  achieved  by  injecting 
at  more  than  one  radius.  If  the  spreading  of  the  propel j ant  streams  is 
neglected,  the  burning  rate  distribution  functibh  can  be  taken  as 


(3.21) 


Where is  the  3.lrac  delta  function,  i.e.,  a  pulse  at  radius  ifj  ,  and 
-  Is  the  fraction  of  the  total  propellant  flo»<  rate  which  Js  Injected  at 


(3.22) 
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Con 


Then 


sicier  an  injector  tor  which  • 

t^^kCn)  <3.22a) 

{  t*i 

In  order  to  obtain  the  most  uniform  injection  density  distribution,  the 
injection  radii  must  be  such  1hat  the  areas  corresponding  to  the  fj  are  all 
equal.  It  may  be  verified  that 

^  V 

For  such  an  injector,  the  variations  with  the  number  of  injection  radii  of  the 

pressure  coefficients  for  the  first  and  second  tangential  modes  and  the  first 

radical  mode  are  given  in  the  following  table: 

j  =  I  Aji  =  I.I42  A21  =  I.IIO  Ao2  =  0.I33 

2  1. 036  1. 009  O.6OO 

:  I. on  1. 003  0.8II 

Thus,  for  the  tangential  modes,  as  few  as  three  injection  circles  can  produce 

a  pressure  sensitivity  substantial!,  equivalent  to  that  of  an  absolutely 

uniform  injection  distribution.  However,  the  first  radial  mode  requires  a 

larger  number  of  circles  in  order  to  accomplish  a  similor  result. 

With  regard  to  pressure  sensi fivity,  nonuniformity  of  the  injection 

density  In  the  tangential  direction  is  not  significant  for  radial  modes  or 

for  the  spinning  forms  of  tangential  modes.  The  reasons  are  tha+  there  are 

no  diametral  nodal  surfaces  for  radial  modes,  aid  the  tangential  modes, 

the  nodal  surfaces  rotate  with  respect  to  the  injector.  Thus,  In  the  defini- 

tion  of  A^,  ^uation  (3.IOa),  the  produci  I  identically  for  radial 

and  spinning  tangential  modes. 

Assuming  that  the  injection  distribution  Is  uniform  radially, 

=jtC9-)  •  Fo'"  ®  standing  mode,  yO  *  so  that 


<3.23) 


The  normalizing  condition  on  Ji(&)  is 

J  =  aTT 

•'o 

If  all  injection  is  concentrated  along  a  single  radius  6  > 

/i.(S)  ~  Sl&O 


(3.24) 


Ai)U&i)-  (3.25) 

In  a  standing  tangential  mode,  the  pressure  nodal  surfaces  art 
located  at  0*^  and  U  ,  and  the  antinodal  surfaces  are  at  ^  =0 

and  J  .  Then  according  to  equation  (3.24),  if  combustion  is  concentrated  at 
a  node,  the  pressure  coefficient  is  zero,  so  that  the  interaction  Index  required 
for  neutral  oscillaiions  is  infinite.  Thus  the  operation  In  unconditionally 
stable.  On  the  other  hand,  if  combustion  is  concentrated  at  an  ant  I  node,  the 
coefficient  attains  its  maximum  value  of  2.  Then  the  interaction 

index  at  a  stability  limit  is  just  half  of  that  for  a  uniform  combustion  dls- 
tribution.  Positions  of  the  combustion  zone  intermediate  between  the  nodes  and 
antinodes  give  Ai*  values  between  zero  and  two.  In  particular,  injection 
at  results  in  a  pressure  sensitivity  equivalent  to  that  of 

uniform. injection.  These  results  correspond  to  those  obtained  by  the  analysis 
of  longitudinal  modes  for  combustion  concentrated  at  arbitrary  uxial  sfafions  (•). 

When  the  propellants  are  injected  at  more  than  one  discrete  angle, 
the  distribution  function  can  be  written  as 


where  is  the  fraction  of  the  total  flow  rate  Injected  at  B  -G'  and 
»CQi)  is  the  Dirac  function,  as  before.  Then 

i 

For  the  case  in  which  «i. 

I  f 

The  experimental  rocket  motors  used  in  this  investigation  had  injection 


(3.26) 


(3.27) 


(3.27a) 
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c...centrate(]  a1  twelve  angular  positions.  It  can  be  checked  that  =  I 

for  this  case,  so  that  the  experimentai  motors  represented  well  the  pressure 

sensitivity  of  combustion  uniformly  distributed  in  the  tangential  direct 'on. 

Finally,  the  case  of  combustion  concentrated  both  radially  and 

tangentially  wi i I  be  considered.  Let 

Jt  (3.28) 

j-*  ‘  * 

This  burning  rate  distribution  represents  that  corresponding  lo  HI  di screte 
injection  points  ( r?  ,  ©j  )•  The  pressure  coefficient  ror  a  standing  tangential 
mode  1 s  then 


which  is  equivaient  to 

(3.291 

J* 

where  ^rj),  4yJ,  (fij;)  are  defined  in  equations  (3.19)  and  (S.25)  The  remarks 
made  above  concerning  the  iocation,  both  radially  and  tangentially,  of  the 
injection  points  apply  to  this  case  also.  If  the  injection  rate  is  the  same 
at  each  point,  then  y-,  ,  and  equation  (3.29)  simplifies  to 

A'a  <E  ("^i)  j{  X  ^  )J 

which  is  just  the  product  of  the  expressions  obtained  previously,  equations 
(3.22a)  and  (3.27a),  by  considering  separately  the  effects  of  radial  and 
tangential  injection  concentration. 

In  order  to  determine  the  velnr- ry  effects  for  nonUnlform  injection 
distributions,  it  is  convenient  to  divide  equation  (3.14)  through  by  4v4  » 


(3.29a) 


after  inserting  the  definitions 


Thus, 


(',+ i&b  ^  srU' 


The  characteristic  equation  is  now  of  the  same  form  as  equation  (2.78).  Then, 
once  the  pressure  effect  of  the  nonuniformity  has  been  taken  Into  account, 
the  velocity  effecrs  mc'dify  tne  stability  limiis  as  discussed  for  the  uniform 
injection  case.  That  is,  if  the  expansion  coefficient  Cf$,  is  of  the  form 

where  ^  is  a  positive,  real  number,  the  tangential  velocity  will 
cause  a  contraction  or  expansion  in  the  /n,  direction  of  the  stability  limit 
curves  on  the  -fv,  T"  -  plane.  The  effect  of  the  radial  velocity,  as  before, 
will  be  a  shifting  of  the  curves  in  the  TT  direction.  The  magnitudes  of  these 
effects  will  depend  on  the  values  taken  bv  and  .  These 

coefficients  will  now  be  examined  for  several  cases  of  nonuniform  Injection 
distri bution. 

From  equations  (3.10),  the  two  velocity  coefficients  for  a  general 
burning  rate  distribution  ^(r,^  )  are 


Byh  _  X 


(3.3!) 


^  «  -J- 

VfA.. 


jf  f  -j.  J/ ^ 

Cj'^cr.e)  rirde 


(3.32) 


When  the  nonuniformity  Is  in  the  radial  direction  only,  the  coefficients 
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become 


1  Tj(syi,r)  refr 

J  <')  'rdr 


^  yulCr) 


^  _ 


(spinning  mode) 

O 

(standing  rrode) 


(3.33) 


i  I 

V.JAvfc  V  (3.34) 


It  is  clear  that  th-e  convergence  problems  encounterec  previously*  with  the 

tangential  velocity  effect  for  a  standing  mode  are  present  for  any  combustion 

distribution  which  is  uniform  in  the  tangential  direction.  The  effects  of 

nonuniformities  can  be  seen  by  considering  the  case  of  combustion  concentrated 

n 

at  a  single  radius.  For yn.  (v  )  =  ^  (  »^  )  the  coefficients  are 

^Vt.  _L  «v  (3.35) 

»"Al  "  « 

~  '  TT  (spinning  mode)  (3.36) 

The  magnitudes  of  these  velocity  effect  coefficients  relative  to  those  for 
uniform  injection  are  shown  fn  Figures  32  and  33. 

For  the  tangential  'nodes,  both  the  radial  and  tangent i a  1  velocity 
effects  increase  in  relation  tc  the  pr^.. .are  effect  as  the  racius  of  i ejection 
approaches  zero.  However,  since  the  p'-essure  effect  coefficieuT  is  also 
decreasing,  and  at  a  more  rapid  rare,  the  overall  result  Is  that  of  increasing 
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stability  as  tha  radius  of  injection  is  deer  ased.  In  the  case  of  the  first 
radial  mode,  similar  conclusions  can  be  drawn  as  the  injection  radius  appioaches 
the  radius  which  cot  responds  fo  the  pressure  node,  r  =  0.6276. 

Whan  the  nonuniformity  in  the  burning  rate  distribution  is  entirely 
in  the  tanyentiai  direction,  the  velocity  effect  coefficients  differ  from  the 
uniform  injection  coefficients  only  for  the  stanuing  mode,  for  which 

fftofhn,  "  iff  ~P~^Z  7 - 

2Jr 


Ql 


X  J' /A C^)  Cos* d& 


Comparison  of  these  equations  with  equations  (2.56)  and  (2.58b)  yields 


Bw,  , 

f  Brf  ) 

*  ’ 

^  Alfk 

Itr 

Qa _ 

(SiL  ) 

f  CosV(9«rfd 

X  cos  V#  J  9 

Thus,  the  radial  velocity  effect  is  insensitive  to  nonuniformities  in  The 
tangential  direction,  as  would  bo  expected.  In  addition,  for  certain  distri¬ 
butions,  tangential  velocity  effects  will  be  present  for  the  standing  mode, 
whereas  there  are  no  such  '=ffects  for  the  uniform  injection  cace. ,  If  /niS  > 

Is  properly  continuous  and  vanishes  at  9-  =  O  ,  TK  ,  Ihen  the  convergence 
problems  encountered  previously  wil!  no  longer  be  present  because  the  function 
which  must  be  expanded  in  the  eigenfunction  series,  -j; StW  > 
does  satisfy  the  same  boundary  conditions  a.-  rhe  eigenfunctions  themselves. 

For  certain  functions  )  the  coefficient  vanishes.  In  order  to 
determine  the  nature  of  such  distribution  functions,  the  integral  in  the 
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nume  ator  of  equatiori  (3.38)  is  written  ■=>s  the  sum  of  integrals  over 
sub  intervals  of  length  , 

J  +  I  “-tie  +  I  - 

^  ZK-l 


(3.39) 


S  i  nee 

~  Sin  -J  (^~&) 

Cos-Jd  =*  ""  V  -0  )  ^ 

the  second  iiiTegral  on  the  right  side  c*  equation  (3.39)  becomes 

- 

J/*d)«*t*'dctf£V6de  =r  _  )Cosy(^~»)<*(^'^) 

Thus,  if  &)  -  ),  +he  first  two  terms  on  fhe  right  cancel. 

Moreover,  if  ju,(  ^  )  is  periodic,  with  period  ^  ,  al 1  pairs  of  terms  in 
equation  (3.39)  cancel.  Therefore,  in  order  that  the  tangential  velocity 
effect  vanish  for  a  standing  tangential  mode,  the  bur^ning  rate  distribution 
must  satisfy  the  conditions 

//■  -&) 

It  is  clear  that  the  uniform  distribution, I,  is  such  a  function. 

/ 

For  functions^f^  )  which  dc  not  satisfy  the  conditions  (3.'t0), 
the  targer.tia!  velocity  coefficient  for  the  standing  mode  is  real,  rather 
than  imaginary  as  for  the  spinning  mode.  Thus,  for  such  injectors  the  langenlia 
velocity  causes  a  shift  of  the  stability  limits  in  the  ^Tdirectlon  for  the 
standing  mode,  but  in  the /redirection  for  a  spinning  mode.  If  the  combustion 
is  concentrated  along  a  single  radius, ^(  &■  )  =  <^ (  ),  and 


{3.^-0) 


KutA^  \.Tfk/A^J  . 

*'**•?«  rwi 


(sianding  mode) 


(3.4i) 
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Although  tgov^-i**  for  Q.  -  ^  .  ( i+'Zie.  ),  K  =  O,  /  V  the  characteristic 

if  ^  2v 

equation. 


t-l 


^  -  if 

ir (a  /Vs/ 


V  /,  -Cojz"  \ 

j/lr^K  i-e  ^ 


v/ith 

angle. 


■^Colvei  ,  cofV«^  /  shows  that/A,-#flO  for  injection  at,.thir 

If  several  injection  angles  are  uneti.  so  that ^  ^  «r <(>') 


^  \  X7:  s;.v<>-cosV^; 

V«*i 


I.*' 


(3.42) 


which  can  vanish  if  )  satisfies  the  periodicify  and  symmerry  conditions 

of  equation  (3.40).  The  tangentiai  velocity  effect  coefficient  for  injection 
at  a  series  of  discrete  points  (i^j^j  )  can  be  obtained,  in  the  same  manner  as 
used  for  tiie  pressure  coefficient  /y/,  ,  by  combining  equations  (3.36)  and 


(3.42). 


B.  Application  of  the  theory  to  a  Variable  Angle  Sector  Motor 

In  the  experimental  verification  of  the  theory  of  longitudinal  mode 
combustion  instability,  stability  limits  were  determined  on  the  mixture  ratio, 
chamber  length  plane  (16).  From  such  limits,  it  was  possible  to  derive  values 
of  the  interaction  index  and  sensitive  time  lag  as  functions  of  the  mixture 
ratio.  This  experimental  technique  was  made  possible  by  the  fact  that  as 
long  as  the  length  was  kept  greater  than  that  necessary  for  substantially 
complete  combustion,  chat'.^ of  length  had  no  effecl  on  the  combustion  zone, 
i.a.,  on /wand  IT  ,  the  parameters  characterizing  '■he  comb'-stlon  process. 

Since  The  chamber  resonance  properties  were  changed  with  a  change  in  length, 
it  was  pc-ssible  to  survey  the  entire  range  of  stability  behavior,  from  stable 
operation  through  several  harmonics  of  oscillatory  operation, 

A  similar  approach  to  the  study  of  purely  transverse  modes  of 
combustion  instability  is  much  mors  difficult.  The  reason  is  that  the  charac- 
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teristic  length  for  transverse  inodes  i;  the  chamber  radius.  Simply  changing 
the  radius  of  the  combustion  chamber  wl I !  result  also  in  changes  In  fhe 
recirculation  pattern  and  in  the  nozzle  contraction  ratio.  If  the  larter  is 
kept  constant,  the  pressure  level  will  fail,  unless  the  injector  hole  size  is 
increased  or  more  holes  are  added  to  the  injector.  Maintaining  the  length 
to  diameter  ratio  of  the  combustion  chamber  carries  with  it  1he  aanger  of 
entering  regions  of  longitudinal  mode  instability.  Although  some  of  these 
variables,  such  as  the  contraction  ratio  and  chamber  length,  can  be  accounted 
for  theoretically,  the  dependence  of  the  combustion  paramefers  on  the  others 
is  not  known  at  the  present  time. 

For  tangential  modes,  which  are  by  far  the  most  commonly  observed, 
the  possibility  exists  of  measuring  stability  limits  by  using  a  variable  angle 
sector  motor.  Such  a  motor  is  shown  in  Figure  34.  Changes  in  the  angle  can 
be  made  by  means  of  sector  shaped  inserts  which  fit  the  internal  contours 
of  the  combustion  chamber  from  the  injector  face  to  beyond  the  throat  of 
the  nozzle.  It  is  not  necessary  to  extend  the  inserts  very  far  into  the 
supersonic  portion  of  the  nozzle.  For  such  a  motor,  the  characteristic 
length  can  be  taken  to  be  . 

The  use  of' a  sector  motor  is  not  without  difficulties.  In  the 
study  of  longitudinal  modes,  the  chamber  length  can  be  varied  continuous!/, 
i.e.,  in  artitrari I /  small  increments,  so  that  quite  pr-cise  measurements  of 
the  stability  limits  can  be  made,  at  least  in  principle.  However,  it  is  not 
physically  possible  to  construct  an  injector  which  is  absolutely  uniform  in 
the  tangential  direction.  Therefore,  the  sector  angle  can  never  be  continuously 
variable.  Even  accepting  such  a  linit.irion,  which  may  be  of  greater  or  lesser 
importance  in  any  particular  investigation,  further  restrictions  are  imposed 
on  an  injector  to  be  used  for  angular  stability  limit  tests. 
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First.  the  injection  oistribution  mu.,t  be  periodic  in  the  tangential 
direction.  The  smallest  increment  by  which  the  sector  angle  can  be  changeo 
is  then  equal  to  the  period  of  the  injection  pattern.  It  is  reasonabie  to 
expect  the  recirculation  pattern  to  have  the  same  period.  However,  the 
further  requirement  is  imposed  that  the  symmetry  (or  periodicity)' of  the 
recirculation  be  maintained  relative  to  meridianal  planes.  That  is,  there 
must  be  no  swirl,  or  angular  distortion,  of  the  recirculal  ion  patt'^'rp.,  Hence, 
it  is  requirea  that  the  net  momentum  of  the  injec+or  sprays  be  directed 
axially.  If  this  condition  is  not  satisfied,  urburned  propellants  may  impinge 
on  the  faces  of  the  inserts.  The  liquid  layers  thus  formed  Riay  alter  the 
combustion  process  to  a  sufficient  degree  that  the  tests  become  worthless. 

Wall  effects  become  increasingly  significant  as  the  sector  angle  is  decreased. 
Beiow  a  certain  angle,  the  equivalence  between  the  sector  and  the  full  chamber 
motors  may  no  longer  be  valid.  However,  it  must  be  said  that  sector  motors 
of  fairly  small  angle  (27°)  have  been  operated  with  steady  state  performance 
very  close  to  that  of  the  corresponding  full  circular  motor  (9). 

In  a  sector  motor,  the  boundary  conditions  in  the  tangential 
direction  are 


ur'  ie  =0)  -Q 
ur'  (,>=  oc  )=£» 


1 

{  (3.-.2) 


Since  the  proper  solution  of  equation  (2.37c)  for  the  standing  mode  is 

=Cos  f 

the  fundamental  component  of  the  tangential  velocity  perturbation  is 


Sir. 
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Then  the  solid  v;a||  condition  requires  t!  a+ 

.t:«  -sr  0 

which  implies  that 

ILJ  rTy=K^.^.^■■  (3.44) 

where  fti  =  i  corresponds  to  the  first  tangential  mode,  WJ  =  2  to  the  second, 
etc.  The  oscillation  pattern  in  a  sector  motor  of  angle  m.  =1T  Is  thus  the 
same  as  that  of  the  standing  mode  in  a  full  chamber.  The  Bessel  functions  of 
the  first  kind,  ,  are  continuous  in  the  order  V  .  Therefore,  the  zeros  -SWh 
of  the  derivatives  of  the  Bessel  functions  are  continuous  in  V  ,  as  shown  in 
Figure  35.  Thus  for  a  given  tangential  mode,  with  each  sector  angle  oj^is 
associated  a  value  of  the  acoustic  oscillation  frequency.  From  Figure  35  and 
equation  (3.44)  it  can  be  seen  that  Si*,  increases  for  decreasing  angle.  Both 
the  nozzle  admittance  and  the  eigenfunction  expansion  coefficients  are 
functions  of  Sv^, .  The  effects  of  the  sector  angle  or.  By/  ,  andCv/ior 
radially  concentrated  combustion  are  shown  in  Figures  36,  37,  and  38. 

The  center  of  the  combustion  chamber,  'T  =  q  ,  is  a  stagnation 
point  for  the  transverse  velocity  components  for  angles  less  than  IT.  However, 
for  sX  >ir  ,  V  <  1  for  the  first  tangential  mode,  so  that  iT  =  0  is  a  point 
of  infinite  velocity  (for  inviscld  flow).  In  the  flow  ci  a  real  fluid  aroi  nd 
such  a  corner  a  vortex  rattern  is  established.  For  osci .atory  flow,  vortices 
will  be  aiternatoly  created  and  destroyed  cyclicly  on  both  sides  of  the 
corner.  Such  a  phenomenon  constitutes  a  disturbance  to  the  oscillation  pattern; 
it  may  have  an  important,  but  incalculable  effect  in  the  stability  behavior. 

In  considering  the  influence  at  tpe  transverse  velocities  and  of  the 
injection  distribution  nonuniformities  in  a  sector  motor,  it  is  convenient  to 
make  the  change  of  variable 


(3.45) 
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so  that  ranges  between  0  andTT  ior  betwe  n  and  .  Then  the  tangential 
eigenfunctions  are 

=  cos  f’  (3.46) 

With  this  modification,  all  of  the  expressions  derived  previously  for  the 
standing  mode  apply  to  the  sector  motor. 

The  effect  of  the  sector  angle  on  ihe  -lability  limits  for  the  f  ■.  rsT 
tangential  mode  with  no  velocity  effects  is  demonstrated  by  Figures  39,  40, 
and  41.  These  curves  were  calculated  for  the  axially  concentrated  steady 
state  combustion  di s-' ; tution  which  was  discussed  in  Chapter  JI..  Figure  39 
shows  the  limits  for  sector  angles  of  60®,  90°,  and  180°  fo'*  a  uniform  trans¬ 
verse  combustion  distribution.  The  slight  differences  in  the  minimum  intei — 
action  index  values  are  caused  by  the  changes  in  the  acoustic  frequency 
through  the  nozzle  admittance  coefficients.  The  solution  for  <*iris  affected 
to  only  a  slight  extent  by  the  value  of  Syj,  .  There  fore, T  is  approximately 
proportional  to  .  As  the  sector  angle  is  decreased,  the  region  of 
instability  on  the /rv, Z  -  piane  for  any  particular  mode  moves  toward  smaller 
T  .  For  a  given  value  of ,  the  range  of  unstable  time  fag  values  also 
decreases  with  decreasing  angle. 

The  effect  of  concentrated  comtrasti en  on  the  stabi I :ty  I imits  is 
shown  in  Figure  40  for  =  8/9  and  in  Figure  4i  for  C  *  5'9.  These 

injection  radii  correspono  -o  those-  used  in  the  experimental  *  •  nch  .di arneter 
rocket  motor.  For  injection  on  the  larger  circle,  the  unstable  zones  extend 
to  lower  'K' as  the  sector  angle  decreases.  On  The  other  hand,  for  the  smaller 
injection  radius,  the  sector  motor  becomes  more  stable  as  the  angle  is 
decreased.  Although  there  is  no  radius  for  cii  the  rnininum  interaction 
index  is  independent  of  the  sector  angle,  from  Figure  35  it  can  be  seen  that 
the  smallest  variation  of  with  o<  occurs  for  y?  ■=  0.76. 
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From  Figures  39  ,  40,  and  41,  it  is  apparent  that  any  point  (./%.  ,'d  ) 
which  is  within  an  unstable  region  lies  on  the  intersection  of  the  stability 
limit  curves  for  two  sector  angles.  Therefore,  there  is  a  range  of  angles 
for  which  the  rocket  motor  represented  by  fhaf /n.  pair  is  unstable  (in 

the  mode  considered).  Thus,  the  possibility  exists  tor  determining  values  of 
the  interaction  index  (for  pressure)  and  the  sensitive  time  lag  by  measuring 
angular  stability  limits.  Since  the  velocity  indices  and  are  not 
presently  calculable,  these  stability  limit  tests  must  be  carried  out  using 
velocity  insensitive  injectors,  such  as  a  tangentially  oriented  spud  in¬ 
jector  (see  Figure  7)  with  a  sufficiently  large  number  of  spuds.  If  the 
pressure  index  and  time  lag  are  known,  then  it  appears  possible  to  determine 
the  radial  velocity  index  j2r  by  noting  the  difference  between  the  pred idled 
and  actual  angular  stability  limits  using  the  radial  spud  orientation.  Then, 

If  and^^  '  where  is  the  spud  orientation  angle, 

as  postulated  in  Chapter  II  ,  the  tangential  velocity  index,  which  is  important 
In  the  instability  of  the  full  circular  motor,  is  determined  at  the  same  time. 

C.  Generalized  Time  tag  Formulation 

In  the  previous  chapter,  an  expression  for  the  burning  rate 
perturbation  was  derived,  following  (!),  on  the  basis  of  the  folinwlng  time 
lag  schematizar ion.  The  actual,  gradual  evolution  of  combustion  products  f-CMB 
the  liquid  propellants  was  approximated  by  an  instantaneous  conversion, 
af+er  a  certain  time  interval  from  injection,  called  the  total  time  lag  (  Ty  ). 
It  was  assumed  that  the  combustion  process  was  not  sensitive  to  fluctuating 
conditions  in  the  combustion  chamber  during  the  first  part  of  the  total  time 
lag,  cailad  the  insensitive  time  lag  (T*  >.  During  the  latter  part,  the 
sensitive  time  lag  (X  ),  the  rates  of  the  processes  occur! ng  during  combustion 
were  assumed  to  depend  on  the  pressure  and  on  the  transverse  components  of 
the  gas  velocity.  The  degree  of  sensii-ivity  of  the  combustion  prucess  to  the 
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chamber  conditions  was  represented  by  the  pres  ure  interaction  index  M/ 
and  the  veiocity  indices  Ar  and  Ji^  .  Thus,  considering  combustion  to 
occur  at  time  t  ,  the  interaction  indices  were  assumed  to  be  zero  from 
t-'Cj  to  t~Z  ,  and  to  have  cons'^ant  values  from  i-Z"  to  ^  .  Stability 
criteria  were  deriveo  on  the  basis  of  this  time  lag  scheme.  The  effect  of 
alternative  time  lag  formulations  on  the  characte;  istic  equation  for  neutra. 
oscillations  will  be  considered  in  this  section. 

Ths  step  function  combustion  model  defining  riie  total  time  lag  will 
be  retained  in  the  foi Sowing  analysis.  However,  the  simple  division  of 
into  a  sensitive  and  insensitive  parr  wiii  be  replaced  by  th«  more  general 
scheme  shown  in  Figure  42  (b).  According  to  this  model,  the  combustion 
process  is  sensitive  to  velocity  effects  only  during  the  time  interval  TJ  to 
7^  and  to  pressure  effects  only  between  “C^and  TJ*.  Therefore  'tymXr'X,  is 
the  veiocity  sensitive  time  lag  and  pressure  sensitive  time 

lag.  A  straightforward  calculation,  following  the  procedure  used  in  Chapter 
II  ,  yields  an  expression  for  the  burning  rate  perturbation. 

The  overall  rate  of  the  processes  involved  in  the  combustion  of  - 
the  liquid  bi propel  1  ants  is  specified  as  follows: 

4'f  y  -fc  s  t  -r^.  f-r, 

^  t-Tr+Zt  i  €  -t 

t-rT+ri.  <  l'<  t-tf 

-t-Cjr  +rj,  £  $  -fr-Tr+xi^ 

where  the  index >^anri  the  perturbation  rrfer  either  to  the  radial  or  to 
the  tangential  velocity,  or  to  both.  The  total  time  lag  is  defined  by  the 


^  <3.47) 
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following  equation 

t 


i-z,  i-r, 

which  expresses  the  condition  that  fhe  t'^ansformation  from  reactants  to 
products  takes  plac-e  when  the  processes  occurring  during  the  time  lag  have 
accumulated  to  a  certain  level. 

Introducing  the  expressions  (3.47)  into  equation  (2.17)  gives 


r  "t 

■h  -f- 


Upon  rearrangement  squat  ion  (3.48)  becomes 

)Ht'  -  -  ffai' 

d-zr  4,,  4^ 

i-£r  ■6-‘Cr+‘4 

Assuming  the  spatial  nonuniformity  of  ^  to  be  negligible,  the  inli-grc!  on 
the  left  side  can  be  written  as 


so  that 


Xt-Xj  -  dt 

<-v*i 

neglecting  higher  order  terms. 


r 

-/V. 


i-Tr+Xi 


(2.17) 


(3.48) 


(3.49) 


(3.50) 
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Tne  burning  rat:  pferturbation  was  found  in  Chapter  II  to  be  related 
to  the  time  rate  of  change  of  the  time  lag  by  the  expression 

<2. 16) 

Differentiation  of  equation  ■;3.50)  wi1h  respect  to  tinie  and  making  use  of 
the  definitions  of  ana  'C^  gives 


n  '  sCi-Vri-x,)  , 

e  ~  *  }  *  *’'3*^®'"  ■hsrms 

Then,  since- the  burning  rate  perturbation  is  of  the  form 


,  I  4u 


c5-52>y 


The  proportionality  factors  are,  from  equation  (3.5i), 

-  ,  _) 

These  expressions  are  to  be  compared  with  those  of  equation  (2.32). 

From  equation  (3.53)  it  can  be  seen  that  the  more  general  for.Tiula- 
tion  of  the  combustion  process  requires  the  use  of  eighv  parameters;  the  total 
time  lag  ^  ;  the  pressure  interaction  index />r>  and  time  lags  ; 

the  volocilf  inaices  j2- ,  -k^  and  time  lags  ,  ■?,  .  If  a  relation  exists 

between  j?,- and  ,  ac  in  the  case  of  the  unlike  impinging  doublet  spud 
injector  considered  in  Chapter  II  ,  then  the  number  of  parameters  required 
i s  reduced  to  seven . 
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is  a  conrpiex  relation  between  the  combustion  parameters  which  must  exist  If 
neutral  oscillations  of  frequency  are  occur.  A  solution  can  be  obtained, 
as  discussed  previously,  only  if  all  but  two  of  the  parameters  are  known.  At 
the  present  time,  it  is  not  possible  to  calculate  any  of  the  sensitive  or 
insensitWe  time  lags  or  interaction  indices.  However,  the  pressure  sensitivity 
theory  hac  bee"  applied  successfully  to  the  analysis  of  longitudinal  modes 
of  instability..  According  to  this  formulation, 

so  that 

(9  =  ^  (i-e  )  (!1.55) 

on  the  basis  of  this  longitudinal  mode  theory,  a  n.e+hod  has  been  developed  for 
measuring  and  experimentally.. 

The  time  lag  model  used  i  previous  sections  of  this  dissertation 
constitutes  the  simplest  extensicr.  of  the  pressure  sensitivity  theory  to 
include  velocity  effects.  Thus,  the  velocity  time  lags  were  assumed  to  be 
identical  to  the  pressure  time  lag.  In  the  preliminary  experimental  stuoi-?s, 
it  was  observed  that  the  velocity  effects  were  more  important  in  the  early 
combustion  region  than  in  the  zone  of  greatest  combustion  Intensity.  There¬ 
fore,  a  better  approximation  to  the  actual  combustion  process  than  the 
simple,  commonH-ime-iag  model  may  be  obtained  by  setting  T,  »  O  ,  so  thi-t 

/  (3.56) 

while  retaining  equation  (3.55)  to  describe  the  pressure  effect. 

It  is  seen  that  the  total  time  lag  is  involved  in  the  separate 
time  lag  formulation.  The  total  time  lag  for  a  propellant  element  is  definec 


0 


(3.57) 
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where  3  is  the  station  at  which  the  elemen"!  undergoes  conversion  to 


combustion  products.  If  combustion  is  concen'T rated  at  a  single  axial  station 
'Em  and  all  elements  have  the  same  velocity,  equation  (3.56)  becomes 

2m 


r  dl  - 
~  J 


where  ULm  is  a  mean  liquid  velocity.  For  this  case  the  functions  fi  and 
^  are  independent  of  i  .  If  the  combustion  is  distributed  axially,  then 
C?=(X(t)and  so  that  the  form  of  the  characteristic  equation  is 

somewhat  modifi-d.  Equation  (2.53b),  for  d  =1^,  q  =  h,  becomes 

^  /  Prui,  \  q;  -St,  ft)  Jr- 


where  ^  ~  1 )  •  Then 


d  f  Pevt» 
(iiV  P«o 


\  I 

;  =  6<  J 


(3.53) 


(3.59a) 


i 


«  f  fCgA 


(3.59b) 


Similar  results  are  obtained  for  .  If  a  function  h,  /  is  oetined  as 

T,  =  r *  -  i  S't.t  8f<>if  ' 


(3.60 


then  the  characteristic  equation  for  the  case  ot  dist'ihuted  combustion  t-xes 


the  form 


(3.61) 
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anci  solutions  are  obtai". 3d  in  the  same  ma-ner  as  for  the  concentrated  case. 

An  imoortant  application  of  either  equation  (3.54)  or  (3.61)  is  the 
determination  of  the  velocity  effect  parameters  once  the  pressure  interaction 
index  AX,  and  time  lag  are  known.  It  is  clear  that  a  value  of  one  of  the 

velocity  indices  ,  Jig  ,  or  a  relation  between  the  two,  must  be  known 
order  to  fine  solutions  of  the  charOwf— equation.  Since  both  rhe 
radiai  and  the  tangential  velocity  effects  are  of  the  same  for'm,  it  is 
necessary  to  discuss  the  solution  for  orly  one  of  them.  Thus,  let  Jl^=  0  In 
the  foi  lowing.  Fc  simplicity,  the  solution  of  equation  (3.54)  ..ill  be 
considered. 

For  this  case,  the  characteristic  equation  can  be  written  in  the 


form 


(3.62) 


where 

Since  equation  (3.62)  is  complex,  it  Is  equivalent  to  the  two  real  equations 
Jif  —  Jlf  4o  s  —  9»- 


<3.63) 


jZr 


=  3i 


The  solution  of  this  sirnultaneous  pair  is 

■2.S*- 


(3.64) 


The  most  important  difference  between  the  "common"  and  "separate" 
+ime  lag  formulations  is  the  presence  of  the  total  time  lag  as  a  parameter  In 
the  latter  model.  Since  ^  is  0(  and  is  0  (1),  the  total  time  lag 
is  0(  ).  Thus,  small  variations  in  the  frequency  tO  .  which  in  the 

method  of  solution  is  regarded  as  the  independent  variable,  causa  large 
changes  in  the  function  9  ('o)  Therefore,  it  is  important  that  the  combustion 
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distribution  tZ  Or)  and  the  liquid  velocity  ?  )  oe  accurately  known  if 
a  valid  solution  for  The  combusfion  parameters  is  to  be  obtained  from  +hs 
separate  time  lag  < ormu i at i on .  It  can  be  concluded,  then,  that  further  study 
of  the  details  of  the  combustion  process  is  needed  in  order  to  refine  the 
time  lag  nodai  properly. 


1 

i 
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CHAPTER  IV 

Exper imenta;  Results  and  Conclusions 

A.  General  Considerations 

The  stability  limit  relations  for  transverse  modes  ui  instability  in 
the  two  p''eceding  chapters  were  based  on  th"  comcept  or  a  sensitive  combustion 
time  lag.  The  validity  of  this  model  as  apolied  to  longitudinal  modes  has 
been  demonstrated  experimentally.  Howevetj  the  extension  io  transverse  modes 
involved  more  than  just  a  geometric  change  in  the  number  of  dimensions.  It 
was  necessary  also  to  include  in  the  combustion  model  a  sensitivity  to  velocity 
as  '.veil  as  pressure  oscillations,  to  treat  nonuniform  injection  distributions, 
and  to  develop  a  ^■etholj  o-  iiieasuring  stability  limits  at  which  the  linearized 
theory  covid  be  applied.  None  of  the  effects  considered  had  been  studied 
systematically  by  any  previous  worker.  Therefore,  the  final  phase  of  the 
present  investigation  consisted  of  a  series  of  experiments  designed  to  test 
tne  ■hneoreticc I  predictions. 

The  experimental  apparatus  used  in  these  tes'. s  was  the  same  as  that 
used  in  the  exploratory  tests  described  in  Chapter  I.  All  tesTs  were  made  at 
a  mean  chamber  pressure  of  150  psia  and  at  a  thrust  level  of  '000  nuunos 
based  on  operation  of  the  full  circular  chamhe'-  The  injector  spuds  were 
designed  to  give  equal  propellant  velocities  and  axiaiiy  directed  spray  fans 
at  a  mixture  ratio  of  1.4.  Most  of  ii>e  effort  was  devoted  to  the  measurement 
of  angular  stability  limits,  using  the  variable  angle  sector  motor.  A  limited 
number  of  full  motor  tests  were  made  it  order  to  observe  the  wave  patterns  of 
the  spi.ining  forms  of  tangential  modes.  In  some  of  tnese  tests,  destruct  ole 
barriers  were  used  to  determine  the  role  of  starting  transients  in  the 
stabli'ty  probiem.  The  number  of  full  motor  tests  was  necessaniiy  restricted 
by  the  extremely  destructive  nature  of  the  spinning  tangential  modes.  Tests 
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were  also  conducTeo  /vi1h  a  variable  I  ngth  motor  in  order  to  use  the  well- 
developed  longitudinal  mode  stability  limit  method  to  measure  the  nressure 
interaction  index  and  sensitive  time  iag.  Correspondence  between  the 
ioiigitudinal  and  transverse  motors  was  achieved  by  making  the  cross  section 
of  1he  variaole  iength  motor  in  the  snape  of  a  30°  sector,  as  shown  m 
Figure  43.  In  this  way  the  nature  oi  ih*  recirculation  pairarn  mintained 
the  same  in  fhe  longitudinal  tests  as  in  the  transverse,  provided  that  the 
injector  satisfies  the  conditions  required  by  the  sector  motor  concept,  which 
were  discL.s::>a  in  Chapter  III. 

In  the  experimentai  program  the  propel iant  r.ombina+lon  and  injector 
type  were  held  constant.  The  comoination  of  liquid  oxygen  and  ethyl  alcohoi 
was  chosen  on  tne  basis  of  the  large  amount  of  experience  gained  in  the  use 
of  this  combination  in  tne  previous  investigation  of  longitudinal  modes  of 
i nstao i i ity..  The  unlike  impinging  doublet  injector  was  selected  because  of 
the  nature  of  its  symmetry.  The  spray  formed  by  this  type  of  injector  can  be 
expected  to  be  sensitive  to  mixing  effects  of  the  velocity  component  parallel 
to  the  line  of  centers  of  the  unlike  propellant  orifices  but  not  of  the 
component  norma!  to  the  line  of  centers.  The  use  of  the  spud  injector,  with 
the  spud  orientation  variable,  facilitated  the  separate  study  of  the  eiects 
of  the  radial  and  tangential  velocity  components.  However,  rne  spi  ay  angle 
with  respect  to  t^c  combustion  chcmbe'"  axis  i<5  a  runction  of  the  mixture  ratio. 
Only  at  the  design  mixture  ratio  is  the  net  spray  momentum  of  the  tangential 
spud  orientation  directed  axially.  Thus,  a  swirl  component  is  introduced  into 
the  recirculation  pattern  for  off-design  operation,  and  serious  trouble  may 
be  encountered  from  wal !  effects,  esnsc.'ally  for  small  sector  angles.  There** 
fore,  fills  type  of  injector  is  not  suitable  for  sector  tests  over  a  wide 
range  of  rpijfture  ratios,  so  that  the  results  of  the  stability  limit  tests 
can  be  regarded  as  valid  only  near  The  design  point. 
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Tie  determination  ot  the  variaL>e  angle  stability  limits  was 
effected  by  means  of  a  large  number  ot  rockcf  test  runs.-  Indivioual  tests 
were  quite  simple  The  motor  was  operated  at  preselected  values  of  the 
mixture  ratio  and  chatriper  pressure,  wnich  were  monitored  ay  steady  state 
gages  and  f lowmete.  ^  Ttie  primary  measurement  was  the  oscillatory  component 
OT  the  chambe*'  pressure.  BoTn  amplitude  and  ri-equency  data  were'  recorded. 

The  determination  of  the  node  of  instability  was  made  on  the  basis  of  the 
frequency  of  the  oscillations  It  has  been  observed  by  many  investigators 
that  tne  freouer.cy  of  purely  transverse  modes  is  always  within  Ten  percent 
of  the  frequency  of  th^^acoustic  modes  in  a  closed  cylifider.  This  result 
is  obtained  even  in  the  C'a.se  of  very  large  amplitude  oscillations, 

Tne  exact '„v^:^t.de  level  of  the  oscillation  in  an  unstable  motor 
is  of  liftle  significance  to  a  linear  theory.  However,  the  variation  ot  the 
amplitude  as  a  stability  limit  is  approached  can  be  an  indication  of  the 
location  of  the  limit  This  fact  is  of  particular  importance  in  view  of  the 
rather  large  sector  angle  increment  required  by  the  experimental  apparatus. 
Figure  4<i  snows  fhe  tvpica!  variation  of  the  oscillation  amplitude  'nith 
sector  angle  for  a  constant  value  of  mixture  ratio  The  transition  from 
tne  first  to  the  second  tangential  modes  was  generally  quite  sharp.  In  1he 
case  shown  in  Figure  44,  a  very  small  mixture  ratio  shift  at  an  angle  of  I80^ 
C3"sed  a  change  ‘n  the  ie  with  little  variation  of  ampliiude.  However,  the 
limits  between  stable  operation  ana  the  first  tangential  mode  of  instability 
were  less  clearly  defined.  Many  cases  of  organized  oscillations  alternating 
irregularly  with  random  combustion  noise  were  observed.  At  the  small  angles 
•associated  with  these  limits,  the  combustion  nc ise  level  was  ouite  high,  often 
exceeding  lO  psi  (rms). 

In  most  of  the  stability  limit  tests  a  single  pressure  transducer 
was  used  Tne  amplitude  of  a  standing  mode  oscillation  varies  approximately 
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as  Cos  -—4  ^  whe!  Q-=  o  and  correspond  to  the  flat  walls  formed 

oc  ^ 

by  the  sector  inserts.  For  large  scale  oscillations,  the  wave  shape  departs 
from  the  sinusoidal,  so  that  tns  amplitude  does  not  vanish  at  the  nodes 
as  shown  iii  Figure  45.  In  order  to  obtain  the  ma.ximum  amplitude,  the  pressure 
transducer  was  placed  as  near  as  practicable  to  one  of  the  flat  wails.  For 
fui  I  motor  tests,  and  in  some  secTOi'  tesl.,  two  or  three  prfe.,sur-  pi'.'.r.ups 
were  inserted  at  90°  intervals  around  the  chamber  circumference.  In  this 
manner  it  was  possible  to  detect  whether  a  particular  oscillation  was  of  the 
spinning  or  sti..iding  type.  By  means  ot  observations  of  the  relative  phases  of 
the  pressure  traces  the  direction  of  travel  of  a  spinning  wave,  which  was  of 
great  significance  in  the  velocity  effect  studies,  could  be  determined. 

The  experiments  were  organized  into  three  series.  First,  sector 
motor  stability  limit  tests  using  the  tangential  orientation  of  the  injector 
spuds  were  made  in  order  to  verify  theoretical  predictions  concerning  the  ef¬ 
fects  of  chamber  diameter  and  injection  diameter  in  the  absence  of  velocity 
effects.  The  bases  for  this  series  were  ihe  theoretical  demonstration  of  the 
absence  of  tangential  velocity  cffsc+s  in  standing  modes  and  the  nature  of  the 
symmetry  of  the  unlike  doublet  type  of  injector.  The  second  series  of 
experiments  utilized  fu'l  circular  chambers  in  order  to  study  the  tar.jenrial 
velocity  effects.  In  particular,  these  tests  sought  +o  determine  ihe  positive 
velocity  direction  lith  respect  to  the  injector  dsubie*  Using  the  Information 
obtained  from  the  second  series,  the  sign  of  the  radial  velocity  index  could 
be  determined.  The  third  series  of  experiments  was  then  devoted  to  measuring 
the  angular  stability  limits  using  the  radial  spud  orientation.  These  tests 
were  designed  to  ascertain  the  validity  of  the  theoretical  prediction  that 
the  radial  velocity  causes  a  shift  of  the  angular  stability  limits  in  a 
direction  which  depends  on  the  sign  of  the  velocity  index.  The  results  of 
the  experimental  program  will  be  discussed  in  connection  with  the  corresponding 
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theoretical  calculations  in  the  following  sections. 

B.  Sector  Motor  Tests  with  Tangential  Spud  Orientation 

The  angular  stability  limits  shown  in  Figure  46  are  typical  of  tne 
results  obtained  from  the  sector  motor  tests.  A  point  representing  each 
test  has  been  plot+sd  on  the  mixture  ratio,  sector  angle  plane  with  a  symbol 
rtiiich  describes  The  nature  of  the  operation,  o.  g.  ,  stable,  firsi  tangwu  ial 
mode,  etc.  At  small  angles,  the  rocket  motors  operate  stably.  As  the  angle 
is  increased,  a  regio"  of  first  tangential  mode  siabllity  is  found.  The 
second  tangential  node  is  observed  at  angles  greater  than  180°.  No  o+her 
modes  have  been  observed  in  these  tests.  Although  the  transition  from  the 
first  to  the  second  mode  is  sharp  in  general,  there  is  some  difficulty  in 
obtaining  reproducibi I i+y.  Possibly  this  effect  is  a  result  of  overlapping  of 
the  unstable  regions  for  the  two  modes  (see  Figuie  21).  At  the  lower  limit 
of  the  first  tangential  mode  the  30°  increment  in  sector  angle  restricts  the 
accuracy  obtainable  in  the  definition  of  the  stability  limit. 

The  variation  of  the  propellant  spray  angle  with  mixture  ratio  results 
in  the  possibility  of  effects  caused  by  distortion  of  the  recirculation  pattern 
and  liquid  impingement  on  the  walls.  Therefore  one  is  advised  to  draw  con¬ 
clusions  only  from  the  tests  made  in  the  vicinity  of  the  design  mixture  ratio. 
However,  it  has  been  found  that  by  surveying  a  rather  wide  mixture  ratio 
range,  it  is  possible  ■’•c  compensate  somewhat  for  the  coarsr-'ess  of  the  sector 
angle  variation.  As  an  example,  in  Figure  46  it  is  apparent  that  the  limit 
between  the  first  and  second  modes  at  a  mixture  ratio  of  i.4  lies  close  to  180°. 
If  testing  had  been  confined  only  to  that  one  mixture  ratio,  the  limit  could 
only  have  been  described  as  being  somewk-’ra  between  180°  and  210°.  For  the 
most  part,  however,  the  variation  of  the  stability  limits  with  mixture  ratio 
is  very  small,  unlike  the  results  obtained  for  longitudinal  mode  stability 
limits  (16). 
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The  tests  -hown  in  Figure  4f  correspond  to  a  rocket  motor  with  a 
9  inch  combustion  chamber  diameter,  for  which  the  chamber  exit  Mach  number 
was  approximately  0.05.  The  propellants  were  injected  near  the  outer 
periphery,  on  an  8  inch  diameter  circle.  The  theoretical  angular  stability 
limits  corresponoing  to  this  configuration  are  shown  5n  Figure  47.  Because 
of  the  theoretical  result  that  tai.i,'.rf''’.!  velocity  effec+-;  are  absent  in 
standing  modes,  and  the  fact  that  radial  velocity  effects  on  the  spray  prod¬ 
uced  by  a  tangentially  oriented  doublet  can  be  expected  to  be  quite  smaii, 
the  iimif  c'.T-es  of  Figure  47  were  calculated  assuming  only  assure  sensitivity 
of  the  combustion  process. 

Considering  the  stability  limits  at  the  design  mixture  ratio  (1.4), 
it  is  seen  that  the  )  point  representing  this  particular  motor  must 

fall  within  the  small  oval  shaded  area.  Thus  the  in+eraction  index  must  lie 
between  0.60  and  0.75  and  the  time  lag  between  0.11  ana  0.14  millisecond.  In 
previous  stability  limit  tests  using  an  unlike  impinging  doublet  injector  with 
a  design  mixture  ratio  of  1.4,  tne  values  obtained  for  the  combustion  parameters 
were:  /W/  =  1.7,  X*"  =0,  5  millisecond.  The  most  sti  iking  difference 

between  these  results  is  the  greatly  reduced  interaction  inaex  as  deiermined 
from  the  transverse  mode  stability  limit  tests.  However,  the  longitudinal 
mode  tests  were  made  with  a  3-inch  diameter  chamber,  at  a  tb’-wst  le/sl  of  25U 
pounds  and  a  chamfer  pressure  of  300  psia.  rhersfo-e,  the  comparison  is  not 
real ly  val id. 

In  order  to  determine  whether  the  combustion  parameters  determined 
from  the  tangential  orientation  sector  tests  agree  with  those  which  can  le 
obtained  from  longituainal  tests,  was  necessary  to  construct  and  operate 
a  variable  length,  30°  sector  motor  tests.  The  use  of  the  sector- shaped  cross 
section  assured  the  matching  of  the  recirculation  pattern  in  the  variable 
angle  and  variable  length  motors.  However,  comparison  of  the  values  of  fiie 


combustion  porameters  wiih  the  fheoretica!  longituciinai  mode  stability  pattern, 
as  shown  ii'  Figure  48,  indicated  that  no  instability  could  be  expected  with 
the  variaDle  lengtn  motor  because  of  the  low  interaction  index  and  sens!'!ive 
tin.e  lag  values.  The  longitudinal  mode  tests  are  presented  in  Figure  49. 

It  can  bo  seen  tha+  the  theoretical  prediction  of  stability  was  indeed  verified, 
l-luwever,  this  result  cannot  be  considered  us  u  positive  demonstration  of  rhe 
identity  of  the  pressure  effect  combustion  parameters  in  both  longitudinal  and 
transverse  modes. 

The  sr.~e  injector  spuds  were  also  used  in  a  6  inch  cii omot''’-  combustion 
chamber  at  the  same  levels  of  pressure  and  thrust,  so  thj't  the  exit  Mact>  number 
was  increased  to  0.10.  Propellant  injection  was  on  a  5  inci'  diameter  circle, 
slightly  smaller  in  relation  to  the  chamber  diameter  than  that  of  the  preceding 
case.  Figure  50  shows  the  theoretical  stabilitv  limit  curves  corresponding 
to  this  motor  configuration.  It  should  be  noted  that,  if  the  dimensicnat  time 
lags  are  the  same  for  the  6  and  9  inch  motors,  the  noi.dimensiona!  time  lags 
are  in  the  inverse  ratio  of  the  diameters,  i.e.,  =  ^  .  However,  the 

experimental  stability  limits  tests  resuitec  in  .a  quite  different  range  of 
possible  values.  From  Figure  51  it  can  be  seen  thaf  the  lower  limit 

remained  between  60°  and  90°  and  the  upper  limit  was  located  between  150°  ar.a 
180°,  Therefore,  the  v-aiues  of  the  combustion  parameters  lor  the  6  inch  motcr 
differ  from  tnose  obta'-ed  with  the  9  inch  motor.  The  in+c  action  index  and 
sensitive  lime  iag  values  corresponding  to  t.he  s'rabiliiy  limits  of  the  smaller 
motor  are;  .65  ,90;  .73<'C**'<1  .93  (milliseconds).  The  primary 

difference  between  the  two  cases  is  see.n  to  be  in  the  time  lag  value. 

Very  little  is  known  concerning  ..le  dependence  of  the  parametei's  /Kv 
and  t:  on  the  physical  conditions  in  the  combustion  chamber.  Since  the 
chamber  pressure,  propel Icint  liow  raie^,  and  injector  spud  geometry  were  the 
same  for  the  9  inch  and  6  inch  motor.s,  the  cone  i  is  ion  can  be  drawn  that  fhe 
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decrease-  in  the  sesi^itive  tiaie  lag  we  j  brought  about  by  the  lateral  contraction 
of  the  recirculation  zone.  It  should  be  pointed  out  that  the  theoretical 
curves  were  calculated  on  the  basis  of  combustion  infinitely  concentrated  in 
the  radial  direction.  However,  it  was  shown  in  Chapter  III  tnat  the  amount 
of  spreading  expected  for  toese  injectors  would  contribute  an  error  of  less 
than  t Ko  percent,  which  is  much  less  thc..  rhe  uncertainty  in  the  e.xper imentr.i 
stability  limits. 

The  third  series  of  tangential  spud  orientation  angular  stability 
limit  tests  ass  designed  to  investigate  the  effect  of  the  injection  diameter 
on  stability.  For  these  experiments  the  9  inch  cor>'*'*.istion  chamber  was  utilized, 
with  the  S2»me  injector  spuds  as  used  previously  arranged  on  a  5  inch  diameter 
circle.  According  to  +he  distribution  function  treatment  of  nonuniform 
Injection,  the  unstable  regions  on  the  /JU  ,  I"  plane  shift  toward  higher 
interaction  index  values  as  the  injection  diameter  is  decreased  relative  to 
the  combustion  chamber  diameter.  As  shown  in  Figure  52,  an  interaction  index 
greater  than  unity  is  required  for  Instability  when  injection  is  concentrated 
at  5/9  of  the  chamber  diameter.  Assuming  that  the  combustion  parameters  are 
not  greatly  affected  by  the  change  in  the  recirculation  patlei  n,  the  prodiction 
is  made  that  the  rocket  motor  win  operate  stably.  The  experimental  i  asults 
shown  in  Figure  53  are  seen  to  be  mostly  stable.  The  combustion  noise  'evel 
on  these  tests  wa-  higher  than  that  observed  in  the  previous  series, 
averaging  15  psi  (rms).  For  some  tests  at  angles  greater  than  180°,  oscil¬ 
lations  at  frequencies  corresponding  to  one  or  another  of  the  transverse 
modes  dominated  the  pressure  fluctuations  intermittently.  Although  these  tests 
are  whown  in  Figure  53  as  " intermlt'‘‘r,.T  'pstabi  1  ity"  they  are  not  an  indicf.+  ion 
of  instability  but  rather  due  to  the  p-'eferential  amplification  associated 
with  the  resonance  properties  of  the  combustion  chamber. 

Although  hhis  test  series  established  qualitatively  the  validity 


of  the  theoretical  appi o^ch  +o  nonuniform  'njection  distributions,  it  was  not 
The  li'jant i tat i ve  proof  one  would  desire.  In  order  to  accomplish  a  more 
positive  ve*- i f i cat i OP  an  injec+ion  diameter  intermediate  between  5  and  8  Inches 
shouid  be  used,  in  order  that  there  be  a  reoion  o!  instability,  from  the 
limits  of  which  values  of  the  combustion  parameters  coulo  be  obtained. 

However,  the  results  of  tiie  6  inch  diomc't-  .T.^'+jr  tests  have  denio.Tstrattd 
that  the  time  lag  and  interaction  index  are  not  insensitive  to  changes  in  the 
injection  geometry.  Therefore,  quantitative  verification  of  the  theory  for 
the  injection  oistrioution  effect  may  not  oe  possible. 

C.  Full  Motor  Tests  with  Tangential  Spud  Orientation 

From  the  linearized  velocity  effect  theory  the  result  is  obtained 
that  the  tangential  velocity  effects  are  absent  in  standing  modes  and  present 
in  spinning  modes.  The  tangential  velocity  perturbation  acts  to  shift  the 
stability  limits  on  the  At  ,  'C  -plane  to  smaller  values  of  the  interaction 
index.  The  effect  of  the  radial  velocity  perturbation  is  to  shift  the 
stability  limits  to  larger  or  smaMer  time  lag  values,  depending  on  the  sign 
of  the  velocity  index,  A  positive  velocity  index  is  associated  with  the 
tangential  velocity  affect.  However,  the  positive  tangential  direction  is 
not  defined  in  terms  of  the  combustion  chamber  georrnetry,  as  in  the  case  cf 
the  radial  velocity.  The  positive  tangential  direction  is  defined  as  me 
direction  in  which  a  spinning  wave  travels. 

The  sign  of  the  velocity  index  depends  on  the  characteristics  of 
the  spray  formed  by  the  injector.  According  to  the  theoretical  model,  a 
posixive  velocity  index  implies  that  the  burning  rate  is  increased  when  the 
velocity  perturbation  is  in  the  positive  direc+ion.  The  symmetry  of  the 
unlike  impinging  doublet  injector  is  such  that  the  radial  and  tangential 
velocity  effects  can  be  separated,  as  discussed  previously.  For  the  radial 
spud  orientation,  the  velocity  effect,  due  solely  to  the  radial  velocity 
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componeiit,  cannot  be  calculated  unles  the  sign  of  the  velocity  index  is  known 
for  a  particular  injector  arrangement.  At  the  present  time  knowledge  of  the 
details  of  the  combustion  process  is  insufficient  to  allow  the  predirlion  of 
tne  velocity  index  sign.  Therefore,  one  must  resort  to  experimental  means, 
based  on  spinning  mode  tests  with  the  tangential  spud  orientation.  By 
observing  the  direction  of  wave  trave.  ri.n  respect  to  the  i..Jec-f':  -puds, 
it  is  possible  to  determine  the  direction  of  the  velocity  perturbation  which 
increases  the  burning  rate.  From  this  information  the  sign  of  the  radial 
velocity  ir.csx.  for  a  given  radial  injector  arrangement"  can  be  found. 

Several  full  motor  tests  were  made  with  the  tangential  spud 
orientation.  Three  pressure  transducers  were  used  to  determine  the  direction 


of  the 

wave  trave! 

with  respect  to  the  injector 

spuds.  The 

results  of 

four 

typical 

1  tests  ace 

shown  in  tne  fol 

lowing  table; 

Test 

Chamber 

Injection 

^design 

Ji 

No. 

D i ameter 

0 i ameter 

•^design 

Mode 

Direction 

1 

9 

8 

1.4 

.90 

IT 

F  -*►  0 

2 

9 

6 

1.4 

.87 

2T 

F  0 

3 

6 

5 

2.2 

IT 

F  0 

4 

9 

8 

1.0 

1.39 

IT 

F  0 

The  wave  travel  direction  was  determined  as  follows.  The  injector  «as  so 
constructed  that  linking  from  the  nozzle,  tne  fuel  m '50! ion  hole  on  each  spud 
lay  on  the  clockwise  side  of  the  oxidizer  hole.  The  transducer  signals  showed 
that  the  wave  traveled  in  the  counterclockwise  direction.  Therefore,  the  wave 
crossed  each  spud  from  the  fuel  side  to  the  oxidizer  side,  or  F  — *r  0,  as 
shown  in  the  table.  For  the  first  ^ongentia!  mode  oscillations,  the  pressure 
transducers  were  arranged  at  90°  intervals,  as  are  No.  I  and  No.  2  in  Figure  10. 
In  the  case  of  the  second  mode,  the  wave  direction  was  determined  by  placing 
a  third  transduce  at  the  downstream  location,  but  angularly  midway  between 
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1ne  other  two.  fhe  resulting  phase  error  was  .ot  large  enough  to  affect  the 
test  results. 

The  angle  of  the  injection  spray  with  respect  to  the  combustion 
chamber  axis  is  a  function  of  the  ratio  of  the  test  mixture  ratio  to  the 
design  mixture  ratio,  A  design.  If  this  ratio  is  less  than  one,  the 
net  lateral  spray  momentum  is  in  the  directioii  r  -  s-  0,  and  if  jf  is  greafo- 
than  one  ,  the  spray  momentum  direction  is  0-^  F.  It  is  clear  from  the 

test  results  that  the  wave  travel  direction  is  not  affected  by  any  rotation 
caused  by  off-design  operaTion  of  the  injector. 

It  should  be  noted  that  although  th '  operating  conditions  for  tests  I 
and  2  were  very  nearly  the  same  the  resulting  oscillations  were  of  different 
modes.  Moreover,  om  Figure  ^6  it  can  be  seen  that  the  unstable  mode  of  the 
standing  type  is  the  second  tangential.  The  velocity  effects  for  this  con¬ 
figuration  are  illustrated  by  Figure  54.  Although  point  "A,"  representing 
the  test  motor,  is  within  only  the  second  mode  unstab 'e  zone  for  the  standing 
mode,  it  lies  within  both  the  first  and  second  mode  zones  for  the  spinning 
mode.  The  difference  between  test  I  and  2  was  in  the  starting  period.  That 
is,  test  I  utilized  the  normal  starting  procedures,  whereas  an  aluclnum 
barrier  was  used  in  test  2.  Tne  effect  of  the  barrier  was  to  stabilize  the 
first  tangential  tirade  during  the  first  second  of  operation  by  constraining 
a  standing  ""ode  osci  I  iat'.c-’  pattern  in  the  chamber,.  Upon  burnout  of  the 
barrier,  the  oscillation  changed  to  the  spinning  type  without  changing  the 
mode.  This  result  is  an  example  of  nonlinear  effects.  Although  for  very 
small  oscillations  both  modes  could  co-exist,  if  the  amplitudes  are  large, 
one  mode  will  dominate  to  the  exclusion  of  ■'ne  other.  In  this  case  the 
dominant  mode  will  be  the  first  tangential,  unless  the  second  tangential  mode 
is  allowed  to  grow  to  large  amplitudes,  as  in  test  number  2. 
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D.  Sector  Motor  Tosts  with  Radial  Toud  Orienlation 

The  full  motor  tangential  orientation  tests  demonstrated  that  for 
the  injector  used  in  this  investigation  tne  burning  rate  is  increased  by  a 
velocity  perturbation  in  the  direction  F- — ►  0.  For  the  radial  orientation 
angular  stab i I i  y  limit  tests,  the  arrangement  of  the  injoctor  was  such  that 
the  direction  0 — ►  F  was  pcsif’/c  (r'l+’-jard) ,  Therefore,  a  npga+i>-5  radial 
velocity  was  required  to  increase  the  burning  ratej  that  is,  the  radial  velocity 
index  Jlr  was  negative.  According  to  the  velocity  effect  theory,  as 
discusseo  ii>  Ciiapter  II,  for  jLf<.  C  ,  't'he  stability  iimits  on  the  /tv/C  “piane 
are  shifted  toward  smaller  'C  .  It  has  also  been  shown  that  the  stability 
limits  are  moved  in  the  direction  of  larger  f  as  the  sector  angle  Is 
increased.  Thus,  if  the  pressure  interaction  index  and  sensitive  time  lag 
are  to  be  independent  of  the  orientation  of  the  injector  spuds,  the  effect 
of  the  radiai  velocity  must  be  counteracted  by  an  increase  in  the  sector 
angle.  That  is,  the  srabiiity  limits  for  the  radial  orientation  should  be 
at  larger  angles  than  those  of  the  tangential  orientation  (no  velocity  effects). 

This  prediction  was  borne  out  by  all  of  the  radial  orienta+ion 
tests.  Figure  55  shows  the  results  for  the  9  inch  chamber  and  6  inch 
injection  diameters.  The  lower  limit  of  the  first  tangential  mode  wa^ 
located  (at  /i  -  1.4)  at  approximately  90°,  and  the  upper  limit  boiween  180° 
and  210°,  rather  than  the  50° -90°  and  180°  resuiis  of  the  previous  tangential 
orientation  series  (Figure  46).  From  Figure  55  it  can  be  seen  that  off-design 
operation  caused  significant  variation  of  the  stability  limits  with  mixture 
ratio,  which  may  be  due  to  distortion  of  the  recirculation  pattern. 

Similar  results  were  obtained  for  the  6  inch  chamber,  5  inch 
injector  (Figure  55).  For  this  case  the  shift  of  the  upper  limit  from  lf0°-iR0° 
to  2IC°  is  greater  than  the  shift  observed  for  the  9  inch  chamber.  However, 
one  hesitates  to  trust  the  test  results  which  show  the  first  tangential  mode 


at  angles  greater  than  180°  because  of  the  unknown  influence  of  the  unsteady 
vortex  pattern  resulting  from  the  flow  past  the  sharp  corner. 

The  experimental  results  for  the  9  inch  motor  v;ith  injection  Ci  a 
5  inch  diameter  circle  are  shown  in  Figure  57.  As  in  the  case  of  the 
corresponding  tangential  orientation  tests,  the  pattern  is  one  of  stability. 
This  result  is  in  agreement  with  the  theoretical  prediction  that  the  eftect 
of  the  raaius  of  injection  is  essentially  to  ciiange  the  interaction  Index 
required  for  instability. 

It  is  possible  in  principle  to  determine  the  msgnilude  of  the 
velocity  index  by  comparing  the  stability  limits  obtained  for  the  radiai  and 
tangential  spud  orientations.  Since  in  the  latter  case  the  velocity  effects 
are  absent,  the  difference  in  the  limits  is  a  measure  of  the  velocity  index. 
However,  the  rate  of  change  of  the  time  lag  with  respect  to  the  radial  velocity 
index  tor  a  constant  value  of  the  pressure  Interaction  index  is  very  small. 
Therefore,  because  of  the  large  uncertainty  in  the  stability  limits,  which  Is 
due  to  the  30°  Increment  in  sector  angle  obtainable  with  the  available 
experimental  apparatus,  it  was  not  possible  to  determine  the  value  of  the 
velocity  index  from  these  radial  orientation  tests  with  any  meaninyfu! 
accuracy. 

E.  Conclusions 

The  primary  g.,-:;  of  this  investigation  ^as  to  evaluate  the  validity 
of  applying  to  the  transverse  modes  of  instability  the  sensitive  time  lag 
theory  which  had  been  used  successfully  to  explain  the  longitudinal  modes. 

In  that  theory  the  combustion  process  was  characterized  by  twc  pai  ameters, 
the  sensitive  time  lag  tT  and  the  inter .iurlon  Index  ^  .  It  was  assured 
that  all  physico-chemical  conditions  which  affected  the  combustion  process 
rates  could  be  correlated  to  the  Iccal  pressure.  The  degree  of  sensitivity 
was  measured  by  the  Infaraction  index.  The  longitudinal  mode  analysis 


showsvj  that  if  the  interaction  ind4.<  ^or  a  given  rocket  motor  system  was 
sufficiently  large,  !i>ere  existed  ranges  of  the  dimensionless  tinie  . 
lag  !  {_*  in  which  the  rocket  operation  would  be  oscillatory.  The 

predictions  of  this  theory  were  verified  for  a  variety  of  injector  types 
and  propellant  coiiib  inat  ions  by  means  of  stability  limit  measurements  using 
a  var  i  ab  1  e  I  ength  (  )  rocket  motor , 

When  the  sensitive  time  lag  theory  was  applied  to  purely  tranr-verse 
modes  of  instabill+y,  similar  results  were  obtained  as  for  the  longitudinal 
case  wiih  rne  combustion  chamber  radius  ^  replacing  the  t-liamber  length 
in  the  nondimensional  ized  time  lag,  /fc  •  However,  a  preliminary 

series  of  exploratory  tests  showed  that  the  pressure  sensitivity  theory  was 
not  sufficient  +o  cxp!a*n  two  experimental  observations:  (1)  the  stability 
behavior  of  the  standing  forms  of  tangential  modes  was  quite  different  from 
that  of  the  spinning  form,  and  (2)  the  orientation  of  the  doublet  injection 
pattern  with  respect  to  the  radial  and  tangential  coordinate  directions  had 
a  significant  effect  on  stability. 

Based  on  the  findings  of  the  preliminary  test's  a  velocity  sertsJ+lvi'ty 
theory  was  developed.  The  radial  and  tangential  velocity  perturbatio  .s  were 
assumed  to  influence  fhe  rates  of  the  processes  involved  i  nr,  the  .combustion, of 
the  liquid  propellants  in  a  manner  analogous  to  the  pressure  perturba  Ion 
in  the  original  t'rory.  A  physical  mechanism  was  s^d9®sted  to  lie  in  the 
enhancement  of  the  mixing  o.  the  vapors  of  the  two  propellants  which  could 
be  effected  by  the  transverse  components  of  the  gas  velocity.  The  analysis 
based  on  these  considerations  predicted  that  the  tangeni iai  velocity 
perturbation  would  have  a  strong  dcitati I izing  effect  on  spinning  tangential 
modes  but  would  have  no  effect  on  standing  modes.  The  radial  velocity  was 
shown  to  cause  a  shifting  of  Hie  unstable  time  lag  ranges  for  a  given  inter¬ 
action  index.  1h-<  orientation  effects  noted  in  the  preliminary  experiments 


-1  16- 


were  seen  +d  be  the  results  of  the  uni-direc .ionai  sensitivity  of  the  doublet 
type  of  injector 

Previous  instability  anaiyses  had  assumed  that  the  steady  sta+e 
flow  in  the  combustion  chamber  was  one-dimensional,  implying  that  the 
propei lants  were  injected  uniformly  across  the  face  of  the  injector.  However, 
practical  considerations  dictated  the  use  of  rocket  motors  in  which  the 
propellants  were  injected  around  the  circumference  of  a  single  c'rcle.  Another 
result  of  the  preliminary  test  phase  was  +Hat  the  radius  of  the  injection 
circle  in  reiation  to  the  combustion  chamber  radius  was  an  important  factor 
in  the  stability  behavior  of  a  rocket  motor.  The  effect  o!  the  nonuniform 
injection  di stirbution  was  treated  in  the  present  analysis  by  means  of 'a 
trans -erse  distribution  function.  For  injection  concentrated  at  a  single 
radius,  it  was  oemonsirated  that  there  is  a  greater  tendency  to  Instability 
as  the  injection  radius  is  increased. 

In  order  to  evaluate  the  instability  theory  in  a  systematic  manner, 
it  was  necessary  to  find  a  method  of  measuring  stability  limits,  as  was  done 
for  the  longitudinal  imode  case.  For  tangential  modes  of  oscillation,  the  sector 
motor  mehhod  of  measuring  angu'ar  limits  was  developed.  In  this  metnod,  the 
sector  angle  replaces  the  chamber  length  as  the  variable  which  governs  the 
resonance  properties  of  the  combustion  chamber.  The  requirement  imposed  on  a 
atebiiif/  llmii  motor  is  rnat  the  recirculation  pattern  be  preserved.  Ji  the 
case  of  the  sector  motor  this  implies  that  the  injection  distribtuion  must  be 
periodic  tn  the  tangential  direction  and  that  the  propellants  sprays  must  be 
directed  a'^ioMy  In  order  to  achieve  proper  simulation  of  the  full  circular 
ciiamber.  A  limitation  of  this  method  was  seen  to  be  present  in  the  case  I'hat 
the  time  lag  of  the  rocket  system  is  such  that  the  first  tangential  mode  is  un¬ 
stable  for  sector  angles  greater  that.  180°.  For  such  angles,  the  oscillating  flow  past 


the  corner  would  set  up  an  unsteady  vortex  pattern  with  unknown  e+fects  on  the 
stab i I ity  1 imits. 

From  the  measurements  of  angular  stability  'imits  made  for  several 
combustion  chamber  conf iguraf ions,  a  qualitativi,  but  not  quantitative, 
confirmation  of  the  theoretical  predictions  concerning  velocity  and  Injection 
distribution  effects  was  obtained.  Tne  interaction  Index  values  derived  tr''m 
the  angular  stability  limits  <-e.  >  so  low  that  longitudinal  mode  instability 
was  not  possible.  Tests  made  with  a  variable  length  30°  sector  motor  did  show 
no  longitudinal  mode  instability,  precluding  a  quantitariva  comparison  of 
combustion  parameters  derived  from  the  two  different  methods.  Because  or  the 
rather  large  sector  angle  increment  (30°)  It  was  Impossible  to  derive  values 
of  the  velocity  interac+ions  indices. 

The  experimental  phase  of  the  investigation  was  not  intended  to  be 
exhaustive.  Tests  were  confined  to  a  single  injector  type  and  one  propellant 
combination.  However,  with  the  limited  number  of  configurations  explored  it 
was  seen  that  the  combustion  parameters  were  strongly  dependent  on  the  Injector 
spray  characteristics  and  on  the  nature  of  the  recirculation  pattern.  Therefore, 
there  is  a  clear  need  for  systematic  investigations  into  the  factors  v;r  ich 
influence  the  sensitive  time  lag  and  interaction  Indices  (both  pressu.'e  and 
velocity).  A  possible  method  for  measuring  the  velocity  parameters  is 
illustrated  by  Fit, -re  58.  The  propellants  are  injected  radially  Inward  from 
the  chamber  wall  in  a  variable  length  rocket  motor.  For  this  manner  of 
injection,  the  velocity  fluctuations  associated  with  longitudinal  modes  will 
act  transversely  to  the  propellant  sprays.  If  the  injector  section  is  placed 
near  the  closed  end  of  the  chamber,  me  pressure  effects  will  dominate,  whereas 
location  of  the  injection  near  the  center  will  maximize  the  velocity  effect. 

By  measuring  stability  limits  for  several  injector  positions,  it  should  be 
possible  to  deduce  values  of  the  velocity  parameters. 
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A  tactor  wnkf,  has  not  been  consMered  in  this  investigation  is  the 
possibiiity  of  an  osciilating  injection  flow  rate.  For  chambers  of  the  sizes 
used  in  the  experimental  tests  the  oscillation  frequency  is  sufficiently  high 
that  the  feed  system  will  not  respond  appreciably  to  the  pressure  fluctuations 
in  the  combustion  chamber.  However,  as  the  diameter  is  increased  the  frequency 
is  reduced,  with  the  rciuir  tnat  ^or  la-gc  — ci'ot  motors  the  in'^crcof io,-- 
the  feed  system  with  the  combustion  chamber  can  no  longer  be  neglected  in  the 
analysis  of  high  frequency  instability.  In  such  problems  the  analyses  of 
Scala  (5)  and  Croke  (29)  should  prove  helpful.  However,  neither  investigator 
con‘'idered  wave  effects  in  the  feed  lines,  which  have  often  been  observed  to 
be  of  considerable  importance.  Thus,  aaditional  complications  due  to  adding 
injection  rate  fluctuations  will  be  involved  in  the  application  of  the 
instability  theory  to  large  liquid  rockets. 

In  the  theoretical  developments,  the  through  velocity  in  the 
combustion  chamber  was  restricted  to  values  sufficiently  small  that  the 
square  of  the  Mach  number  could  be  neglected  compared  to  unity.  The  present 
theory  might  be  carried  to  0.3,  corresponding  to  the  Redstone  engine, 

out  at  tnis  point  terms  of  the  order  of  ten  percent  are  being  neglectea.  In 
order  to  apply  the  analysis  to  more  recent  engines,  the  restriction  on  the 
velocity  will  have  to  be  removed,  thus  adding  a  great  deal  of  compi«;,>ity  to 
tne  prooieiii. 

The  next  major  step  in  the  investigation  of  combusi ion  instabiiity 
must  be  the  consideration  of  nonlinear  effects.  All  of  the  systematic, 
quail Mta five  work  to  date  has  been  concerned  with  linear  instability,  that 
is,  with  oscillations  which  are  spontaneously  initiated.  However,  even  some 
of  the  experiments  made  in  this  investigation  indicated  the  presence  of 
nonlinear  phenomena,  particularly  in  the  starting  transients.  Some  types  of 
injectors  (for  example,  the  showerhead  and  i ike-on- 1  ike)  are  clearly 


A  iso  the  shattering  of  liquid 


susceptible  to  nonlinear  veiocitv  e.-^ects. 
pnopellant  droplets  by  high,  amplitude  waves  has  been  shown  to  exert  an 
important  effect  on  the  combustion  process.  At  the  present  time  The  theo¬ 
retical  approach  to  nonlinear  instability  has  barely  begun.  In  addition  to 
more  sophisticated  theoretical  techniques,  "’uch  more  refinea  instrumentation 
and  experimental  methods  will  be  required.  However,  the  linearized  solution 
0+  the  instability  problem  has  yielded  much  valuable  insight  into  the 
combustion  process  in  liquid  propei.'ant  rockets  and  should  provide  a  firm 
founddtio.,  for  the  future  studies. 


I 
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TABLE  : 


ResuHs  of  Preliminary  Barrier  Tests 


Test  No, 

Cor.f  iguraTion 

Stab  1 1 ity  Behavior 

( 

Pull,  6”  chamber, 

f?*  injector,  spudr.  tangenfiai 

185  psi  peak-to-peak 
osc i 1 1  at  ions 

2 

IcO'^  sector  (insert) 

Stab  1  e 

3 

Barrier  extending  from  injector 
face  to  nozzle  throat 

Stable 

4 

Barrier,  from  injector  to 
l-3/eP  downstrecn 

Stab  1 e 

5 

Barrier,  from  injector  to 

3/4”  downstream 

Stab  1 e 

& 

Barrier,  from  injector  to 

5/8P  downstream 

Stab  1 e 

7 

Barrier,  from  injector  +o 

I/2P  downstream 

40  psi  peak-to-peak 
osci 1 1  at  ions. 

8 

Barrier,  from  injector  to 

I/2P  downstream 

90  psi  peax-to-peak 
osci 1 1  at  ions 

9 

Barrier,  from  injector  to 

1/Z*  downstream 

Stab i e 

10 

Barrier^  from  i/4"  from  injector 
to  l-5/fi?*  ..o'-n stream 

17  psi  oeak-tc“peak 
osci 1  let  ions 

Nwtes.  I.  The  same  injecfor  and  chamber  were  used  on  all  tests. 

2.  All  instabilities  were  of  the  first  tangential  mode. 

3.  The  barriers  showed  evidence  of  burning  on  tests  7  &  3. 
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TABLF  II 


Velocity  Effect  Coetficients  for  Uniform  Injeciion, 
Full  Circuiar  Chamber 


^#,(Spinnin 

Mode 

h 

Sifi, 

mode) 

F  i  rst  tangent i a  1 

1 

1 

0.  |9i 16 

0.77402 

Second  tangent. c 

2 

1 

0.18124 

0.CS253 

F irst  raoiai 

0 

2 

A  1  ?  7  -^A 
—U  .  J 

0.00000 

F  i  rst  comb i ned 

1 

2 

0.02884 

0.22197 

JPR  261 


Transverse  instaMity  motor  schematic 


FiQure 


Back  of  spud  injector 


Transverse  instability  tests 
tangential  spud  orientation 


JPR  446 


Key  !  X  stable  o  second  tangential  mode 

•  first  tangential  mode  •  mixed  IT ,  2T,  I R  modes 


Figure  if 


(a)  Barrier  tests  (b)  Sector  tests 

mporary  and  permanent )  (O  to  360*  30*  increrrents  } 


figure  12 


FlQure  /. 


Exploded  ^lew  of  sector  motor 
using  distributed  injector 


Tangential  orientation  txtrrier  tests 


Figure  14 


Radial  orientation  barrier  tests 
'chamber ,  8" injector ,  iSOpsia,  500  lb  thrust) 


Figure  15 


Mixture  ratio 


JPR  838 


(a)  Spray  produced  by  tangentially  oriented 

injector  spud 


(b)  Spray  produced  by  rotated  spud 


Figure  17 


Rea!  and  imaginary  parts  of  A. 


First  admittance  coefficient,  A 
u  -  0.10,  linear  nozzle 


fO 

CD 


(C 

& 


Figure  Id 


Typical  solution  for  /yu(aj)  and  t((u} 


Figure  20 


JPR  845 


T/feof  atical  stability  limits  for 
several  purely  transverse  modes, 
uniformly  distributed  combustion, 
u^  =  0.10 


U 
i  i 


w 


77 


first  tangential 
second  tt  ngentiol 
first  radh  U 


T*C^ 


Figure  21 


./  2 
Chamber  exit  velocity, 


Figure 


TheoretiCQi  effect  of  chamber  exit  velocity  on  stability  limits, 
first  tangential  modes  uniformly  distributed  combustion 

/S  =  070s  jc=a92,  ^  =1.84 


u^=.05 - 

/  / 

to - 

/  / 

.20 - 

'  /  ' 

>  /  ! 

1 

1  // 

»  §  • 

Figure  23 


Conical 

Linear 


Theoretical  effect  of  /3  on  stability  limits, 
first  tangential  mode,  uniformly  distributed 
combustion 


u.-a!0,  c=a92,  ^-1.84 


Theoretical  radial  velocity  effei 
limits,  first  tangential  mode, 
distributed  comd 

We  =  0.10 


JPR  869 


Presiure  coefficient  A 


Figure  30 


Radial  vetc  city  coefficient 


for  combustion  concentrated  at  r  =  rf 


Figure  32 


Tangential  velocity  coefficient 


for  combustion  concentrated  at  r  =  r^ 

- 1 - p- - 1 - - 

—  First  tangential  mode 

—  Second  tangential  mode 


JPB  712 


Figi 


Vcii'iciblG  ongiG  sGCtor  inotof  ^ 
ith  radially  concentrated  injection 


JPR  732 ‘ 


Acoustic  frequency  number  Sy^ 


JPH  454 


Variable  -  length  sector  motor 


JPR  878 
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Figure  44 


Standing  mo  'e  wave  shapes, 
first  tangential  mode, 
160**  sector  motor 


(a)  Pressure  traces 


(b)  Location  of  pre  "iure  transducers 


Figure  45 


Mixture 


Tangential  orientation  sector  tests 
9"  chamber,  8"  injector,  Kdesign'  14, 

Pg  (nominal)  -  ISO  psia,  F  (nominal)  =  iOOO  lb 


Key:  %  IT  O  2T  €>  IT  a  2T  mixed 
X  Stable  «  n  intermittent  instability 


0  60  120  180  240  300 


Sector  angle,  degrees 


<0 

oc 

a. 


Figure  46 
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Theoretical  sector  motor  stability  limits, 
first  tangential  mode ,  9  "chamber  diameter, 
8"  injection  diameter ,  p.  =150 psia,F-IOOOIb 
no  velocity  effects 


0  12  3 

T 


Figure  47 


Longitudinal  mode  stability  tests 
30^  sector^  =1.4,  tangential  orientation 


X  =  stable 


Theoretical  sector  motor  stability  limits, 
first  fangeritiol  mode  ,6  chamber  diameter 
5“  injection  diameter  ,  p  =l50psia,F=IQ0Qlb 

w 

no  velocity  effects 


ldO\ 


S0\  90\  I50\ 

1  1  V 


mx  range  from 


Figure  50 


Tangentic.  orientation  sector  tests 
6"  chamber,  5”  injector,  e design  *  f 
(nominal)  ~  150  psia,  F  (nominal)  - 1000  lb 

Key  -  %  IT  O  2T  €  IT  3  2T  mixed 

X  Sfnhie  ■  n  intermittent  irztabiUty 


JPR  083 


Theoretical  sector  motor  stability  limits, 
first  tangential  mode ,  9"  chamber  diameter 
Subjection  diameter ,  p  =150 psia ,F=IOOOIb 
no  velocity  effects 


I 


range  fn 
Subjection  di 
tests 


Figure  52 
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Jangeniiol  c  ientafhn  sector  tests 
9  chamber,  5"  injector,  ^design  =  14, 

Pc  (nominal)  -  ISO  psia,  F  (nominal) = 1000 lb 


Key  •  m  IT  O  2r  e  IT  a  2T  mixed 
X  Stable  ■  u  intermittent  instability 


Sector  angle,  degrees 


Figure  53 


Theoreti''al  tangential  velocity  effect, 
combustion  concentrated  at  r.  =  8/ 9, 
=  0.05 

-  no  velocity  effect 

- ^shcity  effect  included 


Radial  orientation  sector  tests 
9“  chamber,  8"  injector,  Kdesign' 

Pc  (nominal)  -  150  psia,  F  (nominal)  - 1000 lb 


Key  :  %  IT  O  2T  ^  IT  Q  2T  mixed 

X  Stable  a  n  intermittent  instability 


0  60  120  180  240  300 

Sector  angle,  degrees 


O 

«0 


a: 
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Figure  55 


Radial  orientation  sector  tests 
6“  chamber,  5"  injector,  Kdesign  - 
pg  (nominal)  =  150  psia,  F  (nominal)  =  1000  lb 

Key  '  %  IT  O  2T  €)  IT  &  2T  mixed 

X  Stable  ■  n  intermittent  hstability 


tf) 

7> 

(t 

a 

“5 


Figure  56 


Radial  orientation  sector  tests 
9"  chamber,  5"  injector,  n. design  ' 

Pp  (nominal)  -  i50  psia,  F  (nominal)  - 1000  lb 


Key:  %  IT  Q  2T  e  IT  a  2T  mixed 
X  Stable  ■  u  iniermittent  instability 


3.0 


2.6 


2.2 

•9 

§  1.6 
.Js 

^  ,A 


1.0 


0.6 


) 

X 

m 

* 

\ 

> 

\ 

m 

1 

D 

i 

X  D 

0  60  120  160  240  300 

Sector  angle,  degrees 


a> 

tf) 

00 


cr 

0. 


Figure  57 


Amplitude 


l£ 
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Figure  58 


Velocity  Effects  Motor 

Variable  length  liquid  rocket  motor,  radial  iniection 


OJO 


u,Ul  0.05 


1.05  , 

...  ^ 

- 5r 

_ 

— 

a 

.2  .4  .6  .8 

axial  distance^  z 


10 


Steady  state  liquid  droplet  velocity  profiles 


Figure  A! 


Linear  nozzle 


Comparison  of  nozzle  profiles 
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Figure  Q 


0.25  0.50  07i 

Distance  from  entrance 
Subsonic  iength 

Comparison  of  subsonic  veiocity  profiles 


for  conical  nozzle 
0 


Range  of - 

interest 
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A  PENDIX  A 

Liquid  Oroolet  Dynamics 


The  steady  state  droplet  aynamics  equation  describing  the  inter- 
cfiange  O'  momentum  between  the  liquid  and  gaseous  phases  has  been  written  as 


Jii 

where  is  the  momer'um  interchange  coefiicient.  An  approximate  expression 
for  v-a.i  ne  obtained  by  consiaer:ng  a  non-burning  spnerica!  droplet.  The 
principle  of  conservation  of  momentum  yields 


(2.5e) 


(A-l) 


where  fi*  's  the  liquid  density  (dimensional),  ol^  the  droplet  diameter, 

ru 

Cp  the  drag  coefficient,  a^and  the  jas  and  liquid  velocities,  thus 


(A-2) 


which  can  be  written  in  terms  of  the  Reynolds  number. 


- 


as 


0  3  /*  O 

*  - 

For  Reynolds  numbers  in  '■he  range  of  rocket  motor  interest  («o<^<iCO  ), 
is  nearly  constant.  In  view  of  the  complexity  of  "Ihe  notion  of  burn¬ 
ing  propel  lar-t  droplets  under  turbuis.T  conditions,  it  is  probably  nest  to 
regard  the  coefficient  as  an  empirical  factor  correlating  the  liquid  and 
gas  velocity  variations.  In  dimensionless  form 

- W 


{A-3) 


A-2 


A  soluticn  of  equation  (2.5e)  can  be  obtained  analytically  by  making 
the  substitution 

which  fcr  3  linear  variation  U  wit*'  g  {i.e.,  U>=C|2  )  gives  the 
integrabie  equation 

The  result  is  th '  cumbersome  expression 


4. 


A  more  convenient  procedure  is  to  solve  equation  (2.5e)  numerically.  For  most 
purposes  the  simple  difference  equation 

with  the  boundary  condition 


yields  sufficiently  good  results.  Figure  A  I  shows  typicai  solutions  for 
different  profiles  2r,d  -fc,  values. 

The  equation  for  the  perturbation  axial  liquid  velocity  is  the 
first  of  equations  (2.14f), 

and  the  boundary  condition  is  Sea  I  a  (5)  has  considered  the 


(A-7) 


solution  of  fhis  equation  in  detail.  His  result  can  be  obtained  by  making 


A-3 


an  order  or  magnitude  a.alysis.  Assuminc  that  the  droplet  momentum 
coefficient  and  the  liquid  droplet  velocity  are  0(  tT*  ),  as  discussed 

above  in  Chapter  II,  consideration  of  the  steady  state  droplet  equation  (2. be'/ 
shows  that  is  also  Of  )-  From  equation  (A-7)  it  can  be  seen 

that  U-L  !s  at  must  Of  ).  Then,  assuming  that  differentiation  Joes 

not  change  the  order  of  magnitude,  i.ic  CL  and 

^  d*  *• 

—■2.  — 

equation  CA-7)  are  Of  Ug  )  while  the  other  terms  are  Of  ).  Therefore, 

the  solution  good  to  terms  of  Of  Ue  )  is 

u:  ^  :^uf  -h  O  (Si)  CA-8) 

However,  tn^s  order  uf  magnitude  procedure  cannot  be  used  in 
gene'al,  since  it  may  not  oe  possible  to  satisfy  the  boundary  conditions 
when  the  term  containing  the  highest  derivative  is  neglected.  Thus,  for 
the  transverse  velocities  I'he  relations 

'0:'  sz  'V''  5=  —  ur' 

car  not  satisfy  the  boundary  conditions  at  the  injector  face,  where 

ur^  =  0 

but 

¥=-o  ^  ^  o 

For  either  the  radial  or  tangential  liquia  velocity  perturbations, 
i  the  governing  equation  is 

S -U-J  +  ^  ’fe.  (qr' -/O-J  ) 


which  can  be  written  as 


An  integrating  factor  for  this  equation  is 

e 


Multiplying  through  by  this  factor  gives 

“‘Sj'  ^  A-rS)  \.l 


nr‘e 


r  Ji' 

o 


It  can  be  seen  that  the  left  side  is  equivalent  to 


3  r  ,  , 

—  !  qT  (O  'i  \ 

3*.  I  ■ 

L 


Thus  integrating  both  sides  with  respect  to  2  yields 


‘’i  j  ^  j./ 


=  -ltl  LSi'e' 

J 


(A-iO) 


The  right  side  of  this  equation  can  be  evaluated  by  inteoratinn  by  parts,  so 
that  equation  (A-!0)  becomes 

'Ole  =  ^  d'^ro)  -I  »  d2  I 


(A-l I) 


It  is  shown  in  Chapter  II  that  is  0(  -v  ).  Neglecting  the 

integral  term  on  the  right  side  gives  an  expression  for  the  liquid  velocity 
perturbation  which  is  good  to  0(  iSe  nr'  )  but  not  to  0(  w?  ): 


'  4 


=  r 


^  V(t^oe)  ~nr'(c^nS)e 


-''W  t' 


(A-12) 


A-5 


Since  i=  0(  uo  Sr'  )  and  z  is  0(i),  fhen  •'/'(«)  ~  ■v(o}  is  ), 

or  \r’(o)  -  ArV^rj+Of%v-).  Also  -rOCtfe^,)  .  Therefore,  the  final 

expression  for  'trj  can  be  written  as 


(A“I3) 


where 


4.‘w  =4(1- 


o 


J 


(A~I4) 


For  nei-.f ' ■  O'C  i '  ia+ions,  s^luj  ,  and  for  s'lial  !  va!"ep  ot  the 
droplet  momentum  coefficient,  ^  ,  Then 


^(ir) 


c 


ti* 


This  function  is  represented  by  a  curve  on  the  complex  plane  which  for 
increasing  axial  distance  2  spirals  inward  around  the  point  (1,0).  As  an 
example,  the  function  is  shown  in  Figure  A2  for  O-OS  , 

l<lw,=  0.oS  ,  .  (the  dimensionless  frequency  of  the  first 

tangential  acoustic  mode),  and  0  S-  ?  ■S  O-S  . 
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APPF  DIX  B 

Peiturbation  Solution  with  Axially  Concell^^ated  Combustic'' 

When  combustion  is  concentrated  near  a  single  axi^i  station,  as 

gIu. 

in  the  case  of  well  "'ixed  propellants,  the  steady  state  velocity  gradient 
is  locally  0(1).  The  analysis  presenteo  !.i  Chapter  II  assumes  I'he 

velocity  gradient  never  exceeds  0(  "--c'  ).  The  case  of  axially  concentrated 

combustion  requires  a  somewhat  different  and  more  complicated  procedure  than 
that  us«d  p''e'.  iousi  y . 

Tie  first  order  perturbation  equations,  with  ^=^0(0  ,  are 

^  o  4.  ^  1  4-  p  =  (B-I) 

1  1  32-  a  r  r  r‘  a  e  o?  i  /  dt  J 


S  u,  -f- 


4*2  /o 


(B-2) 


s i'/ir  f- o  )  j_  4^  -v-  -  4-"*’  \r  = 

s’  °  ‘  Jz  i  •  ®  >r  9  ^ 


(B-3) 


-  -(.i 


dS  , 


'jr_  -r 


pof;  =: 


-± 

'h'r  3  9 


-f,)  ^  si  h)  -i- ^  ^ >1 


(B-S) 


The  entropy  perturbation  is  given  by 


>  =  ^  -r: -  •  (%  + 


(B-6) 


B-2 


At  Ths  chambe"*  e^it,  where  3^  ^  ^  =  O 


,  as  before. 


The  energy  and  contineiTy  equations  are  combined  to  eliminate  the 
density  jo^  .  Then  1he  set  of  linear  partial  differential  equations  to  be 


so ! ved  i s 


“>  +  sl'  +  ^  +  tt'  +J-?2!  =  0  -i“  t 


If?)  *  =  -  2  K  -25  ^(1") 


S  iaT'  -I _ L— —  Ur  ^  ?L?  uf  —  p,  *  ur 

'■  <r2&\V  /  oil  °  ol2  ‘  “ 


(B-IO) 


As  was  done  previously,  the  solutions  ai e  represented  as  sums  of 


partial  solutions. 


-f-^C  +>J> 

U.i  ■=.  +  '^  + 

■*■  Qs  ■*■  Qc 


where  is  to  be  identified  with  the  pressure  depenuencp  Og  'lith  the 

radial,  and  witn  the  tangential,  veloc'+y  deoendence. 

The  right  hand  side  of  equation  (B7)  can  be  simplified  by  noting 


The  otner  equations  can  be  treated  similarly. 


B-3 


'U 


The  equations  for  the  pt'rtial  solution  "A"  are 
'Jfi. 


^  +  5  ( 'i  t-  ^  4-  4  '  — ^  -  CiA  -  - 

33-^  ^  Y  J  3  r  ^  r  r-  9  Q  ^  ' 


2* 


'JiV  t‘)  + 

-(  —  )-  -'UV;) 


'U"  4*  - 

3r 


The  equations  for  the  ”B'*  and  "C**  partial  solutions  are  of  the  same  form. 
For  the  "D'’  soiu+'on,  the  equations  are 

I?  4  S  (t»j  4  4  4  i  If 

2{%)  =-2|,;=^tO 


S'Ts  4  -rr('7  j  =  -4'ft  n; 


Consider  first  tne  set  (BI3).  Since  the  right  side  of  each 
equai ion  contains  only  zeroth  order  terms,  the  "D"  solutions  take  the  form 

P,w  , 

A 


'V'p  -  Vp  {%)  (©) 

1® 


-:bs'*'itution  of  these  expressions  1  ito  (BI3)  gives 

^(Up+yu^)4-S^  —  =  o 

4(1+^445)  4SUp  =o 


sV,  4  S 


(^-12) 


(n-14) 


(B-15) 


B-5 


•?-  2' 
r  r 


5  =  s  '  +s'^i  E '(i-'Jife" -i  B?J^  (B- 

j  j  y 

o  o  o 

since  4  5io)=o  B  iO)-=  .  Thus,  +he  result,  good  to  oC<iij  .  is 

^  =  sfy+i)J  -f  sf j 

Ho  “i  J,,  -i  ^ 

d  O 

v,'hicn  is  identical  to  (2.50). 

Thf.  soi'jticns  “A",  "B“,  and  "C"  are  slightly  more 

compi icated  because  the  terms  on  the  right-hand  side  of  the  eQaatlons  are 
first  order  quantities  rather  than  2eroth.  order. 

All  three  solutions  follow  the  same  procedure;  the  "A"  solution 
will  be  shown  in  detail. 

The  equttio'.s  fo  c.  cati-.'ltd  —  c  I'.e  i»:  or  linear,  non-homogeneous, 
simultaneous  partial  dtf ferenticl  equations  (B''I2),  The  pressure,  veloci+ies, 
and  burning  rato  are  expanded  in  The  transverse  eigenfunctions 

t.  =  £  f  A„  Pa„  (!)  i„(')  14 

P*«  5=  I  'j 

f»0  <PI  \  {Q. 

fto 

pro  2^:ri  ^  *  * 

Then,  +or  each  p  q  pair,  the  equations  to  be  solved  are 

S  (  ^‘‘ri  ^  )  4  s  _s^  y,„  =  c  ') 


(B-22) 


(0-23) 


(3-24; 


s  v„„  +-  +4 


(B-25' 


<  oni  the  first  equation  by  means  of  ihe 


EliminatlriQ 

third  of  equations  (B-25)  gives 


^  ^  ^  ^'Arf  "  ^ 

By  analogy  with  the  '*0"  solution,  ne.;  oe^-ndent  variables  are 


Bp»  =  "If  +  ?  s  U, 


In  terms  of  these  variables,  equations  (B-26)  can  be  written  as  before. 


-t  S  (  1+  =•-  S  E 


t*  ^  r<i 


^  S 


where 


^P?  =^-(H^)2.U  L/Ay<J 

R*}  “  ^  ^AfKj  —  Hap<?  +^cr\4p<} 

Solving  for  Sp,-  in  equations  (B-28)  gives 

1  B  i=r  5 

—  tS  ^SpCfJ  OpQ  —  tp-J 

where  the  bounac,  conditions,  as  for  the  "O'*  solution,  are 


(B-26) 


(B-27) 


(B-28) 


(6-29) 


CB-30) 


For  •  ■i’h®  solution  is  obtained  by  the  method 


of  variation  of  constants  as 


a  ^  S.n(.  ,  P 


(B-31) 


B-7 


However,  the  soi'..rion  of  iniorest  in  determini  jg  the  stability  criteria  is 
that  for  ,  c^=  K  ,  for  which  =C>(ie.)  •  Then  the  equation 

for  Byl,  is 


-h 


s^E 


/ 


where  to  0(  Ug 


),  the  functions 


are 


^  Fvj,  '  -  fHAyt 

4  I 

Now  HavI,  is  defi  'SO  by 

J 

and  neglecting  higher  order  terms, 

H  iZ  P«o 

so  that 

i 

^  +  B^u  (o  ) 

o 

or 

_'il  ru^2')cl2'  + 

Th'S  is  the  same  result  as  presented  in  equations  (2.49)  and  (3.12a).  The 
boiuiions  for  and  exhibit  the  same  behavior,  s'  jwi.-.g  that 

although  the  velocity  gradient  ^  may  be  locally  0(1),  the  overall 
effect  is  the  same  as  if  it  were  of  the  same  order  as  the  velocity  itself. 

i 


B-32) 
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AF  ’LNDIX  C 


Noz2l©  Adcii'M'anc©  Co©f  f  icisni's 


Jn  Chapter  III  a  relation  between  the  pres'^ure,  velocity,  and 
entropy  perruroations  at  the  entrance  of  the  nozzle  was  Introduced  as 
the  boundary  condition  at  fhe  nozzle  end  of  the  combustion  cr.eiiioer.  Thi-‘ 
-"admittance  relation"  for  neutral  oscil.ations  in  a  supercritical  nozzle 
was  derived  by  Crocco,  first  for  longituoinai  modes  of  oscillation  (Refer¬ 
ence  i.  Appendix  B),  and  later  extended  by  him  to  transverse  mooes.  A 
brief  description  of  the  latter  derivation  is  given  in  Reference  25.  The 
purpose  of  this  Appendix  is  to  discuss  the  calculation  of  the  admittance 
coefficients  and  to  present  the  results  obtained  for  two^yplcal  cases. 

Tne  transverse  admittance  relation  for  axisymmetric  flow  is 

p  ^  p  t"  ®  ^  4-  c  ^ 

too  R>,.  R), 


whe'  e 


a  _vg, 

(3=  t 

^  ^  'V®  (  ft  ^  5^  ) 


-v'u. 

are  the  admittance  coefficients  (the  not».tion  here  differs  somewhat  fo 
that  of  Reference  25).  In  the  expressions  (C-2),  u  ,  C  ,  OJ  ,  and  'i* 
refer  to  steady  state  gas  velocity  and  sound  velocity,  oscillation  fre¬ 
quency,  and  specific  heat  ratio,  respectively,  just  as  in  Chapter  HI, 
although  a  different  nondimensional ization  scheme  has  been  used.  The 
auxiliary  functions  %  .  \  ,  and  ^  wi 1 1  be  discussed  below. 


(C-l) 


(C-2) 


C-2 


The  independeni  v=5ric\ble  in  the  an?' lysis  of  the  oscillatory  nozzle  flow 

was  taken  as  the  steady  state  velocity  potential  ,  defined  by  the  equation 


(L  (i) 


ci 


(C-3) 


The  velocity  profile  in  the  nozzle  is  determined  by  the  converge''^  section 
geometry,  although  it  has  bee-  convenient  in  fhis  calculation  to  find 

a  convenient  analytical  expression  for  iTf?)  which  approximates  well  the 
actual  velocity  variation,  which  can  be  determined  only  by  numerical  methods. 

~'ie  admittance  coeff icie^ts  may  be  obtained  for  an  entire  family 
of  nozzles  at  one  time.  This  fatriily  consists  of  all  nozzles  obtained  from 
a  single  nozzle  by  a  change  of  scale  in  the  axial  direction.  That  is, 
is  specif iea  rather  than  ‘Jli)  ,  where  K  is  an  arbitrary  multiplicative 
constant,  and  it  is  clear  that 


oljW) 


The  parameter  of  affinity  K  can  be  defined  in  many  ways.  In  the  present 
calculation  it  has  been  taken  as  the  (dimensionless)  velocity  gradient  at  the 
throat,  i.e.. 


(C-4> 


Then  the  following  transformation  can  be  made  in  order  to  inf'lutie  al !  nozzles 
of  the  family  in  the  calculation  of  the  coof ic  ients: 

X  ~  ZK 

K'  - 


(C-5) 


C-3 


Insertion  these  expressions  inio  the  definitions  of  the  admittance 
coefficients  gives; 


A  = 


)  -Cf' 


(R 


ioj'  ^  ' 

-  il.  V£-A. 


40' 


c  - 


i/i,  = 

-  uf' 


AN  quantities  in  the  nozzle  analysis  have  been  nondimensional ized 
by  corresponding  values  at  the  nozzle  throat  section;  lengths  by  the  throat 
radius,  pressures  by  the  throat  (critical)  pressure,  velocities  oy  the  sonic 
velocity  at  the  throat,  etc.  In  order  to  match  the  nozzle  conditions  to  the 
corresponding  combustion  chamber  conditions  at  the  entrance  of  the  nozzle, 
the  following  relations  are  used, 

k 

^ '  -/W 

u'  = 

/^r/At 

where  -f=f~  )  's  tne  ratio  of  the  actual  chamber  oscillation 

frequency  to  the  frequency  of  acoustic  oscillations,  and  Ae,/a^  is 
the  nozzle  contraction  ratio,  whicn  is  uniquely  related  to  the  Mach 
numuer  at  the  nozzle  entrance  ffe 

Calculations  have  been  carried  out  for  two  convergent  section 


«4|M8£« 


HANie*. 


/ 

\ 

( 


(C-6) 


fC-7) 


geometries,  which  will  be  referred  to  as  the  •’linear”  and  "conical"  nozzles 
respectively.  The  physicc-  profNes  of  the  nozzles  are  shown  in  Figure  Cl, 


C-4 


and  the  correspo..- i ng  veiocity  prof  ies  in  Figure  C2.,  The  linear  nozzle  is 
constructed  so  that  the  velocity  varies  linearly  with  distance  f  long  the 
axis  in  the  converging  section.  The  resulting  expression  for  ^  , 

=  I  +  X  (C-8) 


where  the  throat  section  corresponas  to  K=-0  and  the  positive  X 
direction  is  the  direction  of  flow,  iii-'-its  in  some  simp!  i  /  icctiwi  in  the 
computations.  The  velocity  profile  i  i  the  conical  nozzic,  is  decidedly  not 
linear,  as  shown  in  Figure  C2.  However,  because  of  the  circular  arc  throat 
section,  the  velocity  variation  can  oe  taken  as  linea'"  'n  tno  vicini+y  of  the 
throat  section,  then  bet”een  the  throat  section  ano  a  stahion  a  short  distance 
away,  the  velocity  can  be  represented  by  an  expression  of  the  form  of 
equation  lC-8) .  Beyond  that  station,  a  more  complicated  expression  is 
necessary.  For  the  30°  hair-angle  convergent  section  geometry  used  in  the 
experimental  phase  of  the  transverse  instability  investigation,  the  follow¬ 
ing  expressions  nave  been  found  to  approximate  well  the  subsonic  velocity 
distribution; 


1  t-x  ^  X  <o 

iZs  (o.4of  X  -  ■f-0.94zQ)^  xs  X.,, 


( 


(o-S) 


where  is  the  station  of  changeovei  from  the  linear  to  the  nonf inear 

vfeiociiy  profile,  and  is  tne  limit  of  the  range  of  interest  of  the 

calculations. 

Tfie  caiculation  process  may  be  thought  of  as  having  three  stao'-s: 
first,  the  solurion  of  the  differential  equation  for  the  auxiliary  function 
,  second,  tne  calculation  of  tne  ofher  functions  ?  '  ,  , 

and  fi'  ,  and  third,  the  computation  or  the  admittance  coefficients. 
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The  fjnction  j  satisfies  tne  crnpiex  Riccaii  equation 

A  ( If  t'  f  0 


(C-IO) 


where 


A  .=  u'  (c»-uM 


o  —  if?  ^  J  ^ 
“  ■7-  dx 


=  to'i’ 


i? 

-  ~  »■-'  .t 


^  si; 


C  _  Kl.*  —  ^  * 

4-  dx 

Equation  (C-iO)  is  equivalent  to  two  real  equat'ions 


(C~ll) 


=  i.  a  <•’  -f-1  )  -  2  •t' '  zT/ 

where  fj'=  C+i  it  .  These  differentia!  equations  are  nonlineer,  and 
have  a  singularity  at  the  origin,  which  corresponds  to  the  nozzle  throat, 
where  u-c-^  I  .  The  boundary  conditions  for  the  functions  and 

are  taken  from  the  p'lysical  condition  that  be  regular  at  the 

sonic  tnrcat  (see  Reference  I,  Appendix  B).  Tnerefore,  the  solution  must 
be  s farted  at  X  =  0.  This  can  be  done  by  expanding  and  in 

power  series  about  the  origin.  In  a  certain  range,  >f ,  c  X  g  o  ,2 
suitable  approximation  is  obtained  by  keeping  only  terms  up  to  the  first 
power  in  x.  The  solution  in  fhis  range  is  found  to  be 

[i-  1'“'] 

+  (u'Csif )  +  ( 'V  -  A*)  ♦  P»  1 


(C-12) 


(C-13) 
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‘  \  ^  ^  _  I  "i  "‘^  '6~r^  . ,  *  /^/  ^  ,0  -  '  •  -j  /^■' '}  Vy  X>  > 

'•'  LV-j)->^->-  j-  "C‘=^o“:^o/  •'■  2  (-f  I  <>o|3c  ? 

rV^-/  ^  /uj''^,sl^\~] 
coLT  "S"  Jj 


to 


where 


^'+'  /  ifl!’  -  i’Ao  ^-'  /  ■■ 


+  uj'^  ''  '" 


then 


where 


-  P  . 

Ot  >c  ''')(!■• 

^ff  ht-i  ~  /}  ^‘*w+(  )  =  (A:  ^  ) 

■4.,  =  Ax  -  j 

=  Ax  -  (xo+^  5  \  ) 

4.3  t.  Ax.  t  {X^+^  i  ) 

^  A  X  •  .r  -i-  Ax  J  4  “!<5  j 

A  x  =  Xvi4.i  -Xi, 

Explici;  integral  expressions  can  De  founc  for  the  auxiliary 


vo<  f  c2Lf  i 


es  ■? 

J 

.  1;  .  4' 

:  for 

A  ' 

2C^  -2  >4 

1  d“ 

J  <i  < 

(C-14) 


f>t‘  \-  ■  ■  ~  ^ 

'  ■^  J 

Beyond  the  station  X,  >  the  solution  can  proceed  without  difficulty 
because  there  are  no  other  singularities  in  +he  range  of  interest.  In 
the  prese-’t  computation,  the  di r  +  orofr^ia!  cquaticns  were  integrated  by 
means  of  tne  method  of  Runye-tfutta,  which  evaluates  the  Dependent  variable 
at  each  x-step  by  means  of  the  derivative  calcula'ied  louolly.  That  is, 
g  i  ven 

j 


(C-15) 


Hov-/e''er,  because  of  tne  oscillatory  nature  of  the  inter'-and,  it  was  Decided 


to  evaluate  tne  auxi I iarv  variaoles  by  airecT  numerical  integration  of  the 


V 
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following  differential  equations: 


£  ]  t-  (s:^06^V  ' 


^(1)1  ei^‘  ,  i^c^-'a  / 
“  \ 


35  a.-u,' 


=  ih)  f  r' 

u. 


(C-i6) 


(C-i7) 


(C-!8) 


The  integration  of  equations  (C-!6)  To  (C-:S)  proceeds  in  tne  same  manner 

as  that  discussed  above  for  equations  (C-M)  and  (C-12).  The  fact  tiiat 

I 

this  procedure  leaves  +he  functions  ^  and  ^  undetermined  at  the 
origin  is  unimportant  in  the  calculation  of  tne  admittance  coe. r ic ients. 

Although  it  is  convenient  comceptoally  to  divide  the  conpulation 
process  into  three  stages,  in  practice  it  is  possible  to  carry  out  the 
calculations  of  a!!  stages  in  one  operation,  at  each  x-station.  A  single 
integration  rill  thus  give  the  coefficients  for  a  wide  range  of  chamber 
exit  Mach  number  values,  but  for  only  one  value  of  the  fequency  and  mode 
parameters  srd  •  The  osci !  lator-  nature  of  the 

functions  involved  in  the  integration  process  make  it  necessary  to  use  a 
small  value  of  the  integ-'tion  step  .  In  order  to  cover  a  range 

of  Mach  numbers  dovrn  +o  0.0”,  with  =0.001,  some  I9C0  integration 

steps  are  required.  For  use  in  comparing  the  theoretical  and  experimental 
results,  especially  in  the  case  of  the  sector  motor,  a  rather  wide  range  is 


required  for  the  parameter  Svj, 


/  / 

,  while  file  reduced  frequpncy  /S’y4 


is  restricted  to  a  fairly  narrow  range.  However,  in  order  to  make  the 
results  more  general,  a  much  wider  frequency  range  was  used.  The  calculations 


UJ  ^  t 

were  made  for  values  of  between  I  and  5,  and  for  Ava  between 


C-8 


0.1  and  :,0.  it  is  interesting  to  compare  the  approxiniate  times  required 

by  various  computational  ir.etho<^j  to  perform  the  calcuiatioris  tor  50  pairs 

ot  values  of  '  and 

desk  calculator  -  iOOO  years 

IBM  650  -  SCO  hours 

IBM  704 - 6  hours 

IBM  7090  - 1-1/2  hours 

Typical  results  of  tne  cci-ui  ions  are  shown  in  rig-ies  C3  through 

CIO,  The  coef  * c •  e^ts  plotred  and  , 

where 

E  =  Cr  -  v''-  /  t-  t  (8  <C-I9) 

since  the  nozzle  adm'i lance  reta+ion  is  used  in  the  form 

for  the  analysis  ot  purely  transverse  modes  with  combustion  complete  at  the 
nozzie  entrance  (see  Chapter  III).  FiguresC3  and  C4  show  £>«  and  £c 
for  a  conical  nozzle  with  a  chamber  exit  Mach  number  of  0.05;  Figures  C5  and  C6 
are  plotted  for  a  Mach  number  of  0,10,  Corresponding  values  of  the  admittance 
coefficients  obtained  for  the  I  inear  nozzle  are  presented  in  Figures  C7  to  C-O. 

The  "range  of  interest"  shown  on  the  graphs  of  the  coel f icio.  ts  is  that  coi — 
responding  to  chamber  oscillation  frequencies  within  +10  percent  of  the 
acoustical  oscillation  frequency.  This  restriction  of  tne  chamber  treouency 
range  has  been  observed  experimentally,  and  is  relate-'  to  the  theoretical 
assumption  that  the  axia'  velocity  perturbation  is  an  order  of  magnitude 
smaller  than  the  pressure  and  radial  velocity  perturbations.  It  can  be  seen 
that  the  admittance  coefficients  for  the  cor'ical  nozzle  'ar/  n^re  rapidly 
with  frequency  than  do  those  for  thg  linear  nozzle.  However,  in  the  rang? 
of  interest  of  this  study,  the  coefficients  for  the  two  nozzles  are  of 
approximately  the  same  magnitude. 
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